A Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations
|
|
- Poppy Cummings
- 5 years ago
- Views:
Transcription
1 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
2 Collaborators Tom Lewis, North Carolina Stefan Schnake, Tennessee The work to be presented here has been partially supported by NSF
3 Outline Motivation and Background A Narrow-Stencil Finite Difference Method High Order Extensions Numerical Experiments Conclusion
4 Outline Motivation and Background A Narrow-Stencil Finite Difference Method High Order Extensions Numerical Experiments Conclusion
5 We consider second order fully nonlinear PDEs F(D 2 u, Du, u, x) = 0 Two best known classes of equations: Monge-Ampére equation: det(d 2 u) = f HJB equations: inf ν V (L ν u f ν ) = 0, where L ν u := A ν (x) : D 2 u + b ν (x) u + c ν (x)u, ν V. Both equations arise from many applications such as differential geometry, optimal mass transfer, stochastic optimal control, mathematical finance etc. Remark: c ν 0 in several applications (e.g., stochastic optimal control, Bellman reformulation of Monge-Ampère equation).
6 Example: (Stochastic Optimal Control) Suppose a stochastic process x(τ) is governed by the stochastic differential equation dx(τ) = f (τ, x(τ), u(τ)) dt + σ (τ, x(τ), u(τ)) dw (τ), τ (t, T ] x(t) = x Ω R n, W : Wiener process u : control vector and let [ ] T J (t, x, u) = E t x L (τ, x(τ), u(τ)) dτ + g (x(t )). t Stochastic optimal control problem involves minimizing J (t, x, u) over all u U for each (t, x) (0, T ] Ω.
7 Bellman Principle Suppose u U such that and define the value function u argmin J (t, x, u), u U v (t, x) = J (t, x, u ). Then, v is the minimal cost achieved starting from the initial value x(t) = x, and u is the optimal control that attains the minimum.
8 Bellman Principle (Continued) Let Ω R n, T > 0, and U R m. The Bellman Principle says v is the solution of for with v t = F(D 2 v, v, v, x, t) in (0, T ] Ω, (1) F(D 2 v, v, v, x, t) = inf (L uv h u ), u U n n n L u v = ai,j u (t, x)v x i x j + bi u (t, x)v xi + c u (t, x) v i=1 j=1 i=1 A u := 1 2 σσt b u := f (t, x, u) c u := 0 h u := L (t, x, u)
9 Ellipticity Definition: Let F[u] := F (D 2 u, u, u, x) and A, B SL(n). (a) F is said to be uniformly elliptic if Λ > λ > 0 such that λtr(a B) F(A, p, r, x) F(B, p, r, x) Λtr(A B) A B. (b) F is said to be proper elliptic if F(A, p, v, x) F(B, p, w, x) A B; v, w R d, v w. (c) F is said to be degenerate elliptic if F(A, p, r, x) F(B, p, r, x) A B.
10 Viscosity solutions Definitions: Assume F is elliptic in a function class A B(Ω) (set of bounded functions), (i) u A is called a viscosity subsolution of F[u] = 0 if ϕ C 2, when u ϕ has a local maximum at x 0 then F (D 2 ϕ(x 0 ), Dϕ(x 0 ), u (x 0 ), x 0 ) 0 (ii) u A is called a viscosity supersolution of F[u] = 0 if ϕ C 2, when u ϕ has a local minimum at x 0 then F (D 2 ϕ(x 0 ), Dϕ(x 0 ), u (x 0 ), x 0 ) 0 (iii) u A is called a viscosity solution of F[u] = 0 if u is both a sub- and supersolution of F[u] = 0 where u (x) := lim sup u(x ) and u (x) := lim inf u(x ) are the x x x x upper and lower semi-continuous envelops of u
11 Barles-Souganidis Framework I For approximating viscosity solutions, we first recall the Barles-Souganidis framework. Theorem (Barles-Souganidis ( 91)) Suppose that the elliptic problem F[u](x) = 0 in Ω satisfies the comparison principle. Assume that the (approximation) operator S : R + Ω R B(Ω) R is consistent, monotone and stable (as well as admissible), then the solution u ρ of problem: S(ρ, x, u ρ (x), u ρ ) = 0 in Ω, converges locally uniformly to the unique viscosity of u.
12 Barles-Souganidis Framework II (i) Admissibility and Stability. For all ρ > 0, there exists a solution u ρ B(Ω) to the following problem: S(ρ, x, u ρ (x), u ρ ) = 0 in Ω. Moreover, there exists a ρ-independent constant C > 0 such that u ρ L (Ω) C. (ii) Monotonicity. For all x Ω, t R and ρ > 0 S(ρ, x, t, u) S(ρ, x, t, v) u, v B(Ω), u v. (iii) Consistency. For all x Ω and φ C (Ω) there hold lim sup ρ 0 y x ξ 0 lim inf ρ 0 y x ξ 0 S(ρ, y, φ(y) + ξ, φ + ξ) ρ S(ρ, y, φ(y) + ξ, φ + ξ) ρ F (D 2 φ(x), φ(x), φ(x), x), F (D 2 φ(x), φ(x), φ(x), x).
13 Wasow-Motzkin Theorem Theorem (Wasow-Motzkin ( 53)) Any monotone and consistent method has to be a wide stencil scheme. Theorem (Bonnans and Zidani ( 03)) If A ν in the HJB equation is not diagonally dominant, then wide stencils are required to preserve monotonicity. Remark: The difficulty is the directional resolution. Wider stencils are used to increase the resolution.
14 Remarks on Monotone Schemes and Wide-stencils Why monotonicity? For the numerical scheme to identify the correct viscosity solution, it needs to respect ordering in some sense. The monotonicity provides such an ordering, which is the best known one so far. A drawback" of monotonicity is that one must use wide stencils according to Wasow-Motzkin ( 53), which could be very problematic for anisotropic problems because very fine directional resolution is required, besides the difficulty from handling boundary conditions.
15 Remarks on Monotone Schemes and Wide-stencils Why monotonicity? For the numerical scheme to identify the correct viscosity solution, it needs to respect ordering in some sense. The monotonicity provides such an ordering, which is the best known one so far. A drawback" of monotonicity is that one must use wide stencils according to Wasow-Motzkin ( 53), which could be very problematic for anisotropic problems because very fine directional resolution is required, besides the difficulty from handling boundary conditions. Consequently, in order to avoid using wide stencils, one must relax (or abandon) the concept of monotonicity (in the sense of Barles and Souganidis).
16 Outline Motivation and Background A Narrow-Stencil Finite Difference Method High Order Extensions Numerical Experiments Conclusion
17 Goals: To construct finite difference methods (FDMs) whose solutions converge to viscosity solutions of the underlying fully nonlinear 2nd order PDE problems, especially, to go beyond the domain of Barles-Souganidis framework and to be more suitable for FDMs and DG methods. Remark: A few existing narrow-stencil" methods are Glowinski et al. ( 04-12): Mixed FE for MA eqns (H 2 solns). Brenner et al. ( 09-13): DG for MA eqns (classical solns). Jensen-Smears ( 12): Linear FE for isotropic HJB eqns. Smears-Süli ( 14): DG-FE for (Cordes-) HJB (H 2 solns). F.-Neilan ( 07-11): FE and DG based on the vanishing moment approach.
18 Finite Difference Operators Let {e j } d j=1 denote the canonical basis of Rd. Define δ + x k,h k v(x) v(x + h ke k ) v(x) h k, δ x k,h k v(x) v(x) v(x h ke k ) h k δ µν x k,h k v(x) δ ν x k,h k δ µ x k,h k v(x), δ µν x k,h k ;x l,h l v(x) δ ν x l,h l δ µ x k,h k v(x) Discrete Gradients: Two natural "sided" choices [ ] ± h k δ± x k,h k Discrete Hessians: Four natural "sided" choices [ D µν h ]k,l δµν x k,h k ;x l,h l, µ, ν {, +} Remark: Low-regularity can be resolved by using sided" gradient and Hessian approximations.
19 Ideas Used for 1st Order Hamilton-Jacobi Equations (Crandall and Lions, 84) FD schemes with the form Ĥ( h U α, + h U α, U α, x α ) = 0 converge to the viscosity solution of a Hamilton-Jacobi equation assuming Consistency: Ĥ(q, q, u, x) = H(q, u, x), Monotonicity: Ĥ(,, u, x). Ĥ is called a numerical Hamiltonian. Ĥ is a function of both h U α and + h U α. The monotonicity requirement is compatible with a discrete first derivative test.
20 Vanishing Viscosity and Numerical Viscosity H( u, u, x) = 0 ɛ u ɛ + H( u ɛ, u ɛ, x) = 0 (E. Tadmor, 97) Every convergent monotone finite difference scheme for HJ equations implicitly approximates the differential equation βh u + H( u, x) = 0 for sufficiently large and possibly nonlinear β > 0, where βh u is called a numerical viscosity. Note: If h k h, then 1 ( d + h U α h U α) = h h U α h δx 2 k,hu α. Lax-Friedrichs numerical Hamiltonian k=1 ( q Ĥ(q, q + + q + ), u, x) H, u, x b (q + q ) 2 }{{} Numerical Viscosity
21 FD Approximations of F(D 2 u, u, u, x) = 0 Since u xi x j may be discontinuous at x α, it needs be approximated from multiple directions.
22 FD Approximations of F(D 2 u, u, u, x) = 0 Since u xi x j may be discontinuous at x α, it needs be approximated from multiple directions. There are 4 possible FD approximations of u xi,x j (x α ): and u xi,x j (x α ) δ µ x i δ ν x j u(x α ) D 2 u(x α ) D µν h u(x α) µ, ν {+, }. µ, ν {+, }.
23 FD Approximations of F (D 2 u, u, u, x) = 0 Inspired by the above observations, we propose the following form of FD methods for F(D 2 u, u, u, x) = 0: F (D ++ h U α, D + h U α, D + h U α, D h U α, h U α, + h U α, U α, x α ) = 0 F is called a numerical operator, which needs to satisfy
24 FD Approximations of F (D 2 u, u, u, x) = 0 Inspired by the above observations, we propose the following form of FD methods for F(D 2 u, u, u, x) = 0: F (D ++ h U α, D + h U α, D + h U α, D h U α, h U α, + h U α, U α, x α ) = 0 F is called a numerical operator, which needs to satisfy Consistency: F(P, P, P, P, q, q, u, x) = F(P, q, u, x), Generalized Monotonicity: F(,,,,,,, x), uses the natural (partial) orderings for symmetric matrices, vectors, and scalars Solvability/Admissibility and Stability: h 0 > 0 and C 0 > 0, which is independent h, such that F [U α, x α ] = 0 has a (unique) solution U and U l (T h ) < C 0 for h < h 0.
25 Remarks on Numerical Operator F F depends on all four "sided" Hessians and both "sided" gradients. The generalized monotonicity requirement is an extension of the Crandall and Lions framework that enforces the elliptic structure of the PDE with respect to the mixed discrete Hessians D ± h. If F is not continuous, F may not be continuous either. In this case the consistency definition should be replaced by lim inf F (P 1, P 2, P 3, P 4, r, s, y) F (P, q, t, x), P k P,r q s t,y x lim sup F (P 1, P 2, P 3, P 4, r, s, y) F (P, q, t, x). P k P,r qn s t,y x Above consistency and monotonicity are different from Barles-Souganidis definitions.
26 Examples of Numerical Operators F I Lax-Friedrichs-like schemes: 1-D (F-Lewis-Kao, 14) p1 + p 2 + p ) 3 F 1 (p 1, p 2, p 3, x) := F(, x + α(p 1 2p 2 + p 3 ) 3 F 2 (p 1, p 2, p 3, x) := F(p 2, x) + α(p 1 2p 2 + p 3 ) p1 + p ) 3 F 3 (p 1, p 2, p 3, x) := F(, x + α(p 1 2p 2 + p 3 ) 2
27 Examples of Numerical Operators F II Godunov-like schemes: 1-D (F-Lewis-Kao, 14): F 4 (p 1, p 2, p 3, x) := ext F(p, x) p I(p 1,p 2,p 3 ) where I(p 1, p 2, p 3 ) := [p 1 p 2 p 3, p 1 p 2 p 3 ] and ext := p I(p 1,p 2,p 3 ) min p I(p 1,p 2,p 3 ) max p I(p 1,p 2,p 3 ) if p 2 max{p 1, p 3 }, if p 2 min{p 1, p 3 }, min if p 1 < p 2 < p 3, p 1 p p 2 min if p 3 < p 2 < p 1. p 3 p p 2
28 Examples of Numerical Operators F III Godunov-like schemes: 1-D (F-Lewis-Kao, 14) (continued): F 5 (p 1, p 2, p 3, x) := extr F(p, x) p I(p 1,p 2,p 3 ) where I(p 1, p 2, p 3 ) := [p 1 p 2 p 3, p 1 p 2 p 3 ] and extr := p I(p 1,p 2,p 3 ) min p I(p 1,p 2,p 3 ) max p I(p 1,p 2,p 3 ) if p 2 max{p 1, p 3 }, if p 2 min{p 1, p 3 }, max if p 1 < p 2 < p 3, p 2 p p 3 max if p 3 < p 2 < p 1. p 2 p p 1
29 Higher Dimensional Lax-Friedrichs-like Schemes Central Discrete Hessian: D 2 hu α 1 4 Central Discrete Gradient: ( D ++ h + D + h + D + h h U α 1 2 ( + h + h ) Uα Lax-Friedrichs-like Numerical Operator: ) F [U α, x α ] F (D 2 hu α, 2 h U α, U α, x α + D ) h Uα + A(U α, x α ) : ( D ++ h U α D + h U α D + h b(u α, x α ) ( + h U α h U ) α U α + D h U ) α
30 Numerical Moment I The discrete operator A(U α, x α ) : ( D ++ h U α D + is called a numerical moment. Observation: h U α D + h U α + D h U ) α ( δ + x k,h k δ + x l,h l δ + x k,h k δ x l,h l δ x k,h k δ + x l,h l + δ x k,h k δ x l,h l ) Uα = h k h l δ 2 x k,h k δ 2 x l,h l U α = h k h l δ 2 x l,h l δ 2 x k,h k U α, an O ( h 2 k + h2 l ) approximation of uxk x k x l x l (x α ) scaled by h k h l. 1 d d : ( D ++ h U α D + h U α D + h U α + D h U α) h 2 2 u(x α ),
31 Numerical Moment II Remark: Here the numerical moment plays a similar role for the 2rd order fully nonlinear PDEs to what the numerical viscosity does for the 1st order fully nonlinear PDEs. The introduction of numerical moments is very much consistent with the vanishing moment method (at the PDE level) proposed and analyzed by F. and Neilan ( 07-12): ε 2 u ε + F(D 2 u ε, u ε, u ε, x) = 0, with the special choice of the parameter ε = αh 2. Is there an analogue of E. Tadmor s result for 2rd order PDEs?
32 Convergence I Assumption: F is uniformly elliptic and Lipschitz continuous in the Hessian argument, and A and b are uniformly bounded for the family of linear operators defining the HJB problem with Dirichlet boundary data. Assume c 0. Theorem (F-Lewis, 16) The Lax-Friedrichs-like scheme is admissible, consistent, stable, and generalized-monotone. Main ingredients of proof: Consistency is trivial. Generalized monotonicity is also easy" by choosing the numerical moment large enough (related to the Lipschitz constant of F).
33 Convergence II Admissibility and stability are the most difficult to verify. The idea is to show that the mapping M ρ : U α Ûα defined by h Û α = h U α ρ F[U α ] is monotone and a contraction for small enough ρ > 0. h is an enhanced discrete Laplacian. Stability follows from applying Crandall-Tartar lemma to M ρ which commutes with additive constants (this is the reason to assume c 0).
34 Convergence III Theorem (F-Lewis, 16) In addition, suppose F[u] = 0 satisfies the comparison principle. Let U be the solution to the Lax-Friedrichs-like scheme and u h be its piecewise constant extension. Then u h converges to u locally uniformly as h 0 +. Main ingredients of proof Follow and modify Barles-Souganidis proof to show that u(x) := lim inf u h (x) and u(x) := lim sup u h (x) are respectively viscosity super- and sub-solution. First work with quadratic test functions, then with general test functions. Use the numerical moment to control the off-diagonal entries in the discrete Hessians (this is the most difficult and technical step). Use the consistency and comparison principle to get u = u.
35 The Local Stencil (2D) The extra nodes in the Cartesian directions smooth the approximation through the diagonal components of the vanishing moment. Directional resolution is no longer necessary avoiding the need for a wide-stencil.
36 Overcoming a Lack of Monotonicity v v x 0 x k x 0 x k δ 2 x k,h k v(x 0 ± h k e k ) < δ 2 x k,h k v(x 0 ) < 0 δ 2 x k,h k v(x 0 ) δ 2 x k,h k v(x 0 ) < 0 δ 2 x k,h k v(x 0 ± h k e k ) > 0 and δ 2 x k,h k v(x 0 ) < 0 δ 2 x k,h k v(x 0 ) δ 2 x k,h k v(x 0 ) > 0
37 Outline Motivation and Background A Narrow-Stencil Finite Difference Method High Order Extensions Numerical Experiments Conclusion
38 Extensions Goals: To develop high order methods and to use unstructured meshes Bad" news: It is not possible to construct higher than 2nd order monotone FD schemes (we have yet given up since g-monotonicity is a weaker requirement)
39 Extensions Goals: To develop high order methods and to use unstructured meshes Bad" news: It is not possible to construct higher than 2nd order monotone FD schemes (we have yet given up since g-monotonicity is a weaker requirement) Good" news: Inspired by a work of Yan-Osher (JCP, 11) for fully nonlinear 1st order Hamilton-Jacobi equations, we are able to construct high order MDG (mixed DG) and LDG (local DG) methods. Ideas of MDG: write F(u xx, x) = 0 as F ( p R, p M, p R, x) = 0 p L = u xx p M = u a xx p R = u + xx and use different numerical fluxes on the linear equations.
40 Extensions (continued) Ideas of LDG: write F(u xx, u x, u, x) = 0 as F( p 1, p 2, p 3, p 4, q 1, q 2, u, x) = 0, q 1 = u x (x ), q 2 = u x (x + ), p 1 = q 1x (x ), p 2 = q 1x (x + ), p 3 = q 2x (x ), p 4 = q 2x (x + ), where u x (x ) (resp. q jx (x )) and u x (x + ) (resp. q jx (x + )) denote the left and right limits of u x (resp. q jx ) at x. They dictate how numerical fluxes should be chosen in the discretization. and stands for monotone increasing and decreasing, respectively. r = 0 is allowed!
41 Outline Motivation and Background A Narrow-Stencil Finite Difference Method High Order Extensions Numerical Experiments Conclusion
42 1-D simulations Test 1: Consider the problem uxx = 0, 0 < x < 1 u(0) = 0, u(1) = 0.5 This problem has the pointwise solutions u + (x) = 0.5x 2 (convex), u (x) = 0.5x 2 + x (concave) Computed using Lax-Friedrichs-like scheme F 1 with α = 1 (left) and α = 1 (right)
43 Test 2: Consider the problem min A θu xx S(x) = 0, 1 < x < 1 θ {1,2} u( 1) = 1, u(1) = 1 { 12x 2, if x < 0 A 1 = 1, A 2 = 2, S(x) = 24x 2, if x 0 This problem has the exact solution u(x) = x x 3
44 Test 3: Consider the problem u 3 xx + 8 sign(x) = 0, 1 < x < 1 u( 1) = 1, u(1) = 1. This problem has the exact solution u(x) = x x C 1 ([ 1, 1]) which is not classical
45 2-D simulations Test 4: Let Ω = (0, 1) 2. Consider Monge-Ampére equation det(d 2 u) = 1 u = 0 in Ω on Ω Remark: No explicit solution formula and solution is not classical.
46 ( α < 0: concave solution) ( α > 0: convex solution)
47 Test 5: Let Ω = ( 1, 1) 2. Consider Monge-Ampére equation det(d 2 u(x, y)) = 0 u = g in Ω on Ω Choose g such that the exact solution is u(x, y) = x.
48 h x L norm order L 2 norm order 2.50E E E E E E E E E Top: α = I, r = 1, and h y = 1/3 fixed. Bottom: α = I, r = 1, and odd number of intervals in the x-direction.
49 3-D Simulations Test 6: Let Ω = (0, 1) 3. Consider Monge-Ampére equation det(d 2 u) = f u = g in Ω on Ω f (x, y, z) = (1 + x 2 + y 2 + z 2 )e 3(x2 +y 2 +z 2 ) 2 g(x, y, z) = e y2 +z 2 2 if x = 0 e x2 +z 2 2 if y = 0 e x2 +y 2 2 if z = 0 e y2 +z if x = 1 e x2 +z if y = 1 e x2 +y if z = 1 The exact solution is u 0 (x, y, z) = e x2 +y 2 +z 2 2
50 (computed solution) (error function)
51 Test 7: Let Ω = (0, 1) 3. Consider Monge-Ampére equation det(d 2 u) = 1 u = 0 in Ω on Ω Remark: No explicit solution formula and solution is not classical.
52 computed solution (α > 0)
53 Outline Motivation and Background A Narrow-Stencil Finite Difference Method High Order Extensions Numerical Experiments Conclusion
54 Concluding Remarks The generalized monotone structure allows us to deal with low regularity functions by considering multiple gradient and Hessian approximations. The key new concept is that of a numerical moment which can be used to design generalized monotone methods such as the Lax-Friedrichs-like scheme. By using narrow stencil schemes that can be expressed using only forward and backward difference quotients, the methods can be easily extended to higher-order and non-cartesian domains/grids using the DG techniques (or DG Finite Element Calculus). Formulation easily extends to Monge-Ampère type equations either directly or based on its HJB reformulation, although the analysis is not yet. Many open problems/issues: such as extension to degenerate PDES, rate of convergence in C 0 -norm, boundary layers, parabolic PDEs, nonlinear solvers, etc.
55 References XF, R. Glowinski and M. Neilan, Recent Developments in Numerical Methods for Fully Nonlinear 2nd Order PDEs", SIAM Review, 55:1-64, 2013 M. Neilan, A. Salgado and W. Zhang, Numerical Analysis of Strongly Nonlinear PDEs", arxiv.org/abs/ , to appear in Acta Numerica. Check arxiv.org for more (and new) references.
56 Thanks for Your Attention!
Fast convergent finite difference solvers for the elliptic Monge-Ampère equation
Fast convergent finite difference solvers for the elliptic Monge-Ampère equation Adam Oberman Simon Fraser University BIRS February 17, 211 Joint work [O.] 28. Convergent scheme in two dim. Explicit solver.
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
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 informationBuilding Solutions to Nonlinear Elliptic and Parabolic Partial Differential Equations
Building Solutions to Nonlinear Elliptic and Parabolic Partial Differential Equations Adam Oberman University of Texas, Austin http://www.math.utexas.edu/~oberman Fields Institute Colloquium January 21,
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 informationThe Kuratowski Ryll-Nardzewski Theorem and semismooth Newton methods for Hamilton Jacobi Bellman equations
The Kuratowski Ryll-Nardzewski Theorem and semismooth Newton methods for Hamilton Jacobi Bellman equations Iain Smears INRIA Paris Linz, November 2016 joint work with Endre Süli, University of Oxford Overview
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 informationLecture 3: Hamilton-Jacobi-Bellman Equations. Distributional Macroeconomics. Benjamin Moll. Part II of ECON Harvard University, Spring
Lecture 3: Hamilton-Jacobi-Bellman Equations Distributional Macroeconomics Part II of ECON 2149 Benjamin Moll Harvard University, Spring 2018 1 Outline 1. Hamilton-Jacobi-Bellman equations in deterministic
More informationCONVERGENT DIFFERENCE SCHEMES FOR DEGENERATE ELLIPTIC AND PARABOLIC EQUATIONS: HAMILTON JACOBI EQUATIONS AND FREE BOUNDARY PROBLEMS
SIAM J. NUMER. ANAL. Vol. 44, No. 2, pp. 879 895 c 26 Society for Industrial and Applied Mathematics CONVERGENT DIFFERENCE SCHEMES FOR DEGENERATE ELLIPTIC AND PARABOLIC EQUATIONS: HAMILTON JACOBI EQUATIONS
More informationNumerical schemes of resolution of stochastic optimal control HJB equation
Numerical schemes of resolution of stochastic optimal control HJB equation Elisabeth Ottenwaelter Journée des doctorants 7 mars 2007 Équipe COMMANDS Elisabeth Ottenwaelter (INRIA et CMAP) Numerical schemes
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
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 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 informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
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 informationOn the approximation of the principal eigenvalue for a class of nonlinear elliptic operators
On the approximation of the principal eigenvalue for a class of nonlinear elliptic operators Fabio Camilli ("Sapienza" Università di Roma) joint work with I.Birindelli ("Sapienza") I.Capuzzo Dolcetta ("Sapienza")
More informationA STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI EQUATIONS *
Journal of Computational Mathematics Vol.8, No.4, 010, 55 568. http://www.global-sci.org/jcm doi:10.408/jcm.1003-m0016 A STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI
More informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
More informationarxiv: v2 [math.na] 20 Jan 2016
Piecewise Constant Policy Approximations to Hamilton-Jacobi-Bellman Equations C. Reisinger and P.A. Forsyth arxiv:1503.05864v2 [math.na] 20 Jan 2016 January 21, 2016 Abstract An advantageous feature of
More informationHJ equations. Reachability analysis. Optimal control problems
HJ equations. Reachability analysis. Optimal control problems Hasnaa Zidani 1 1 ENSTA Paris-Tech & INRIA-Saclay Graz, 8-11 September 2014 H. Zidani (ENSTA & Inria) HJ equations. Reachability analysis -
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 informationMean field games and related models
Mean field games and related models Fabio Camilli SBAI-Dipartimento di Scienze di Base e Applicate per l Ingegneria Facoltà di Ingegneria Civile ed Industriale Email: Camilli@dmmm.uniroma1.it Web page:
More informationHIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS
June 6, 7 :7 WSPC/Lecture Notes Series: 9in x 6in chapter HIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS Chi-Wang Shu Division of Applied Mathematics, Brown University Providence,
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 informationarxiv: v1 [math.ap] 12 Aug 2016
arxiv:1608.03682v1 [math.ap] 12 Aug 2016 Viscosity solutions for junctions: well posedness and stability Pierre-Louis Lions 1 and Panagiotis Souganidis 2,3 October 15, 2018 Abstract We introduce a notion
More informationR. Abgrall 18/9/2008
order HJ on Institut de Mathématiques de Bordeaux et INRIA Bordeaux Sud Ouest 18/9/2008 order 1 order 2 3 order Collaboration with J.D. Benamou, former Otto project, INRIA Roquencourt, now at Rice D. Ribereau
More informationSeismic imaging and optimal transport
Seismic imaging and optimal transport Bjorn Engquist In collaboration with Brittany Froese, Sergey Fomel and Yunan Yang Brenier60, Calculus of Variations and Optimal Transportation, Paris, January 10-13,
More informationAN OVERVIEW OF STATIC HAMILTON-JACOBI EQUATIONS. 1. Introduction
AN OVERVIEW OF STATIC HAMILTON-JACOBI EQUATIONS JAMES C HATELEY Abstract. There is a voluminous amount of literature on Hamilton-Jacobi equations. This paper reviews some of the existence and uniqueness
More informationAn adaptive RBF-based Semi-Lagrangian scheme for HJB equations
An adaptive RBF-based Semi-Lagrangian scheme for HJB equations Roberto Ferretti Università degli Studi Roma Tre Dipartimento di Matematica e Fisica Numerical methods for Hamilton Jacobi equations in optimal
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 informationA VARIATIONAL METHOD FOR THE ANALYSIS OF A MONOTONE SCHEME FOR THE MONGE-AMPÈRE EQUATION 1. INTRODUCTION
A VARIATIONAL METHOD FOR THE ANALYSIS OF A MONOTONE SCHEME FOR THE MONGE-AMPÈRE EQUATION GERARD AWANOU AND LEOPOLD MATAMBA MESSI ABSTRACT. We give a proof of existence of a solution to the discrete problem
More informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationOn Aleksandrov-Bakelman-Pucci Type Estimates For Integro-Differential Equations
On Aleksandrov-Bakelman-Pucci Type Estimates For Integro-Differential Equations Russell Schwab Carnegie Mellon University 2 March 2012 (Nonlocal PDEs, Variational Problems and their Applications, IPAM)
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationA positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations
A positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations Sashank Srinivasan a, Jonathan Poggie a, Xiangxiong Zhang b, a School of Aeronautics and Astronautics,
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More 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 informationConstrained Optimal Stopping Problems
University of Bath SAMBa EPSRC CDT Thesis Formulation Report For the Degree of MRes in Statistical Applied Mathematics Author: Benjamin A. Robinson Supervisor: Alexander M. G. Cox September 9, 016 Abstract
More informationA Nonlinear PDE in Mathematical Finance
A Nonlinear PDE in Mathematical Finance Sergio Polidoro Dipartimento di Matematica, Università di Bologna, Piazza di Porta S. Donato 5, 40127 Bologna (Italy) polidoro@dm.unibo.it Summary. We study a non
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 informationFast Marching Methods for Front Propagation Lecture 1/3
Outline for Front Propagation Lecture 1/3 M. Falcone Dipartimento di Matematica SAPIENZA, Università di Roma Autumn School Introduction to Numerical Methods for Moving Boundaries ENSTA, Paris, November,
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 informationGENERAL EXISTENCE OF SOLUTIONS TO DYNAMIC PROGRAMMING PRINCIPLE. 1. Introduction
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
More informationGradient Estimate of Mean Curvature Equations and Hessian Equations with Neumann Boundary Condition
of Mean Curvature Equations and Hessian Equations with Neumann Boundary Condition Xinan Ma NUS, Dec. 11, 2014 Four Kinds of Equations Laplace s equation: u = f(x); mean curvature equation: div( Du ) =
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
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 informationarxiv: v2 [math.na] 2 Aug 2013
A VISCOSITY SOLUTION APPROACH TO THE MONGE-AMPÈRE FORMULATION OF THE OPTIMAL TRANSPORTATION PROBLEM arxiv:1208.4873v2 [math.na] 2 Aug 2013 JEAN-DAVID BENAMOU, BRITTANY D. FROESE, AND ADAM M. OBERMAN Abstract.
More informationExistence and uniqueness results for the MFG system
Existence and uniqueness results for the MFG system P. Cardaliaguet Paris-Dauphine Based on Lasry-Lions (2006) and Lions lectures at Collège de France IMA Special Short Course Mean Field Games", November
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition Sukjung Hwang CMAC, Yonsei University Collaboration with M. Dindos and M. Mitrea The 1st Meeting of
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationDuality and dynamics in Hamilton-Jacobi theory for fully convex problems of control
Duality and dynamics in Hamilton-Jacobi theory for fully convex problems of control RTyrrell Rockafellar and Peter R Wolenski Abstract This paper describes some recent results in Hamilton- Jacobi theory
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 Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
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 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 informationHeterogeneous Agent Models in Continuous Time Part I
Heterogeneous Agent Models in Continuous Time Part I Benjamin Moll Princeton University of Chicago, 28 September 2015 What this mini-course is about Many interesting questions require thinking about distributions
More informationStochastic Homogenization for Reaction-Diffusion Equations
Stochastic Homogenization for Reaction-Diffusion Equations Jessica Lin McGill University Joint Work with Andrej Zlatoš June 18, 2018 Motivation: Forest Fires ç ç ç ç ç ç ç ç ç ç Motivation: Forest Fires
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
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 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 informationOptimal transport for Seismic Imaging
Optimal transport for Seismic Imaging Bjorn Engquist In collaboration with Brittany Froese and Yunan Yang ICERM Workshop - Recent Advances in Seismic Modeling and Inversion: From Analysis to Applications,
More informationA Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang 3, and Hongkai Zhao 4 ABSTRACT In this paper, we construct a second order fast
More informationHopfLax-Type Formula for u t +H(u, Du)=0*
journal of differential equations 126, 4861 (1996) article no. 0043 HopfLax-Type Formula for u t +H(u, Du)=0* E. N. Barron, - R. Jensen, and W. Liu 9 Department of Mathematical Sciences, Loyola University,
More informationarxiv: v1 [math.ap] 10 Apr 2013
QUASI-STATIC EVOLUTION AND CONGESTED CROWD TRANSPORT DAMON ALEXANDER, INWON KIM, AND YAO YAO arxiv:1304.3072v1 [math.ap] 10 Apr 2013 Abstract. We consider the relationship between Hele-Shaw evolution with
More informationThe Fast Sweeping Method. Hongkai Zhao. Department of Mathematics University of California, Irvine
The Fast Sweeping Method Hongkai Zhao Department of Mathematics University of California, Irvine Highlights of fast sweeping mehtod Simple (A Gauss-Seidel type of iterative method). Optimal complexity.
More informationJean-David Benamou 1, Brittany D. Froese 2 and Adam M. Oberman 2
ESAIM: MA 44 (010) 737 758 DOI: 10.1051/man/010017 ESAIM: Mathematical Modelling and umerical Analysis www.esaim-man.org TWO UMERICAL METHODS FOR THE ELLIPTIC MOGE-AMPÈRE EQUATIO Jean-David Benamou 1,
More informationHomogenization of a Hele-Shaw-type problem in periodic time-dependent med
Homogenization of a Hele-Shaw-type problem in periodic time-dependent media University of Tokyo npozar@ms.u-tokyo.ac.jp KIAS, Seoul, November 30, 2012 Hele-Shaw problem Model of the pressure-driven }{{}
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 informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
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 informationConvergence Analysis of Fast Sweeping Method for Static Convex Hamilton-Jacobi Equations
Convergence Analysis of Fast Sweeping Method for Static Convex Hamilton-Jacobi Equations Songting Luo 1, Department of Mathematics, Iowa State University, Ames IA 50011. Hongkai Zhao 2 Department of Mathematics,
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 informationA stochastic particle system for the Burgers equation.
A stochastic particle system for the Burgers equation. Alexei Novikov Department of Mathematics Penn State University with Gautam Iyer (Carnegie Mellon) supported by NSF Burgers equation t u t + u x u
More informationHigh order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes
High order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes Feng Zheng, Chi-Wang Shu and Jianxian Qiu 3 Abstract In this paper, a new type of high order Hermite
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationOn continuous time contract theory
Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem
More informationMulti-dimensional Stochastic Singular Control Via Dynkin Game and Dirichlet Form
Multi-dimensional Stochastic Singular Control Via Dynkin Game and Dirichlet Form Yipeng Yang * Under the supervision of Dr. Michael Taksar Department of Mathematics University of Missouri-Columbia Oct
More informationSOLUTION OF STEADY-STATE HAMILTON-JACOBI EQUATION BASED ON ALTERNATING EVOLUTION METHOD
ITALIAN JOURNAL OF PURE AND APPLIED MATHEMATICS N. 38 2017 (184 193 184 SOLUTION OF STEADY-STATE HAMILTON-JACOBI EQUATION BASED ON ALTERNATING EVOLUTION METHOD A. Tongxia Li Department of Mathematics Hulunbeier
More informationA Monotone Scheme for Hamilton-Jacobi Equations via the Nonstandard Finite Difference Method
A Monotone Scheme for Hamilton-Jacobi Equations via the Nonstandard Finite Difference Method Roumen Anguelov, Jean M-S Lubuma and Froduald Minani Department of Mathematics and Applied Mathematics University
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 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 informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
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 informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS
ALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS HAILIANG LIU AND MICHAEL POLLACK Abstract. In this work, we propose a high resolution Alternating Evolution Discontinuous
More informationGERARD AWANOU. i 1< <i k
ITERATIVE METHODS FOR k-hessian EQUATIONS GERARD AWANOU Abstract. On a domain of the n-dimensional Euclidean space, and for an integer k = 1,..., n, the k-hessian equations are fully nonlinear elliptic
More informationarxiv: v2 [math.ap] 28 Nov 2016
ONE-DIMENSIONAL SAIONARY MEAN-FIELD GAMES WIH LOCAL COUPLING DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES arxiv:1611.8161v [math.ap] 8 Nov 16 Abstract. A standard assumption in mean-field game
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 informationREGULARITY RESULTS FOR THE EQUATION u 11 u 22 = Introduction
REGULARITY RESULTS FOR THE EQUATION u 11 u 22 = 1 CONNOR MOONEY AND OVIDIU SAVIN Abstract. We study the equation u 11 u 22 = 1 in R 2. Our results include an interior C 2 estimate, classical solvability
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
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 informationHigh-order filtered schemes for time-dependent second order HJB equations
High-order filtered schemes for time-dependent second order HJB equations Olivier Bokanowski, Athena Picarelli, and Christoph Reisinger November 15, 2016 Abstract In this paper, we present and analyse
More informationExam February h
Master 2 Mathématiques et Applications PUF Ho Chi Minh Ville 2009/10 Viscosity solutions, HJ Equations and Control O.Ley (INSA de Rennes) Exam February 2010 3h Written-by-hands documents are allowed. Printed
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationNumerical schemes for short wave long wave interaction equations
Numerical schemes for short wave long wave interaction equations Paulo Amorim Mário Figueira CMAF - Université de Lisbonne LJLL - Séminaire Fluides Compréssibles, 29 novembre 21 Paulo Amorim (CMAF - U.
More information