LECTURE 3: DISCRETE GRADIENT FLOWS
|
|
- Donna Bishop
- 5 years ago
- Views:
Transcription
1 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 Related Topics Isaac Newton Institute, January 2014
2 Outline Subdifferentials and Coercivity Subgradient Operators Angle-Bounded Operators Space Discretization Numerical Experiments for TV Flow The Allen-Cahn Equation
3 OUTLINE Subdifferentials and Coercivity Subgradient Operators Angle-Bounded Operators Space Discretization Numerical Experiments for TV Flow The Allen-Cahn Equation
4 Subdifferentials Space: H separable Hilbert space. Functional: φ : H R proper lower semicontinuous (l.s.c) convex function with (nonempty, convex) domain D(φ). Subdifferential: F = φ is the subdifferential of φ: w F(w) w, v w φ(v) φ(w) v D(φ) with D(F) densily contained in D(φ). The operator F is in general multivalued. Monotonocity: F(v) F(w), v w 0 v, w D(F).
5 Energy Solutions Strong Solution: Given u 0 D(F) and f L 1 (0, T ; H), a function u C 0 ([0, T ]; H) is a strong solution of u (t) + F(u(t)) = f(t) a.e. t (0, T ) if it is absolutely continuous in (0, T ) and satisfies the equation. Energy Solutions: Multiply the equation by u(t) v to obtain the variational inequality u (t) f(t), u(t) v + φ(u(t)) φ(v) 0 v D(φ). Energy Identity: If f L 2 (0, T ; H) then t φ(u(t)) is absolutely continuous in (0, T ) and u (t) 2 + d dt φ(u(t)) = f(t), u (t) a.e. t (0, T ). Dissipation Inequality: If u is a strong solution and f is absolutely continuous 1 d 2 dt u (t) 2 f (t), u (t).
6 Example 1: Parabolic Obstacle Problem Space: H = L 2 (Ω) K = {v H 1 0 (Ω) : v(x) ψ(x) a.e. x Ω}. Functional: If I K is the indicator function of K, then φ(w) = 1 Z w 2 + I K(w), D(φ) = K 2 Ω Variational Inequality: u solves Z Z tu(u v) + u (u v) f(u v) v K. Ω Ω
7 Example 2: Degenerate Parabolic PDE Space: H = H 1 (Ω) with w = Gw L 2 (Ω), G = ( ) 1 with homogeneous Dirichlet data. Functional: Let β : R R be continuous and monotone, and Z Z w(x) φ(w) = β(r)dr dx, D(φ) = {v L 2 (Ω) : j(v) L 1 (Ω)}. Ω 0 {z } =j(w(x)) Subdifferential: F(w) = β(w), D(F) = {v D(φ) : β(v) H 1 0 (Ω)} Equation: u solves t β(u) = f.
8 Example 3: Total Variation Flow φ(v) = R v D(φ) := Ω L2 (Ω) BV (Ω). u t div = 0 F(u) = φ(u) = div u u u u Characterization of the subdifferential φ [Andreu et al. 2004]: Let X = z L (Ω) : divz L 2 (Ω), z n Ω = 0 B 1(X) unit ball of X Z divz φ(u) z B 1(X) : z u = φ(u) Ω We exploit this characterization to restate the TV flow: given u D(φ) find u H 1 (0, T ; L 2 (Ω) L (0, T ; D(φ)) and z L (0, T ; X) Initial condition: u t=0 = u 0 TV flow: u t, w + R z w = 0, w D(φ) Ω R divz φ(u): z B 1(X) (q z) u 0 Ω q B1(X).
9 Coercivity σ(w; v) := φ(v) φ(w) F(w), v w 0 ρ(w, v) := σ(w, v) + σ(v; w) = F(w) F(v), w v Variational Inequality: u f, u v + φ(u) φ(v) + σ(u; v) 0 v D(φ). If ρ is p-coercive, p 2, ρ(w; v) [w v] p w, v D(F), then σ(w, v) 1 [w v]p p w D(F), v D(φ). We will ignore coercivity in the rest of the discussion.
10 OUTLINE Subdifferentials and Coercivity Subgradient Operators Angle-Bounded Operators Space Discretization Numerical Experiments for TV Flow The Allen-Cahn Equation
11 Hilbert Theory: Subgradient Operators Energy Solutions: (with vanishing forcing term) u t + F(u) = 0 in H (Hilbert space) F = φ (φ convex) = u (t) 2 + d φ(u(t)) = 0. dt u, u v + φ(u) φ(v) v D(φ) Backward Euler Method: 1 U n U n 1, U n V + φ(u n) φ(v ) V D(φ) Un = E n := Un φ(un) φ(un 1) 0 (Numerical Dissipation)
12 A Posteriori Error Analysis (Steps 1 and 2) We follow [N, Savaré, Verdi 00]: 1. ODE for U(t) = t t n 1 U n + tn t U n 1 (linear interpolant of U n 1 and U n) U, U v + φ(u) φ(v ) U, U U n + φ(u) φ(un) {z } R(t)=residual 2. Estimate of R(t) U U n = t tn `Un U n 1 φ(u) tn t φ(u n 1) + t tn 1 φ(u n) (φ convex) = R(t) = (t n t)e n
13 A Posteriori Error Analysis (Steps 1 and 2) We follow [N, Savaré, Verdi 00]: 1. ODE for U(t) = t t n 1 U n + tn t U n 1 (linear interpolant of U n 1 and U n) U, U v + φ(u) φ(v ) U, U U n + φ(u) φ(un) {z } R(t)=residual 2. Estimate of R(t) U U n = t tn `Un U n 1 φ(u) tn t φ(u n 1) + t tn 1 φ(u n) (φ convex) = R(t) = (t n t)e n
14 A Posteriori Error Analysis (Steps 1 and 2) We follow [N, Savaré, Verdi 00]: 1. ODE for U(t) = t t n 1 U n + tn t U n 1 (linear interpolant of U n 1 and U n) U, U v + φ(u) φ(v ) U, U U n + φ(u) φ(un) {z } R(t)=residual 2. Estimate of R(t) U U n = t tn `Un U n 1 φ(u) tn t φ(u n 1) + t tn 1 φ(u n) (φ convex) = R(t) = (t n t)e n
15 A Posteriori Error Analysis (Steps 3 and 4) 3. Error Equation v = U V = u u, u U + φ(u) φ(u) 0 U, U u + φ(u) φ(u) R(t) = 1 d 2 dt u U 2 R(t) = (t n t)e n 4. Integration in Time max u 0 t T U 2 u 0 U NX τne 2 n n=1
16 A Posteriori Error Analysis (Steps 3 and 4) 3. Error Equation v = U V = u u, u U + φ(u) φ(u) 0 U, U u + φ(u) φ(u) R(t) = 1 d 2 dt u U 2 R(t) = (t n t)e n 4. Integration in Time max u 0 t T U 2 u 0 U NX τne 2 n n=1
17 A Posteriori Error Analysis (Steps 3 and 4) 3. Error Equation v = U V = u u, u U + φ(u) φ(u) 0 U, U u + φ(u) φ(u) R(t) = 1 d 2 dt u U 2 R(t) = (t n t)e n 4. Integration in Time max u 0 t T U 2 u 0 U NX τne 2 n n=1
18 A Posteriori Error Analysis (Steps 3 and 4) 3. Error Equation v = U V = u u, u U + φ(u) φ(u) 0 U, U u + φ(u) φ(u) R(t) = 1 d 2 dt u U 2 R(t) = (t n t)e n 4. Integration in Time max u 0 t T U 2 u 0 U NX τne 2 n n=1
19 A Posteriori Error Analysis via Residual Monotone Terms Continuous Dissipation Inequality: If u is a strong solution and f = 0, then Discrete Dissipation Inequality: 1 τ 2 n Residual Monotone Term: because implies E n = τ 1 n d dt u (t) 2 0. U n U n U τn 1 2 n 1 U n n N E n τn 1 F(U n) F(U n 1), U n U n 1 = D n τ 1 n φ(u n) φ(u n 1) F(U n 1), U n U n 1 φ(u n 1) φ(u n) τn U 2 n U n 1 2 F(U n 1), U n U n 1 + τ 1 F(U n), U n U n 1. n
20 A Priori Rate of Convergence We follow [N, Savaré, Verdi 00]. Similar results in [Baiocchi 89, Rulla 96]. Weak Regularity: 0 φ(u 0) C NX n=1 τ 2 ne n φ(u 0) max 1 n N τn {z } =τ U 0 = u 0 D(φ) max 0 t T u U p τφ(u 0) Strong Regularity: F(U 0) C NX τne 2 n 1 2 F(U0) 2 max 2 1 n N n=1 {z } =τ 2 U 0 = u 0 D(F) max 0 t T u U τ 2 F(U 0)
21 A Priori Rate of Convergence: Proofs Weak Regularity: Since 0 E n τn `φ(un 1) 1 φ(u n), and φ can be assumed to be nonnegative after an affine modification, we deduce U 0 = u 0 NX NX τne 2 n τ φ(u n 1) φ(u n) τφ(u 0) n=1 with τ := max 1 n N, whence max u 0 t T U 2 n=1 NX τne 2 n φ(u 0) τ. n=1 Strong Regularity: Since τn 1 `Un U n 1 = F(Un) we have 0 D n = F(U n) F(U n 1), F(U n) 1 F(U n 1) 2 F(U n) 2 2 whence for U 0 = u 0 max u 0 t T U 2 NX τnd 2 n 1 2 F(u0) 2 τ 2. n=1
22 Examples of Subgradient Operators 1. Example 1: Parabolic Obstacle Problems (H = L 2 (Ω), D(φ) = H 1 0 (Ω), D(F) = H 2 (Ω)) u ψ : u t, u v + u, (u v) 0 v ψ Z n := U n Un Un 1 E n 1 (U n U n 1) Z n Z n 1, U n U n 1 2. Example 2: Degenerate Parabolic Problems (H = H 1 (Ω), D(φ) = L 2 (Ω), D(F) = H 1 0 (Ω)) u t β(u) = 0 (β non-decreasing) E n 1 β(u n) β(u n 1), U n U n 1
23 Example 3: Time-Discrete TV Flow Continuous Problem: Z u t, w + z w = 0, w D(φ) Ω Z z B 1(X) (q z) u 0 q B 1(X) Ω Semidiscrete Variational Scheme: [Bartels, N, Salgado 13] Given τ > 0 find {U n } h 1 (L 2 (Ω)) and {Z n } l (X) such that Z 1 τ U n+1 U n, w + Z n+1 w = 0, w D(φ) Ω Z `q Z n+1 U n+1 0, q B 1(X). Ω A Priori Error Estimates (H = L 2 (Ω)) u 0 D(φ) u U 2 L (H) u 0 U 0 2 H + τφ(u 0 ). u 0 D( φ) u U 2 L (H) u 0 U 0 2 H + τ 2 φ(u 0 ) 0 2 H.
24 OUTLINE Subdifferentials and Coercivity Subgradient Operators Angle-Bounded Operators Space Discretization Numerical Experiments for TV Flow The Allen-Cahn Equation
25 Hilbert Theory: Angle-Bounded Operators (N, Savaré, Verdi 00) u + F(u) = 0, U + F(U n) = 0 1 d 2 dt u U 2 = F(u) F(U n), u U = tn 1 t F(u) F(U n), u U n + tn t F(U n) F(u), u U n 1 γ > 0 : F(v) F(w), w z γ F(v) F(z), v z
26 Examples of Angle-Bounded Operator Example 1: Linear Advection-Diffusion PDE F(v) = ν 2 v + b v + cv, H = L 2(Ω) ( b(x) b 0 d(x) = 1 2 divb(x) + c(x) d2 0 γ 2 = 1 4 b ν 2 (d ν2 p 2 0 ) u U 2 L (L 2 ), ν 2 (u U) 2 L 2 (L 2 ) u 0 U 0 2 L 2 X N Z + 2γ 2 ν 2 (U n U n d(x) U n U n 1 2 n=1 τ 2 n Ω Example 2: Obstacle Problem for Advection-Diffusion PDE F(v) = ν 2 v + b v + cv, v ψ, H = L 2(Ω)
27 OUTLINE Subdifferentials and Coercivity Subgradient Operators Angle-Bounded Operators Space Discretization Numerical Experiments for TV Flow The Allen-Cahn Equation
28 Space Discretization: Abstract Theory (Bartels, N, Salgado 13) Let {H h } h>0 be a family of finite dimensional subspaces of H and φ h : H h R convex and lower semicontinuous functionals such that: Monotonicity: For every v h H h φ h (v h ) φ(v h ). Approximation: For every h > 0, there is an operator C h : W H h with W D(φ) dense such that C h w w H ε 1(h) w W, w W D(φ). Energy diminishing: For all w W φ h (C h w) φ(w) ε 2(h) w W.
29 A Priori Error Estimate Semidiscrete problem Fully discrete problem 1 τ U n+1 U n, U n+1 w + φ(u n+1 ) φ(w) 0. 1 n+1 Uh Uh n, U n+1 h W h + φ h (U k+1 h ) φ h (W h ) ρ n+1, U n+1 h W h, τ with ρ n+1 a convenient perturbation (inexact solution). Theorem (abstract a priori error estimate) If U 0 h W uniformly in h and ρ l (H) Cτ, then U U h 2 l (H) U 0 U 0 h 2 H + C`ε 1(h) + ε 2(h) U l (W) + Cτ.
30 Example: Fully Discrete TV Flow Direct discretization: find {U n h } H h H = L 2 (Ω) such that 1 n+1 Uh Uh n, U n+1 h W h + φ h (U n+1 h ) φ h (W h ) ρ n+1, u n+1 h W h. τ This is equivalent to: find {Uh n } H h and {Z n h} X such that Z 1 n+1 Uh Uh n, W h + Z n+1 h W h = ρ k+1, W h W h V h, τ Ω Z `q Z n+1 h U n+1 h 0 q B 1(X). Ω Questions: The dual variable Z n h is not discrete. How to discretize it? If we discretize it, how do we solve the variational inequality?
31 Fully Discrete TV Flow: Quadrature Finite element space: Let H h be constructed using P 1 or Q 1 continuous finite elements. Scheme: The TV flow with quadrature is to find {U n h } H h such that 1 n+1 Uh Uh n, U n+1 h W h +φ h (U n+1 h ) φ h (W h ) ρ k+1, U k+1 h W h, τ Quadrature: the quadrature rule is φ h (v h ) = Q h ( v h ) = X QX ω q T v h (ξ q) φ(v h ). T T h q=1
32 Fully Discrete TV Flow: Approximation Properties Monotonicity: If Q h is supported at the vertices of the cells, then Z φ h (w h ) = Q h ( w h ) w h = φ(w h ), Approximation: If C h w = Π h w ε is the Clement interpolant of w ε, a C approximation of function w L (Ω) BV (Ω), we get h 2 C h w w L 2 (Ω) C Ω 1 ε + ε 2 {z } =ε 1 (h) w 1 2 L (Ω) Dw (Ω) 1 2. Energy diminishing: If Q h is exact for linears, then φ h (w) φ(w) = Q h ( C h w ) Dw (Ω) h 1 C ε + ε 2 Dw (Ω) {z } =ε 2 (h) w BV (Ω).
33 Fully Discrete TV Flow: A Priori Error Analysis Theorem (a priori error estimates) Let u 0 L (Ω) BV (Ω) and U 0 h = Π h u 0. Then U U h l (L 2 (Ω)) Ch 1 6. Assume we can construct a TV-diminishing interpolation operator C h, i.e. Z Z C h w w, Ω Ω The error becomes O(h 1 4 ). This suboptimal rate is related to lack of time regularity; compare with heat equation. The error for total variation minimization (image processing) becomes O(h 1 2 ). This is optimal given the regularity L (Ω) BV (Ω) B 1 2 (L 2 (Ω)).
34 Construction of the TV Diminishing Operator (Bartels, N, Salgado 13) Let Ω = S 1 S 1 be a periodic domain. The construction extends to non-periodic domains and dimension d > 2. V h p.w. Q 1, continuous. For a vertex z N h W z := 1 h xh y Z Q w v Let Π h : L 1 (Ω) V h be Π h w(z) := W z hy hx
35 TV Diminishing Operator: Stability Properties Theorem (stability) The operator Π h is TV diminishing Z Π h w Dw (Ω) Ω and is L p stable for 1 p, i.e. Πw L p (Ω) w L p (Ω), w BV (Ω) w L p (Ω) The key is to exploit the symmetry of the construction Notice that the constant is also = 1 in the L p estimate
36 TV Diminishing Operator: Approximation Properties Theorem (error estimates) Let h = max{h x, h y}. If w BV (Ω) w Πw L 1 (Ω) Ch Dw (Ω) If w BV (Ω) L p (Ω) w Πw L q (Ω) Ch 1 s Dw (Ω) 1 s w s L p (Ω), 1 q = 1 s + s 1 p If w W 2 1 (Ω) w Πw L 1 (Ω) Ch 2 w W 2 1 (Ω) The operator preserves linears locally = second order
37 OUTLINE Subdifferentials and Coercivity Subgradient Operators Angle-Bounded Operators Space Discretization Numerical Experiments for TV Flow The Allen-Cahn Equation
38 Experiment 1: Characteristic of a Ball If u 0 = χ BR, then [Andreu et. al. 2004] u(x, t) = (1 λt) + χ BR (x), λ = 2 R. The computational convergence rate is O(h 1 2 ) with τ = τ = h 10. No oscillations. No diffusion of the support B R (no regularization) Excellent agreement with the extinction time (not covered by theory).
39 Experiment 2: Effect of Global Forcing (Condensation) We consider f = constant, Neumann BC, u 0 = χ BR and u t div = f. u u characteristic of ball B R Forcing f prevents χ BR to decrease: no diffusion; Characteristic χ Ω/BR grows with constant speed until it merges with χ BR, and next χ Ω grows with the same speed.
40 Experiment 3: Characteristic of a Disconnected Set Let u 0 = P 3 i=1 χ B(x i,r i ), where the centers lie on an equilateral triangle of length 1 and R i 0.2. Then [Andreu et.al. 2004] 3X u(x, t) = λ it) i=1(1 + χ B(xi,R i )(x), λ i = 2 R i
41 Experiment 4: Characteristic of an Annulus Let u 0 = kχ E with E = B R \ B r and k R. Then [Andreu et al. 04] u(x, t) = ( sgn(k) ( k λ Et) + χ E(x) + λ rtχ Br, t < T 1, sgn(κ) ( κ λ r(t T 1)) + χ BR, t > T 1 with λ E = 2t R r, λr = 2 r T 1 = k λ E + λ r, κ = T 1λ r
42 Experiment 4: Characteristic of an Annulus No diffusion of the support. Although not covered by the theory, there is good agreement for the meeting time. Although not covered by the theory, there is good agreement for the extinction time.
43 Experiment 5: Comparison with Regularized TV Flow [Feng-Prohl 03]: If u ε is the solution of the regularized TV flow u ε tu ε div p = 0, ε2 + u ε 2 then u u ε L (L 2 ) Cε 1 2. If u ε L 2 (H 2 ) ε 1, our analysis yields provided ε h dt. u U ε,h L (L 2 ) Ch 1 2
44 Experiment 5: Comparison with Regularized TV Flow 0.5 ε=0 ε=h ε=h u x The mesh size is h = 2 5, the time-step t = 2 10, and T=5. Notice that the requirement t = O(h 2 ) needed for the L 2 -convergence of the regularized flow (Feng-Prohl) is valid.
45 OUTLINE Subdifferentials and Coercivity Subgradient Operators Angle-Bounded Operators Space Discretization Numerical Experiments for TV Flow The Allen-Cahn Equation
46 The Allen-Cahn Equation as a Gradient Flow Space: H = L 2 (Ω) with norm u, v = ε R Ω uv Functional: Z φ(v) = Ω ε 2 u ` u ε Gradient Flow in L 2 : If f(u) = u` u 2 1, then ε tu ε u + 1 ε f(u) = 0 Error Analysis: Since φ is a Lipschitz perturbation of a convex functional, the a priori error analysis applies but with a stability constant C = e T /ε2. This exponential deterioration in the small parameter ε is in general unavoidable as corresponds to the logistic ODE ε t + 1 f(u) = 0. ε Developed Interfaces: Once diffuse interfaces form, the exponential behavior can be replaced a priori by a polynomial dependence on ε 1 [Feng-Prohl 02].
47 The Allen-Cahn Equation: A Posteriori Error Analysis (Kessler-N, Bartels) Singularly Perturbed Parabolic PDE: ε tu, v + ε u, v + 1 ε f(u), v = 0, f(u) = u`u 2 1 Finite Element Discretization: U n V n : for all V V n ε U n I nu n 1, V + ε U n, V + 1 ε f(inun 1), V = 0 Piecewise Linear Interpolant: Residual: U(t) = t tn 1 U n + tn t U n 1 R(t), v := ε tu, v + ε U, v + 1 f(u), v ε
48 A Posteriori Error Estimate Without Exponential Spectral Estimate [DeMottoni-Schatzman, X. Chen] ε v 2 L 2 (Ω) + 1 ε f (u)v, v λ 0ε v 2 L 2 (Ω) v H 1 (Ω) E 0, E 1 computable estimators in terms of the residual R(t) Tolerance: δ Cε 4 Conditions on E 0, E 1: λ T, λ 0 T, λ 1 T e 2λ 0T u 0 U 0 L2 (Ω) λ T δ, E 0 λ 0 T εδ, E 1 λ 1 T ε 2 δ A posteriori error estimate (with polynomial dependence on ε): u U L2 (0,T ;L 2 (Ω)) δ
A Posteriori Error Estimates for Variable Time-Step Discretizations of Nonlinear Evolution Equations
A Posteriori Error Estimates for Variable Time-Step Discretizations of Nonlinear Evolution Equations RICARDO H. NOCHETTO University of Maryland GIUSEPPE SAVARÉ Universitá di Pavia AND CLAUDIO VERDI Università
More informationOn some nonlinear parabolic equation involving variable exponents
On some nonlinear parabolic equation involving variable exponents Goro Akagi (Kobe University, Japan) Based on a joint work with Giulio Schimperna (Pavia Univ., Italy) Workshop DIMO-2013 Diffuse Interface
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
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 informationNonlinear Evolution Governed by Accretive Operators in Banach Spaces: Error Control and Applications
Nonlinear Evolution Governed by Accretive Operators in Banach Spaces: Error Control and Applications Ricardo H. Nochetto Department of Mathematics and Institute for Physical Science and Technology, University
More informationNumerical Solution of Nonsmooth Problems and Application to Damage Evolution and Optimal Insulation
Numerical Solution of Nonsmooth Problems and Application to Damage Evolution and Optimal Insulation Sören Bartels University of Freiburg, Germany Joint work with G. Buttazzo (U Pisa), L. Diening (U Bielefeld),
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
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 informationContents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping.
Minimization Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping. 1 Minimization A Topological Result. Let S be a topological
More informationA POSTERIORI ERROR ANALYSIS FOR HIGHER ORDER DISSIPATIVE METHODS FOR EVOLUTION PROBLEMS
A POSTERIORI ERROR ANALYSIS FOR HIGHER ORDER DISSIPATIVE METHODS FOR EVOLUTION PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. We prove a posteriori error estimates for time discretizations
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationTD 1: Hilbert Spaces and Applications
Université Paris-Dauphine Functional Analysis and PDEs Master MMD-MA 2017/2018 Generalities TD 1: Hilbert Spaces and Applications Exercise 1 (Generalized Parallelogram law). Let (H,, ) be a Hilbert space.
More informationEntropy-dissipation methods I: Fokker-Planck equations
1 Entropy-dissipation methods I: Fokker-Planck equations Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Introduction Boltzmann equation Fokker-Planck equations Degenerate parabolic
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
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 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 informationNUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE XIAOBING FENG AND ANDREAS PROHL Abstract. In this
More informationStability and Geometric Conservation Laws
and and Geometric Conservation Laws Dipartimento di Matematica, Università di Pavia http://www-dimat.unipv.it/boffi Advanced Computational Methods for FSI May 3-7, 2006. Ibiza, Spain and Introduction to
More informationCalculus of Variations. Final Examination
Université Paris-Saclay M AMS and Optimization January 18th, 018 Calculus of Variations Final Examination Duration : 3h ; all kind of paper documents (notes, books...) are authorized. The total score of
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationGALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS
GALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND CHARALAMBOS MAKRIDAKIS Abstract. We consider discontinuous as well as continuous Galerkin methods for the time discretization
More informationUna aproximación no local de un modelo para la formación de pilas de arena
Cabo de Gata-2007 p. 1/2 Una aproximación no local de un modelo para la formación de pilas de arena F. Andreu, J.M. Mazón, J. Rossi and J. Toledo Cabo de Gata-2007 p. 2/2 OUTLINE The sandpile model of
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationDaniel Kessler 1, Ricardo H. Nochetto 1, 2 and Alfred Schmidt 3
ESAIM: MAN MAN, Vol. 38, N o 1, 004, pp. 19 14 DOI: 10.1051/man:004006 ESAIM: Mathematical Modelling and Numerical Analysis A POSTERIORI ERROR CONTROL FOR THE ALLEN CAHN PROBLEM: CIRCUMVENTING GRONWALL
More informationarxiv: v1 [math.ap] 31 May 2007
ARMA manuscript No. (will be inserted by the editor) arxiv:75.4531v1 [math.ap] 31 May 27 Attractors for gradient flows of non convex functionals and applications Riccarda Rossi, Antonio Segatti, Ulisse
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
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 informationA POSTERIORI ERROR ESTIMATES FOR THE BDF2 METHOD FOR PARABOLIC EQUATIONS
A POSTERIORI ERROR ESTIMATES FOR THE BDF METHOD FOR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND PANAGIOTIS CHATZIPANTELIDIS Abstract. We derive optimal order, residual-based a posteriori error estimates for
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 informationKrein-Rutman Theorem and the Principal Eigenvalue
Chapter 1 Krein-Rutman Theorem and the Principal Eigenvalue The Krein-Rutman theorem plays a very important role in nonlinear partial differential equations, as it provides the abstract basis for the proof
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
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 informationt y n (s) ds. t y(s) ds, x(t) = x(0) +
1 Appendix Definition (Closed Linear Operator) (1) The graph G(T ) of a linear operator T on the domain D(T ) X into Y is the set (x, T x) : x D(T )} in the product space X Y. Then T is closed if its graph
More informationLocal strong convexity and local Lipschitz continuity of the gradient of convex functions
Local strong convexity and local Lipschitz continuity of the gradient of convex functions R. Goebel and R.T. Rockafellar May 23, 2007 Abstract. Given a pair of convex conjugate functions f and f, we investigate
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
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 informationWeak Convergence Methods for Energy Minimization
Weak Convergence Methods for Energy Minimization Bo Li Department of Mathematics University of California, San Diego E-mail: bli@math.ucsd.edu June 3, 2007 Introduction This compact set of notes present
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More informationOn second order sufficient optimality conditions for quasilinear elliptic boundary control problems
On second order sufficient optimality conditions for quasilinear elliptic boundary control problems Vili Dhamo Technische Universität Berlin Joint work with Eduardo Casas Workshop on PDE Constrained Optimization
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationNumerical Modeling of Methane Hydrate Evolution
Numerical Modeling of Methane Hydrate Evolution Nathan L. Gibson Joint work with F. P. Medina, M. Peszynska, R. E. Showalter Department of Mathematics SIAM Annual Meeting 2013 Friday, July 12 This work
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 informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More informationJuan Vicente Gutiérrez Santacreu Rafael Rodríguez Galván. Departamento de Matemática Aplicada I Universidad de Sevilla
Doc-Course: Partial Differential Equations: Analysis, Numerics and Control Research Unit 3: Numerical Methods for PDEs Part I: Finite Element Method: Elliptic and Parabolic Equations Juan Vicente Gutiérrez
More informationTheory of PDE Homework 2
Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an
More informationOblique derivative problems for elliptic and parabolic equations, Lecture II
of the for elliptic and parabolic equations, Lecture II Iowa State University July 22, 2011 of the 1 2 of the of the As a preliminary step in our further we now look at a special situation for elliptic.
More informationarxiv: v1 [math.na] 15 Aug 2007
A POSTERIORI ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF THE CAHN-HILLIARD EQUATION AND THE HELE-SHAW FLOW XIAOBING FENG AND HAIJUN WU arxiv:78.116v1 [math.na] 15 Aug 7 Abstract. This paper develops
More informationOn the shape of solutions to the Extended Fisher-Kolmogorov equation
On the shape of solutions to the Extended Fisher-Kolmogorov equation Alberto Saldaña ( joint work with Denis Bonheure and Juraj Földes ) Karlsruhe, December 1 2015 Introduction Consider the Allen-Cahn
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differential Equations Eric de Sturler University of Illinois at Urbana-Champaign Read section 8. to see where equations of type (au x ) x = f show up and their (exact) solution
More informationSome PDE s in Image Processing: The Total Variation Flow
Some PDE s in Image Processing: The Total Variation Flow José M. Mazón Abstract We summarize in this lectures some of our results about the Minimizing Total Variation Flow, which have been mainly motivated
More informationWeak Formulation of Elliptic BVP s
Weak Formulation of Elliptic BVP s There are a large number of problems of physical interest that can be formulated in the abstract setting in which the Lax-Milgram lemma is applied to an equation expressed
More informationPriority Programme Density of Convex Intersections and Applications
Priority Programme 1962 Density of Convex Intersections and Applications Michael Hintermüller, Carlos N. Rautenberg, Simon Rösel Non-smooth and Complementarity-based Distributed Parameter Systems: Simulation
More informationarxiv: v1 [math.ap] 13 Mar 2017
1 Mathematical Modeling of Biofilm Development arxiv:173.442v1 [math.ap] 13 Mar 217 Maria Gokieli, Nobuyuki Kenmochi and Marek Niezgódka Interdisciplinary Centre for Mathematical and Computational Modelling,
More informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationFROM VARIATIONAL TO HEMIVARIATIONAL INEQUALITIES
An. Şt. Univ. Ovidius Constanţa Vol. 12(2), 2004, 41 50 FROM VARIATIONAL TO HEMIVARIATIONAL INEQUALITIES Panait Anghel and Florenta Scurla To Professor Dan Pascali, at his 70 s anniversary Abstract A general
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 informationA posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model
A posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model Jan Giesselmann joint work with Ch. Makridakis (Univ. of Sussex), T. Pryer (Univ. of Reading) 9th DFG-CNRS WORKSHOP
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 information************************************* Applied Analysis I - (Advanced PDE I) (Math 940, Fall 2014) Baisheng Yan
************************************* Applied Analysis I - (Advanced PDE I) (Math 94, Fall 214) by Baisheng Yan Department of Mathematics Michigan State University yan@math.msu.edu Contents Chapter 1.
More informationON SOME SINGULAR LIMITS OF HOMOGENEOUS SEMIGROUPS. In memory of our friend Philippe Bénilan
ON SOME SINGULAR LIMITS OF HOMOGENEOUS SEMIGROUPS P. Bénilan, L. C. Evans 1 and R. F. Gariepy In memory of our friend Philippe Bénilan Mais là où les uns voyaient l abstraction, d autres voyaient la vérité
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationRegularity Theory a Fourth Order PDE with Delta Right Hand Side
Regularity Theory a Fourth Order PDE with Delta Right Hand Side Graham Hobbs Applied PDEs Seminar, 29th October 2013 Contents Problem and Weak Formulation Example - The Biharmonic Problem Regularity Theory
More informationVariational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 2012, Northern Arizona University
Variational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 22, Northern Arizona University Some methods using monotonicity for solving quasilinear parabolic
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationOn the p-laplacian and p-fluids
LMU Munich, Germany Lars Diening On the p-laplacian and p-fluids Lars Diening On the p-laplacian and p-fluids 1/50 p-laplacian Part I p-laplace and basic properties Lars Diening On the p-laplacian and
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationParameter robust methods for second-order complex-valued reaction-diffusion equations
Postgraduate Modelling Research Group, 10 November 2017 Parameter robust methods for second-order complex-valued reaction-diffusion equations Faiza Alssaedi Supervisor: Niall Madden School of Mathematics,
More informationThree remarks on anisotropic finite elements
Three remarks on anisotropic finite elements Thomas Apel Universität der Bundeswehr München Workshop Numerical Analysis for Singularly Perturbed Problems dedicated to the 60th birthday of Martin Stynes
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationASYMPTOTIC BEHAVIOUR OF NONLINEAR ELLIPTIC SYSTEMS ON VARYING DOMAINS
ASYMPTOTIC BEHAVIOUR OF NONLINEAR ELLIPTIC SYSTEMS ON VARYING DOMAINS Juan CASADO DIAZ ( 1 ) Adriana GARRONI ( 2 ) Abstract We consider a monotone operator of the form Au = div(a(x, Du)), with R N and
More informationBlow-up directions for quasilinear parabolic equations UNIVERSITY OF TOKYO
UTMS 2006 20 August 18, 2006 Blow-up directions for quasilinear parabolic equations by Yukihiro Seki, Ryuichi Suzuki and Noriaki Umeda T UNIVERSITY OF TOKYO GRADUATE SCHOOL OF MATHEMATICAL SCIENCES KOMABA,
More informationDiscrete entropy methods for nonlinear diffusive evolution equations
Discrete entropy methods for nonlinear diffusive evolution equations 1 Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Joint work with E. Emmrich (TU Berlin), M. Bukal (Zagreb),
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationError estimates for the finite-element approximation of a semilinear elliptic control problem
Error estimates for the finite-element approximation of a semilinear elliptic control problem by Eduardo Casas 1 and Fredi röltzsch 2 1 Departamento de Matemática Aplicada y Ciencias de la Computación
More informationGradient estimates and global existence of smooth solutions to a system of reaction-diffusion equations with cross-diffusion
Gradient estimates and global existence of smooth solutions to a system of reaction-diffusion equations with cross-diffusion Tuoc V. Phan University of Tennessee - Knoxville, TN Workshop in nonlinear PDES
More informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationA POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS
A POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS D. LEYKEKHMAN AND L. B. WAHLBIN Abstract. This paper considers a posteriori error estimates by averaged gradients
More informationOptimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms
Optimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms Marcus Sarkis Worcester Polytechnic Inst., Mass. and IMPA, Rio de Janeiro and Daniel
More informationA Necessary and Sufficient Condition for the Continuity of Local Minima of Parabolic Variational Integrals with Linear Growth
A Necessary and Sufficient Condition for the Continuity of Local Minima of Parabolic Variational Integrals with Linear Growth E. DiBenedetto 1 U. Gianazza 2 C. Klaus 1 1 Vanderbilt University, USA 2 Università
More informationWell-posedness and asymptotic analysis for a Penrose-Fife type phase-field system
0-0 Well-posedness and asymptotic analysis for a Penrose-Fife type phase-field system Buona positura e analisi asintotica per un sistema di phase field di tipo Penrose-Fife Salò, 3-5 Luglio 2003 Riccarda
More informationA posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem
A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem Larisa Beilina Michael V. Klibanov December 18, 29 Abstract
More informationP(E t, Ω)dt, (2) 4t has an advantage with respect. to the compactly supported mollifiers, i.e., the function W (t)f satisfies a semigroup law:
Introduction Functions of bounded variation, usually denoted by BV, have had and have an important role in several problems of calculus of variations. The main features that make BV functions suitable
More informationPERTURBATION THEORY FOR NONLINEAR DIRICHLET PROBLEMS
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 28, 2003, 207 222 PERTURBATION THEORY FOR NONLINEAR DIRICHLET PROBLEMS Fumi-Yuki Maeda and Takayori Ono Hiroshima Institute of Technology, Miyake,
More informationElliptic Kirchhoff equations
Elliptic Kirchhoff equations David ARCOYA Universidad de Granada Sevilla, 8-IX-2015 Workshop on Recent Advances in PDEs: Analysis, Numerics and Control In honor of Enrique Fernández-Cara for his 60th birthday
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationFinite Element Method for Ordinary Differential Equations
52 Chapter 4 Finite Element Method for Ordinary Differential Equations In this chapter we consider some simple examples of the finite element method for the approximate solution of ordinary differential
More informationQuasistatic Nonlinear Viscoelasticity and Gradient Flows
Quasistatic Nonlinear Viscoelasticity and Gradient Flows Yasemin Şengül University of Coimbra PIRE - OxMOS Workshop on Pattern Formation and Multiscale Phenomena in Materials University of Oxford 26-28
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
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 informationChapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached
More informationThe heat equation in time dependent domains with Neumann boundary conditions
The heat equation in time dependent domains with Neumann boundary conditions Chris Burdzy Zhen-Qing Chen John Sylvester Abstract We study the heat equation in domains in R n with insulated fast moving
More informationA posteriori estimators for obstacle problems by the hypercircle method
A posteriori estimators for obstacle problems by the hypercircle method Dietrich Braess 1 Ronald H.W. Hoppe 2,3 Joachim Schöberl 4 January 9, 2008 Abstract A posteriori error estimates for the obstacle
More informationEvolution Under Constraints: Fate of Methane in Subsurface
Evolution Under Constraints: Fate of Methane in Subsurface M. Peszyńska 1 Department of Mathematics, Oregon State University SIAM PDE. Nov.2011 1 Research supported by DOE 98089 Modeling, Analysis, and
More information