R. Abgrall 18/9/2008
|
|
- Dale Tate
- 6 years ago
- Views:
Transcription
1 order HJ on Institut de Mathématiques de Bordeaux et INRIA Bordeaux Sud Ouest 18/9/2008
2 order 1 order 2 3
3 order Collaboration with J.D. Benamou, former Otto project, INRIA Roquencourt, now at Rice D. Ribereau and S. Ballereau, SNPE St Médard en Jalles S. Augoula former PhD student Financial support from SNPE, CNES
4 order order
5 Motivations order Not done before ( 93 94). Work of C.W. Shu and S.J. Osher on Cartesian, J.D. Benamou (rapport d activité INRIA Projet Ondes, 93?) Multivalued solutions for Eikonal equation, Computation of multivalued solutions, Big Ray method (Benamou) Computation for Ariane V boosters : what is the flight time lenght of a booster (SNPE CNES ALENIA) General domain, accuracy.
6 Cauchy problem for HJ in R d order u + H(x,u, u) = 0 t u(x,t = 0) = u 0 (x,t) Hamiltonian H uniformly continuous u 0 uniformly continuous Seek for continuous solutions
7 order Link with hyperbolic problems (2D) Differenciate u t set (p,q) T = (u x,u y ) T t ( ) ( p Hp 0 + q 0 H p + H(x,u, u) = 0 )( ) ( p Hq 0 + q 0 H x q )( ) p = 0 q y Problem :Weakly hyperbolic system (one eigenvalue) Problem : Which entropy condition(s)?
8 Non uniqueness order Solutions u(x,t) = u 0 (x) etc u t + (u ) 2 2 { = 0 x R x if x > 0 u 0 (x) = x else { u0 (x) if x > t u(x,t) = t x2 2 else
9 Viscosity solutions, 1 order Let ε > 0, u ε solution of ( uε ) t + H(x,u ε, u ε ) = ε u ε,u ε (x,0) = u 0 (x). Let ε 0. Take φ C 2 (R R + ) such that u ε ϕ maximal at some point x 0,t 0. Then u ε (t 0,x 0 ) = φ(x 0,t 0 ) supplementary assumption, free u ε t (x 0,t 0 ) = φ t (x 0,t 0 ) u ( ε (x 0,t 0 ) = φ(x 0,t 0 ) u ε φ (x 0,t 0 ) 0
10 Viscosity solutions, 1 order Hence ( φ ) t + H(x,φ, φ) 0 at (x 0,t 0 ).
11 Viscosity solution, 3 order Similarly, take φ such that u ε φ is minimal at some (x 0,t 0 ) Then ( φ ) t + H(x,φ, φ) 0 at (x 0,t 0 ).
12 u Sub solution order For any φ C 2 (R d R + ). If u φ has a local maximum at (x 0,t 0 ) then ( φ ) t + H(x,φ, φ) 0 at (x 0,t 0 ).
13 u Super solution order For any φ C 2 (R d R + ). If u φ has a local minimum at (x 0,t 0 ) then ( φ ) t + H(x,φ, φ) 0 at (x 0,t 0 ).
14 u solution order if u is a sub- and super solution
15 order Standard results Exsitence and uniqueness results if Assumptions H(x, u, p) H(x, v, p) γ R (u v) γ R > 0, R v u R, p H(x, u, p) H(y, u, p) C x y (1 + p ) Assumption H convex there exists a strict sub solution H(x,u,p) + when p + uniformly in x,u + H1
16 Uniqueness principle order Under any of the previous assumptions sub solution super solution
17 Numerical approximation order u t + H( u) = 0 u(x,0) = u 0 (x) {u 0 i } vertices u n+1 i = u n i th(u n i, {u n j } j Vi )
18 Constraints order Consistancy If ui 0 = C + p 0 x i, then H(ui n, {un j } j V i ) = H(p 0 ) Monotonicity If i, ui n vi n then i, u n+1 i v n+1 i under possible CFL, Accuracy Convergence
19 Numerical approximation, Idea 1 order Take the hyperbolic system ( ) ( p Hp 0 + t q 0 H p )( ) ( p Hq 0 + q 0 H x q )( ) p = 0 q y discretise (finite difference/finite volume) Unknown (u x,u y ) T Integrate (u x,u y ) T u i Possible for Cartesian, difficult for
20 idea 1 order
21 Idea 2 : Hopf formulas order define Ψ convex Ψ (p) = sup x R d ( ) x p Ψ(x) define Ψ concave Ψ (p) = inf x R d ( ) x p + Ψ(x)
22 H uniformly Lipschitz order u 0 convex, u 0 concave u(x,t) = sup p R d u(x,t) = inf p R d (x p u 0 (p) th(p) ) ( ) x p + u0(p) th(p)
23 u 0 uniformly continuous order H concave u(x,t) = inf y R d ( ( )) x y u 0 (y) + th t H convex u(x,t) = sup y R d ( ( )) y x u 0 (y) + th t
24 First formula (H Lispchitz) order Osher Bardi, Siam J. Numer Anal., 1989 Decompose u 0 = u convex + u concave u 0 u convex x p + ( u concave ) (p) = v p (x) Solve the Cauchy problem u t + H( u) = 0, Comparison principle u(x,t) sol of above problem p := V p (x,t) u(x,t) inf p V p (x,t) IC v p (x)
25 Osher Bardi order If u 0 defined on vertices of a Cartesian mesh, then uniqueness of decomposition, V p (x,t) easily computable at vertices
26 Case of Unstructured order u 0 defined at vertices (piecewise linear interpolation) u 0 = u convex + u concave Uniqueness only for cartesian Leads to unconsistant Hamiltonians
27 Idea 2 : Huygens principle order Assumption H = H 1 + H 2, H 1 convex in p, H 2 concave in p Dp 2 H bounded,h 2 = α p 2,H 1 = H + α p 2. H 2 (p) p p + (H 2 ) (p ) p H H 1 + p p + (H 2 ) (p ) p Comparison principle?
28 Comparison principle order If H(p) = H 1 (p) + H 2 (p), H 1 convex in p, H 2 concave in p. H, H 1, H 2 uniformly Lipschitz continuous Cauchy problem u t + H( u) = 0 u(x,0) = u 0 (x) Solution u(x,t) Then Φ 2 (x,t) u(x,t) Φ 1 (x,t)
29 Comparison principle order with Φ 1 (x,t) = inf sup q R d y R d [ ( u 0 (y) + th1(q) + th2 q + y x )] t Φ 2 (x,t) = inf sup q R d y R d [ ( u 0 (y) + th2(q) + th1 q + x y )] t
30 Idea of proof order H(p) H 1 (p) p p + H 2 (p ) = G p (p) G p is convex in p u t + H( u) u t + G p ( u) u sub-solution of u t + G p ( u) = 0 u(x,t) V p (x,t) take the inf.
31 Application order {ui 0 } vertex values Construct piecewise linear interpolation u 0 use one of the formula to approximate u 1 i repeat
32 Application order i i Ω l u Ω 0 l (M) = ui 0 + U i M i M
33 Application order Φ(M i, t) = ui 0 t min max set H i = min max q R 2 0 l k i q R 2 0 l k i ) sup (U l (z q) H 1 (z) H 2 (q) z Ω l +q ( ) sup U l (z q) H1(z) H2(q) z Ω l +q Consistant, Monotone (L 1 + L 2 ) t h 1 2
34 Case of H convex smooth order ( ) H i = min inf U i z H (z) 0 l k l z Ω l H ( x(p) ) + H ( p(x) ) = x(p) p(x) p(x) = x H ( x(p) ) x(p) = p H ( p(x) ) ) 1 x(u l ) Ω l, sup x Ωl (U i z H (z) = H(U i ) 2 x(u l ) Ω l, sup reached on boundary of ( Ω l )
35 Eikonal order ( H i = max max 0, u k u i, u ) j u i 0 l k i ki ji k j i
36 Properties order Consistant, Monotone Intrinsic k j i
37 Result, Crandall Lions extended order Let H : R d R uniformly continuous, u 0 bounded uniformly continuous, lipschitz constant L. Consider T triangulation, h = max{radius}, uniform. Numerical Hamiltonian H (1) consistant, (2)monotone, (3) intrinsic Initial condition u 0, ui 0 = u 0 (M i ), u n+1 i = ui n th i. u the viscosity solution of u t + H( u) = 0 u(x,0) = u 0 (x) then for all T > 0 and n, n t T max i un i u(m i,n t) C t
38 Other examples : Lax Friedrichs order H i =H ε ( ) u(m) ui dm ρ C ρ =H ε u ndm ρ C ρ with ǫ L π ( udx ) D ρ H = H πρ 2 or D ρ H( u)dx πρ 2
39 Properties order Consistant, Monotone under CFL Intrinsic by construction k j i
40 Numerical examples order distance functions various geometry
41 dist1 order
42 mesh order
43 dist2 order
44 lens order
45 order
46 order Overview What is the problem for boundary? Viscosity solutions and boundary Numerical approximation
47 The problem order Take Solution u 1 = εu xx x [0,1], u(0) = 0 u(1) = 2 > 1 u ε (x) = x + exp ( ) ( ) x 1 ε exp 1 ε 1 exp ( 1 ) ε
48 order
49 More serious situation sol.di.desc order
50 condition order Define H(x,u(x), u(x)) = 0 x Ω F(x,u(x), u(x) = 0 x Ω { G H(x,t,p) if x Ω (x,t,p) = max(f(x,t,p),h(x,t,p)) if x Ω { H(x,t,p) if x Ω G (x,t,p) = min(f(x,t,p),h(x,t,p)) if x Ω
51 Viscosity sub-solution order u super semi continuous function (semi continue supérieurement) φ C 2 (Ω), if x 0 Ω is a local maximum of u φ, then G (x 0,u ( x 0 ), φ(x 0 )) 0.
52 Viscosity sub-solution order u lower semi continuous function φ C 2 (Ω), if x 0 Ω is a local minimum of u φ, then G (x 0,u ( x 0 ), φ(x 0 )) 0.
53 Dirichlet order G (x,t,p) = G (x,t,p) = { H(x,t,p) if x Ω max(h(x,t,p),u(x) g(x)) if x Ω { H(x,t,p) if x Ω min(h(x,t,p),u(x) g(x)) if x Ω
54 Numerical approximation order H(x,u(x), u(x)) = 0 x Ω F(x,u(x), u(x) = 0 x Ω F(x,t,p) = t g(x) (Dirichlet) or F(x,t,p) = p n g(x) Neuman
55 Numerics order S H (x i,u(x i ), {u j ;j V i }) = 0 x i Ω S(ρ,x i,u(x i ),u) = max(s Hb (x i,u(x i ), {u j ;j V i }, S F (x i,u(x i ), {u j ;j V i }) = 0 x i Ω Under which, does the scheme converge to the solution of the problem.
56 General convergence result order Barles and Souganidis result for S(ρ,x,t,u) Monotonicity: if u v, for all ρ 0, x Ω, t R and u,v L (Ω), then S(ρ,x,t,u) S(ρ,x,t,v). Stability: for all ρ > 0, there exists a solution u ρ L (Ω) to (1) with a bound independant of ρ, Consistency: for all x Ω and φ Cb (Ω), lim sup S(ρ,y,φ(y) + ξ,φ + ξ) G (x,φ(x),dφ(x)) ρ 0,y x,ξ 0 lim inf S(ρ,y,φ(y) + ξ,φ + ξ) G (x,φ(x),dφ(x)) ρ 0,y x,ξ 0
57 General convergence result order and Strong uniqueness principle: if u L (Ω) is a subsolution and v L (Ω) is a supersolution, then u v on Ω.
58 Barles and Souganidis order Theorem (Barles, Souganidis) Assuming the monotonicity, the consistency, the stability and the strong uniqueness properties, then, when ρ 0, the solution u ρ of converges locally uniformly to the unique continuous viscosity solution of H(x,u(x), u) = 0.
59 Application order Assume that 1 H b H, S H, S Hb and S F are monotone and stable, 2 for all φ Cb (Ω), we have, For any x Ω or in a neighborhodd of Ω lim S H(ρ, y, ϕ(y) + ξ, ϕ + ξ) = H(x, ϕ(x), Dϕ(x)) ρ 0,y x,ξ 0 lim S H b (ρ, y, ϕ(y) + ξ, ϕ + ξ) = H b (x, ϕ(x), Dϕ(x)) ρ 0,y x,ξ 0 lim S F(ρ, y, ϕ(y) + ξ, ϕ + ξ) = F(x, ϕ(x), Dϕ(x)) ρ 0,y x,ξ 0 3 Problem has a uniqueness principle,
60 Application order Then the family u ρ defined by { 0 = S(ρ,x,u ρ (x),u ρ SH (ρ,x,u ) = ρ (x),u ρ ) max(s Hb (ρ,x,u ρ (x),u ρ ),S F (ρ,x,u ρ (x) converges locally uniformly to the solution of H(x,u, u) = 0 in Ω F(x,u, u) = 0 in Ω
61 Visual proof order Convex Hamiltonians + Dirichlet + Dirichlet H(x,u,p) = sup { b(x,v).p + λu f (x,v)} v V 0 = u(x) inf v(.) [ min(t,τ) 0 f (y x (t), v(t))e λt dt + 1 {T<τ } u(y x (T))e λt + 1 {T τ } ϕ(y x (τ))e λτ ]
62 Examples order H b = impose strongly boundary condition, Case of convex hamiltonians ( H b = max max Ui y H (y) ) angular sectors y ang. sect i.e. do nothing, etc
63 Numerical examples order u = 1 in Ω u = 0 for x = 0.5 u = C for x = 1. The constant C takes the values displayed in Table 1 Case C Table: condition for the Dirichlet b.c.s
64 Examples order (a) (b) Figure: Dirichlet, case A, (a) weak, (b) comparaison weak ( ) strong ( )
65 order (a) (b) Figure: Dirichlet, case B, (a) weak, (b) comparaison weak ( ) strong ( )
66 order (a) (b) Figure: Dirichlet, case C, (a) weak, (b) comparaison weak ( ) strong ( )
67 order (a) (b) Figure: Dirichlet, case D, (a) weak, (b) comparaison weak ( ) strong ( )
68 Improved convergence result order For convex Hamiltonians, take the Godunov boundary Hamitonian. No boundary layers : control of u.
69 Soner Bcs sol.soner.desc order
70 Dirichlet Bcs order sol.di.desc
71 order Non convex Hamiltonians Situation much less clear...
72 order
73 Scheme : 2 points of view order solve unsteady problem : u t + H(x,u, u) = 0 u(x,0) = u 0 (x) du dt + H(x i,u(x i ), {u j } j Vi ) and Runge Kutta type integration Solve steady problem directly H(x,u, u) = 0 unsteady one : u n+1 u n t u n+1 u n + H(x,u n+1, u n+1 ) = 0 t ( ) H(x,u n+1, u n+1 ) + H(x,u n, u n ) =
74 Overview order 1 ENO type scheme (WENO) 2 Tuning of artificial dissipation 3 Schemes for steady problems
75 order ENO WENO type du dt + H(x i,u(x i ), {u j } j Vi ) Runge Kutta (Shu Osher for preservation of stability properties of H) Take first order scheme H i = H(x i,u i, {u j } j Vi ) H(x i,u i, u T1,, u Tk ) i i Ω l Provide high order accuracy approximations of u T1,, u Tk
76 order Algorithm need to define a family of stencils, compute the one that minimise oscillation
77 approximation of u T order start from node values (Lagrange) k-th order approximation of u T : polynomial of degree k + 1, need (k + 1) ( (k + 1) + 1 ) /2 points construction of stencils with (k + 1) ( (k + 1) + 1 ) /2 need M i inside (k + 1) ( (k + 1) + 1 ) /2 1 additional points
78 Construction of potential stencils order All directions needed, As few stencils as possible, What is oscillation
79 Hierarchical algorithm order Step 1 Consider all triangles that surround M i (about 6 for regular ): u T Choose T min that minimize u T Step 2 : T min known construct 3 more points
80 Hierarchical algorithm, 2 order Centered
81 Facts and justifications order Compute stencils : not always possible (think of 3 aligned points) In practice always a solution P P k interpolation for f smooth max f (r) P (r) h p+1 x K CM p+1 ρ r f not smooth in the stencil of P P = r a ij (x x 0 ) i (y y 0 ) j l=0 i+j=r i+j=r ] a ij C [f (r) h p r
82 Numerical results order H(p,x,y) = p n(x,y) { 1 n(x,y) = ( ) 1 cos 2πd if d 2 = ( ) x 2 ( x 2 0.2) 1 1 else
83 Problems order Non compact stencil, not obvious to increase order (even for WENO) Not obvious for 3D, computational not very efficient
84 order Construction of Compact, problems Discontinuous Galerkin like (see Shu et al., JSC, SIAM SISC, etc), How to have provable convergent (no Lax Wendroff like theorem )? i.e How to have automatic entropy Do this directly in the scheme, i.e. start from a viscosity solution framework. Different from Shu et al.
85 Solution 1 order S. Augoula s thesis. Available at abgrall/these.tar.gz J. Scientific Computing, 2000, vol 2
86 order Solution 2, steady problems, SIAM SISC, in revision Compute H i H(x,u(x), u) as H i = l ( H i ) L + (1 l) ( Hi ) H ( H i ) (L) first order, monotone, approximation of H(x, u(x), u) ( H i ) H high order approximation of H(x,u(x), u)
87 Properties of l, 1 order Steady problems H(x, u(x), u(x)) = 0 u smooth ( H i ) H = O(h p ), p 2 l(u) = O(h p ) then H i = O(h p )
88 Properties of l, 2 order Iterative scheme u n+1 i = u n i th i ( H i ) (1) monotone under t Ch ( H i ) (H) L p stable for t C h l [0,1] same stability range for blended scheme
89 Compute l order l ( ) ) L ( ) H (Hi ) L H i + (1 l) Hi = (l + (1 l)r so ask r = H(H) H (L) l + (1 l)r 0
90 order ell 1 range 1 r
91 accuracy order l ( ( ) ) L ( ) H l (Hi H i + (1 l) Hi = r + (1 l) ) (H) r = H(H) H (L) If l [0,1] and l r C for some C > 0, then l( ) L ( ) H H i + (1 l) Hi C H (H). and high order accuracy.
92 order Solution Hyperbola l + (1 l)r = 0 slope 1 at (0,0), C 1 or l = { 0 if r > 0 min(1,c r ) if r 0 l = min(1,c r ), C 1
93 Minmod like order ell 1 range 1 r
94 order Proof for convergence for H(x, u(x), u) = 0 Barles and Souganidis result for S(ρ,x,t,u) H ρ (x,u(x), {u j } j Vi ) = 0 (1) 1 Monotonicity: if u v, for all ρ 0, x Ω, t R and u, v L (Ω), then S(ρ, x, t, u) S(ρ, x, t, v). 2 Stability: for all ρ > 0, there exists a solution u ρ L (Ω) to (1) with a bound independant of ρ, 3 Consistency: for all x Ω and φ C b (Ω), lim S(ρ, y, φ(y) + ξ, φ + ξ) = H(x, φ(x), Dφ(x)) ρ 0,y x,ξ 0 4 Strong uniqueness principle: if u L (Ω) is a subsolution and v L (Ω) is a supersolution, then u v on Ω.
95 Barles and Souganidis order Theorem (Barles, Souganidis) Assuming the monotonicity, the consistency, the stability and the strong uniqueness properties, then, when ρ 0, the solution u ρ of converges locally uniformly to the unique continuous viscosity solution of H(x,u(x), u) = 0.
96 order Sketch, 1 Assume H (L) monotone l as before assume iterative convergence reached
97 Sketch, 2 order l [0,1] L : bound on the solution (CFL condition) Take l = min(1, r ), same CFL number as low order scheme.
98 Sketch, 3 order u ϕ maximum at x 0 : u(x 0 ) ϕ(x 0 ) ( ) 0 = lh (L) + (1 l)h (H) = l + (1 l)r (u)h (L) (u) ( ) l + (1 l)r (u) H (L) (ϕ) so H (L) (ϕ) 0 if ( l + (1 l)r ) (u) 0 i.e. H (H) (u) 0 then u = limu h super solution of H(x,u(x), u) = 0
99 order Example u 1 = 0 x Ω u = 0 x Γ Ω
100 Examples order first example : Godunov scheme H (H) i = max i u ang sec (x i ) K (k) i = angular sectors K(k) max(0, u p u i ) else pi ( ) if u ang sec (x i ) angular sector
101 order Godunov :Practical implementation H=-1.e+30! Initialisation of the numerical Hamiltonian DO jt=1,nt! loop over element is(:)=nu(:,jt)! index of the 3 vertices of the element jt compute gradient of u compute H T H(is(:)) = max(h(is(:)),h T ) ENDOO
102 Lax Friedrich order H LF ρ ( u Ω1,, u Ωki ) = D h H( u) πh 2 ǫ h C h [u(m) u(m i )]dl.
103 Implementation order H=0! Initialisation of the numerical Hamiltonian DO jt=1,nt! loop over element is(:)=nu(:,jt)! index of the 3 vertices of the element jt compute gradient of u compute the angles θ 1 jt,θ2 jt,θ3 jt compute H T s and D T s H(is(1)) = H(is(1)) + H 1 T + D1 T H(is(2)) = H(is(2)) + H 2 T + D2 T H(is(3)) = H(is(3)) + H 3 T + D3 T
104 Lax Friedrichs 2 order H LF ( u Ω1,, u Ωki ) = T H( u T ) + α (u i u j ) T M i M j T T T and α h T max p p H where h T is the largest edge of T.
105 order Implementation H=0! Initialisation of the numerical Hamiltonian air=0! Initialisation of the area DO jt=1,nt! loop over element is(:)=nu(:,jt)! index of the 3 vertices of the element jt compute gradient of u compute the angles θ 1 jt,θ2 jt,θ3 jt compute H T s and D T s ENDDO H(is(1)) = H(is(1)) + T (H T + DT 1 ) H(is(2)) = H(is(2)) + T (H T + DT 2 ) H(is(3)) = H(is(3)) + T (H T + DT 3 ) air(is(1)) = air(is(1)) + T air(is(2)) = air(is(2)) + T air(is(3)) = air(is(3)) + T
106 order any of the above with P 2 interpolation i p
107 or SUPG like order H σ (σ,u σ,u ξξ Vσ ) := + h Ω ( H(x, u h )N σ dx Ω ] [ p H(x, u h ) p H(x, u h ) N σ / N σ dx Ω H(x, u h )dx )
108 One dimensional examples order u n(x) = 0 x [0,1] u(0) = u(1) = 0 where n(x) = 3x 2 + a with a = 1 2.x3 0 2x 0 1 and x 0 = The solution is u(x) = { x 3 + ax if x [0,x 0 ] 1 + a ax x 3 if x [x 0,1]
109 numerics order First order Blended scheme Second order Exact Comparison of the first order, second order and blended, 50 mesh points, zoom in [0.4,0.8]
110 ... order First order Second order, blended Second order, unlimited Exact solution 0.75 First order Second order, blended Second order, unlimited Exact solution
111 Errors order First order Second order unlimited Second order limited Slope 1 Slope 2 5 First order Second order unlimited Second order limited Slope 2 slope L 1 norm L norm
112 Non homogeneous Hamiltonian order The problem reads H(x,p) = cos(p) 2 + p. H(x,u ) = 0 in x [0,1],with u(0) = u(1) =
113 order Numerics numericaly checked that the gradient of the solution may be larger that 1 in absolute value, so that non convex effect are real First order Second order limited Second order unlimited
114 Distance function order
115 numerics order First order solution (blue isolines), second order unlimited solution (gree isolines) and second order limited solution (red isolines).
116 Non homogeneous case order Ω = [0,1] [0,1] C where C = {(x,y) [0,1/2] [0,1/2] and (x 1/2) 2 + (y 1/2) 2 1}, ψ( u) 1 = 0 in Ω and u = 0 on Ω. p ψ(p) = 2 p if p 1 else.
117 numerics order
118 order numerics First order Second order
119 numerics order
120 Distance order
121 order First order solution
122 order Second order solution, no limitation
123 Second order, no limitation order
124 Second order, no limitation order
125 Second order solution, limitation order
126 3D case order
127 second order, no limitation order
128 second order, limitation order
129 Error order
130 Error order
HIGH 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 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 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 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 informationSolution Sheet 3. Solution Consider. with the metric. We also define a subset. and thus for any x, y X 0
Solution Sheet Throughout this sheet denotes a domain of R n with sufficiently smooth boundary. 1. Let 1 p
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 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 informationRECENT DEVELOPMENTS IN VERY HIGH ORDER RESIDUAL DISTRIBUTION SCHEMES
RECENT DEVELOPMENTS IN VERY HIGH ORDER RESIDUAL DISTRIBUTION SCHEMES FOR INVISCID AND VISCOUS PROBLEMS. R. Abgrall Team Bacchus INRIA Bordeaux Sud Ouest and UniversitédeBordeauxTalence,France Roscoff,
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 informationFinite difference Hermite WENO schemes for the Hamilton-Jacobi equations 1
Finite difference Hermite WENO schemes for the Hamilton-Jacobi equations Feng Zheng, Chi-Wang Shu 3 and Jianian Qiu 4 Abstract In this paper, a new type of finite difference Hermite weighted essentially
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 informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
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 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 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 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 informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
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 informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
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 central discontinuous Galerkin method for Hamilton-Jacobi equations
A central discontinuous Galerkin method for Hamilton-Jacobi equations Fengyan Li and Sergey Yakovlev Abstract In this paper, a central discontinuous Galerkin method is proposed to solve for the viscosity
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 informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
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 informationIMA Preprint Series # 1993
AN ADAPTIVE METHOD WITH RIGOROUS ERROR CONTROL FOR THE HAMILTON-JACOBI EQUATIONS. PART I: THE ONE-DIMENSIONAL STEADY STATE CASE By Bernardo Cockburn and Bayram Yenikaya IMA Preprint Series # 993 ( September
More informationSemi-discrete central-upwind schemes with reduced dissipation for Hamilton Jacobi equations
IMA Journal of Numerical Analysis 005 5, 113 138 doi: 10.1093/imanum/drh015 Semi-discrete central-upwind schemes with reduced dissipation for Hamilton Jacobi equations STEVE BRYSON Programme in Scientific
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 informationFast 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 informationA Central Discontinuous Galerkin Method for Hamilton-Jacobi Equations
J Sci Comput (00) 45: 404 48 DOI 0.007/s095-009-9340-y A Central Discontinuous Galerkin Method for Hamilton-Jacobi Equations Fengyan Li Sergey Yakovlev Received: 6 August 009 / Revised: November 009 /
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 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 informationDiscontinuous Galerkin and Finite Difference Methods for the Acoustic Equations with Smooth Coefficients. Mario Bencomo TRIP Review Meeting 2013
About Me Mario Bencomo Currently 2 nd year graduate student in CAAM department at Rice University. B.S. in Physics and Applied Mathematics (Dec. 2010). Undergraduate University: University of Texas at
More informationAn adaptive mesh redistribution method for nonlinear Hamilton-Jacobi equations in two- and three-dimensions
An adaptive mesh redistribution method for nonlinear Hamilton-Jacobi equations in two- and three-dimensions H.-Z. Tang, Tao Tang and Pingwen Zhang School of Mathematical Sciences, Peking University, Beijing
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 informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationInterpolation. 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter
Key References: Interpolation 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter 6. 2. Press, W. et. al. Numerical Recipes in C, Cambridge: Cambridge University Press. Chapter 3
More informationSemi-Discrete Central-Upwind Schemes with Reduced Dissipation for Hamilton-Jacobi Equations
Semi-Discrete Central-Upwind Schemes with Reduced Dissipation for Hamilton-Jacobi Equations Steve Bryson, Aleander Kurganov, Doron Levy and Guergana Petrova Abstract We introduce a new family of Godunov-type
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 informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
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 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 informationMath Partial Differential Equations 1
Math 9 - Partial Differential Equations Homework 5 and Answers. The one-dimensional shallow water equations are h t + (hv) x, v t + ( v + h) x, or equivalently for classical solutions, h t + (hv) x, (hv)
More informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationResidual Distribution. basics, recents developments, relations with other techniques
Residual Distribution basics, recents developments, relations with other techniques MARIO RICCHIUTO July 19, 2012 INTRODUCTION with historical perspective t u+ F(u) = 0 In the 80 s numerical techniques
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
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 informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationConvergence of a large time-step scheme for Mean Curvature Motion
Outline Convergence of a large time-step scheme for Mean Curvature Motion M. Falcone Dipartimento di Matematica SAPIENZA, Università di Roma in collaboration with E. Carlini (SAPIENZA) and R. Ferretti
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationPropagation d interfaces avec termes non locaux
Propagation d interfaces avec termes non locaux P. Cardaliaguet Univ. Brest Janvier 2008 Joint works with G. Barles (Tours), O. Alvarez (Rouen), O. Ley (Tours), R. Monneau (CERMICS), A. Monteillet (Brest).
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 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 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 informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
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 informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
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 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 informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationAn Introduction to Variational Inequalities
An Introduction to Variational Inequalities Stefan Rosenberger Supervisor: Prof. DI. Dr. techn. Karl Kunisch Institut für Mathematik und wissenschaftliches Rechnen Universität Graz January 26, 2012 Stefan
More informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationMaster Thesis. Nguyen Tien Thinh. Homogenization and Viscosity solution
Master Thesis Nguyen Tien Thinh Homogenization and Viscosity solution Advisor: Guy Barles Defense: Friday June 21 th, 2013 ii Preface Firstly, I am grateful to Prof. Guy Barles for helping me studying
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 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 informationDevelopment and stability analysis of the inverse Lax-Wendroff boundary. treatment for central compact schemes 1
Development and stability analysis of the inverse Lax-Wendroff boundary treatment for central compact schemes François Vilar 2 and Chi-Wang Shu 3 Division of Applied Mathematics, Brown University, Providence,
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
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 Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationHybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation
Numerical Analysis and Scientific Computing Preprint Seria Hybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation W. Guo J. Qiu Preprint #21 Department of Mathematics University
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
More informationA POSTERIORI ERROR ESTIMATES FOR TIME-DEPENDENT HAMILTON-JACOBI EQUATIONS
A POSTERIORI ERROR ESTIMATES FOR TIME-DEPENDENT HAMILTON-JACOBI EQUATIONS A DISSERTATION SUBMITTED TO THE FACULTY OF THE GRADUATE SCHOOL OF THE UNIVERSITY OF MINNESOTA BY IVAN GEORGIEV MEREV IN PARTIAL
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 Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations
Int. Conference on Boundary and Interior Layers BAIL 6 G. Lube, G. Rapin Eds c University of Göttingen, Germany, 6 A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations. Introduction
More informationAccepted Manuscript. A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
Accepted Manuscript A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang, Hongkai Zhao PII: S-()- DOI:./j.jcp... Reference: YJCPH To
More informationDiscontinuous Galerkin methods Lecture 2
y y RMMC 2008 Discontinuous Galerkin methods Lecture 2 1 Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu y 1 0.75 0.5 0.25 0-0.25-0.5-0.75 y 0.75-0.0028-0.0072-0.0117 0.5-0.0162-0.0207-0.0252
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 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 informationNon-linear Methods for Scalar Equations
Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 3, 04 / 56 Abstract
More informationLois de conservations scalaires hyperboliques stochastiques : existence, unicité et approximation numérique de la solution entropique
Lois de conservations scalaires hyperboliques stochastiques : existence, unicité et approximation numérique de la solution entropique Université Aix-Marseille Travail en collaboration avec C.Bauzet, V.Castel
More informationVISCOSITY SOLUTIONS OF HAMILTON JACOBI EQUATIONS WITH DISCONTINUOUS COEFFICIENTS
Dept. of Math. Univ. of Oslo Pure Mathematics No. 2 ISSN 86 2439 January 25 VISCOSITY SOLUTIONS OF HAMILTON JACOBI EQUATIONS WITH DISCONTINUOUS COEFFICIENTS GIUSEPPE MARIA COCLITE AND NILS HENRIK RISEBRO
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 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 informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
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 information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More information