Numerical Solution of Nonsmooth Problems and Application to Damage Evolution and Optimal Insulation
|
|
- Jeremy Benson
- 5 years ago
- Views:
Transcription
1 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), M. Milicevic (U Freiburg), R. Nochetto (U Maryland), A. Salgado (U Tennessee), M. Thomas (WIAS Berlin) FoCM 2017, Barcelona, 17 July / 32
2 Nonsmooth Models Obstacle problem Minimize 1 2 Ω subject to u χ Placticity (friction) models t p I K (σ), div C ( ε(u) p ) = f Sparse reconstruction u 2 dx fu dx Ω Minimize x l 1 subject to Ax = b Image denoising Minimize Du + α u g 2 2 Ω 2 / 32
3 Outline 1 Introduction 2 Thin Insulating Films 3 Total Variation Minimization 4 BV Regularized Damage Evolution 5 Computing Insulating Films 6 Summary 3 / 32
4 Thin Insulating Films 4 / 32
5 Thin Insulation Ω Model of insulation: ε \ Ω Ω R d conducting body Ω εl(s) 0 insulating layer of width εl(s) conductivity of layer is ε u = f ε u = 0 u = 0 in Ω c ε PDE system: u C(Ω ε ) with u = f ε u = 0 u = 0 in Ω in Ω ε \ Ω on Ω ε Limit ε 0 leads to Robin BC: u = f u + l n u = 0 in Ω on Ω See: Brézis, Caffarelli & Friedman 80, Acerbi & Buttazzo 86 Operator form of limit model A l u = f with A l : H 1 (Ω) H 1 (Ω) A l u, v = u v dx Ω + l 1 uv ds Ω A l coercive & symmetric 5 / 32
6 Loss of Heat Associated heat equation t u + A l u = 0, u(0) = u 0 Prinicipal eigenvalue determines long-time behaviour λ m,l = inf A lu, u λ m,l u 2 A l u, u u =1 Decay of variance of enthalpy 1 d 2 dt u 2 = A l u, u λ m,l u 2 = u(t) u 0 e λ m,lt λ m,l smallest eigenfrequency Dependence on available mass m = Ω l ds? Optimal distribution l? Class of admissible mass distributions { I m = l L 1 ( Ω) : l 0, Ω } l ds = m 6 / 32
7 Optimal Distribution For fixed m optimize λ m,l over l I m and note interchange λ m = inf l I m λ m,l = inf u =1 Given u H 1 (Ω) optimal thickness l I m is Problem reduces to λ m = l(s) = m inf J m(u), J m (u) = u =1 u(s) u L 1 ( Ω) Ω inf A l u, u l I m u 2 dx + 1 m ( Ω ) 2 u ds Existence of eigenfunctions u m Nonlocal, nondifferentiable, constrained problem Thickness l proportional to trace of u 7 / 32
8 Symmetry Break Recall λ m = min u =1 J m(u), J m (u) = u m u 2 L 1 ( Ω) m λ m continuous, strictly decreasing, λ m 0 as m Dirichlet and Neumann eigenvalues for λ D = min u =1, u Ω =0 u 2 > λ N = min u =1, Ω u dx=0 u 2 λ m λ D for m 0 and λ m0 = λ N for some m 0 > 0 λ D λ N m λ m Theorem (Bucur, Buttazzo & Nitsch 16). Nonradial eigenfunctions u m for Ω = B R (0) iff m < m 0, particularly, optimal mass distributions l leave gaps. Gap symmetric and/or connected? m 0 8 / 32
9 Proof (sketched) radial u m solves 1D Sturm Liouville problem and has no roots u m satisfies PDE ( u m, v) + 1 m u m L ( Ω)( 1 σ(um ), v ) Ω = λ m(u m, v) with σ(u m ) u m, i.e., σ(u m ) = 1 since u m > 0 testing u N with v(x) = c 1 x 2 c 2 gives (u N, 1) Ω = ( u N, v) = λ N (u N, v) = 0 u m L 2,H u 1 N since λ m (u m, u N ) = ( u m, u N ) = λ N (u N, u m ) for ε small u m + εu N L1 ( Ω) = u m L1 ( Ω) implies ( um + εu ) N λ m J m = J m(u m ) + ε 2 u N 2 u m + εu N 1 + ε 2 = λ m + ε 2 λ N 1 + ε 2 contradicts λ m > λ N hence u m not radial 9 / 32
10 Total Variation Minimization 10 / 32
11 BV Model Problem For noisy image g : Ω R find u BV (Ω) minimal for I(u) = Du + α u g 2 2 Proposed by Rudin, Osher & Fatemi 92 Convex, nondifferentiable minimization problem Captures discontiniuties, staircasing effect Ω = Contributions: Elliott & Mikelić 91, Chambolle & P.-L. Lions 97, Feng & Prohl 03, Hintermüller & Kunisch 04, Tsai & Osher 05, Osher, Burger, Goldfarb, Xu & Yin 05, Wu & Tai 09, Chambolle & Pock 11, Wang & Lucier 11, Stamm & Wihler 15, / 32
12 Properties of BV Function v BV (Ω) if v L 1 (Ω) and Dv = Dv (Ω) = sup Ω ξ C 0 (Ω;Rd ), ξ 1 Examples: (i) v(x) = sign(x), Dv = 2δ 0, Dv (Ω) = 2 (ii) v(x) = χ A (x), Dv (Ω) = Per Ω (A) v div ξ dx < Ω Existence of solutions by convexity and compactness W 1,1 (Ω) BV (Ω) strictly and v L 1 (Ω) = Dv (Ω) C (Ω) BV (Ω) weakly but not strongly dense Intermediate convergence in BV (Ω): v j v in L 1 (Ω) & Dv j (Ω) Dv (Ω) 12 / 32
13 Aspects and Tools Choice of finite element space V h BV (Ω) Convergence rates unter appropriate regularity conditions A posteriori error estimates and local mesh refinement Effective iterative numerical solution Strong convexity: For arbitrary v BV (Ω) L 2 (Ω) α 2 u v 2 I(v) I(u) σ I (v, w) Strong duality: Sharp relation I(v) D(q) with q H N (div; Ω) and D(q) = 1 2α div q + αg 2 + α 2 g 2 I K1 (0)(q) Combination: For minimizer u and arbitrary v and q α 2 u v 2 I(v) I(u) I(v) D(q) v w 13 / 32
14 Convergence Rates P0 FE fails as perimeter approximated incorrectly: h u h L 1 χ {x<0} = Du h (Ω) 2 P1 FE quasi-interpolant via regularization and interpolation Lemma (B., Nochetto & Salgado 12). There exists ũ h,ε S 1 (T h ) with ũ h,ε L 1 (1 + cε + chε 1 ) Du (Ω), u ũ h,ε L 1 c(h 2 ε 1 + ε) Du (Ω). Implies rate O(h 1/4 ) with ε = h 1/2 : α 2 u u h 2 I(u h ) I(u) inf I(v h ) I(u) ch 1/2 v h S 1 (T h ) 1.5 Optimal rate is O(h 1/2 ) for generic u BV (Ω) TVD property ũ h,ε Du (Ω) yields optimality u P h u h h = 1/10 h = 1/20 h = 1/ / 32
15 Adaptivity Computable localized error control via approximation of dual (B. 15): α 2 u u h 2 I(u h ) D(p h ) = u h L 1 = Ω T T h η T (u h, p h ), p h u h dx + 1 2α div p h α(u h g) 2 provided p h H N (div; Ω) with p h 1; note η T (u h, p h ) 0 Nearly optimal rate h 1/2 N 1/4 Dual solution in P1 vector fields Better: H(div; Ω) FE spaces 15 / 32
16 Solution via ADMM Decouple nondifferentiability from gradient via p = u L τ (u, p; λ) = p dx + α 2 u g 2 + (λ; u p) H + τ 2 u p 2 H Ω Formally, e.g., on discrete level, for arbitrary τ 0 inf sup u,p λ L τ (u, p; λ) = inf u I(u) Choice of Hilbert space H: weighted L 2 norm on FE spaces Iterative solution via decoupled minimization: (Glowinski & Morocco 76, Gabay & Mercier 76, He & Yuan 12, Deng & Yin 16) Algorithm. Choose (u 0, λ 0 ) and τ > 0; set j = 1. (1) Compute minimizer p j for p L τ ( u j 1, p; λ j 1). (2) Compute minimizer u j for u L τ ( u, p j ; λ j 1 ). (3) Update λ j 1 via λ j = λ j 1 + τ( u j p j ). 16 / 32
17 Properties Unconditional convergence if τ > 0 Proof holds for sequence of decreasing step sizes τ j+1 τ j Monotonicity of residuals Error control via residual R j = λ j λ j 1 2 H + τ 2 j (u j u j 1 ) 2 H u u j 2 + u j u j 1 2 C 0 τ 1 j R 1/2 j Linear convergence of residuals R j+1 γr j, 0 < γ < 1 Motivates adaptive step size adjustment (B. & Milicevic 17): Start with τ 1 and γ (0, 1) Accept τ if γ-reduction satisfied, decrease τ otherwise If τ = 1 increase γ and re-start 17 / 32
18 Experiment Comparison to ADMM and accelerated version (Goldfarb et al. 13) Hyphens indicate more than 10 4 steps Variable step size variant best for large initial step sizes Restricted to convex problems 18 / 32
19 Subdifferential Flow Gradient flows most general and robust tool to minimize Formal PDE for TV flow and backward Euler method t u = div u u, d tu k = div uk u k Unconditional energy stability by convexity d t u k 2 + d t Ω u k dx d t u k 2 u k + Ω u k d tu k dx = 0 Practical variant semiimplicit and regularized d t u k u k = div u k 1 ε Violates monotonicity property and might depend critically on ε 19 / 32
20 Unconditional Stability Continuous differentiation formula d dt a(t) = a(t)a (t) a(t) Discrete product and quotient rule = 1 2 d dt a(t) 2 a(t) d t a k b k = (d t a k )b k + a k 1 (d t b k ), d t 1 c k = d tc k c k 1 c k Binomial formula Crucial formula d t a k = d t a k 2 a k 2a k d t a k = d t a k 2 + τ d t a k 2 = d t a k 2 a k 1 + ak 2 d t 1 a k = d t a k 2 a k 1 ak d t a k a k 1 = 1 d t a k 2 2 a k 1 τ d ta k 2 2 a k 1 1 d t a k 2 2 a k 1 Implies unconditional convergence of semi-implicit scheme 20 / 32
21 Experiments Change of height proportional to mean curvature of level set Evolution of characteristic function of disk u 0 = χ B1 Exact solution given by u(t, x) = (1 2t) + χ B1 (x) Error u u h L (L 2 ) depends on ε External Movie implicit ε = h 1/2 ε = h ε = h / 32
22 BV Regularized Damage Evolution 22 / 32
23 BV Regularized Damage Model Energy and dissipation for displacement ϕ and damage variable z: E(ϕ, z) = f (z)ε(ϕ) : Cε(φ) dx + Dz (Ω) + I [0,1] (z) dx Ω Ω { ϱ v if v (, 0], R(v) = R(v) dx, R(v) = + if v > 0. Ω Evolution of q = (ϕ, z) via Biot s equation 0 q E(q) + q R( q) Energetic solutions are globally stable and fulfill energy balance E ( q(t) ) E ( q ) + R ( z z(t) ) E ( q(t) ) ( ) ( ) t + Diss R z, [s, t] = E q(s) + r E ( q(r) ) dr with Diss R (z, [s, t]) = sup { N j=1 R( z(t j ) z(t j 1 ) ), t 0 <... < t N }. derivative-free formulation existence theories: Mielke 05; partial damage: Thomas 13 s 23 / 32
24 Experiment Time-stepping via decoupled minimization External Link 24 / 32
25 Computing Insulating Films 25 / 32
26 Gradient Flow Minimization problem J m (u) = u m u 2 L 1 ( Ω) subject to u = 1 Regularized gradient flow with u ε = (u 2 + ε 2 ) 1/2 ( t u, v) + ( u, v) + 1 m u L 1 ε ( u ), v = µ(u, v) u ε Ω u(t) = 1 implies t u L 2 u; RHS disappears if v L 2 u Algorithm (B. 17). Given u k 1 compute u k with d t u k L 2 u k 1 and (d t u k, v) + ( u k, v) + 1 ( u k ) m uk 1 L 1 ε u k 1, v ε for all v with v L 2 u k 1 with quotient d t u k = (u k u k 1 )/τ. Iteration (essentially) unconditionally stable Ω = 0 26 / 32
27 Spatial Discretization For triangulation T h of Ω use P1 FE-space S 1 (T h ) = { } v h C(Ω) : v h T P 1 (T ) for all T T h. and discretized functional J m,ε,h (u h ) = u h ( ) 2 β z u h (z) ε m z N h Ω Stability estimate carries over to J m,ε,h Lemma (B. 17). If eigenfunction u = u m satisfies u H 2 then 0 min u h S 1 (T h ) J m (u h ) J m (u) ch u 2 H 2 (Ω), Jm,ε,h (u h ) J m (u h ) c( uh H1 (Ω) + 1) 2 (ε + h 1/2 ). Consistency estimates imply Γ convergence via density 27 / 32
28 Experiments Dirichlet and Neumann eigenfunctions on unit disk Optimal and constant l with m = 0.4 on unit disk Optimal and constant l with m = 0.1 on unit square (λ D = 2π 2, λ N = π 2 ) 28 / 32
29 Iteration Numbers Ω = B 1 (0) R 2, m = 1/2, h = 2 l, ε = h/10, τ = 1 l #N h #T h K stop λ m,h Ω = (0, 1) 2 R 2, m = 1/8, h = 2 l, ε = h/10, τ = 1 l #N h #T h K stop λ m,h Fixed stopping criterion: d t u k h ε stop = / 32
30 Rotational Shapes Symmetrical ellipsoids E aa with radii (a, a, r a ), a [1, 1.5] Optimal assembled ellipsoid E ab with radii (a, b, r ab ) Fixed volume c 0 = B 1 and mass m = 6.0 λ m(e aa) λ m(e opt ab ) (i) (ii) (i) a (ii) High resolution via two-dimensional reduction Egg shape optimal among convex shapes? 30 / 32
31 Summary 31 / 32
32 Summary Reduced convergence rates for nondifferentiable problems Iterative solution via regularized gradient flows Semi-implicit discretization unconditionally stable Future topics Convergence of adaptive schemes Optimal isolating bodies Error estimates for damage evolution Further information 32 / 32
SÖREN BARTELS AND GIUSEPPE BUTTAZZO
NUMERICAL SOLUTION OF A NONLINEAR EIGENVALUE PROBLEM ARISING IN OPTIMAL INSULATION arxiv:1708.03762v1 [math.na] 12 Aug 2017 Abstract. The optimal insulation of a heat conducting body by a thin film of
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 informationSymmetry breaking for a problem in optimal insulation
Symmetry breaking for a problem in optimal insulation Giuseppe Buttazzo Dipartimento di Matematica Università di Pisa buttazzo@dm.unipi.it http://cvgmt.sns.it Geometric and Analytic Inequalities Banff,
More informationALTERNATING DIRECTION METHOD OF MULTIPLIERS WITH VARIABLE STEP SIZES
ALTERNATING DIRECTION METHOD OF MULTIPLIERS WITH VARIABLE STEP SIZES Abstract. The alternating direction method of multipliers (ADMM) is a flexible method to solve a large class of convex minimization
More informationA posteriori error control for the binary Mumford Shah model
A posteriori error control for the binary Mumford Shah model Benjamin Berkels 1, Alexander Effland 2, Martin Rumpf 2 1 AICES Graduate School, RWTH Aachen University, Germany 2 Institute for Numerical Simulation,
More informationNUMERICAL SOLUTION OF NONSMOOTH PROBLEMS
NUMERICAL SOLUTION OF NONSMOOTH PROBLEMS Abstract. Nonsmooth minimization problems arise in problems involving constraints or models describing certain discontinuities. The existence of solutions is established
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 informationANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE INTO CARTOON PLUS TEXTURE. C.M. Elliott and S.A.
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 6, Number 4, December 27 pp. 917 936 ANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE
More informationSubdiffusion in a nonconvex polygon
Subdiffusion in a nonconvex polygon Kim Ngan Le and William McLean The University of New South Wales Bishnu Lamichhane University of Newcastle Monash Workshop on Numerical PDEs, February 2016 Outline Time-fractional
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 informationAdaptive approximation of eigenproblems: multiple eigenvalues and clusters
Adaptive approximation of eigenproblems: multiple eigenvalues and clusters Francesca Gardini Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/gardini Banff, July 1-6,
More informationVariational Image Restoration
Variational Image Restoration Yuling Jiao yljiaostatistics@znufe.edu.cn School of and Statistics and Mathematics ZNUFE Dec 30, 2014 Outline 1 1 Classical Variational Restoration Models and Algorithms 1.1
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationIterative regularization of nonlinear ill-posed problems in Banach space
Iterative regularization of nonlinear ill-posed problems in Banach space Barbara Kaltenbacher, University of Klagenfurt joint work with Bernd Hofmann, Technical University of Chemnitz, Frank Schöpfer and
More informationControl of Interface Evolution in Multi-Phase Fluid Flows
Control of Interface Evolution in Multi-Phase Fluid Flows Markus Klein Department of Mathematics University of Tübingen Workshop on Numerical Methods for Optimal Control and Inverse Problems Garching,
More informationAdaptive discretization and first-order methods for nonsmooth inverse problems for PDEs
Adaptive discretization and first-order methods for nonsmooth inverse problems for PDEs Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Barbara Kaltenbacher, Tuomo Valkonen,
More informationDual methods for the minimization of the total variation
1 / 30 Dual methods for the minimization of the total variation Rémy Abergel supervisor Lionel Moisan MAP5 - CNRS UMR 8145 Different Learning Seminar, LTCI Thursday 21st April 2016 2 / 30 Plan 1 Introduction
More informationParameter Identification in Partial Differential Equations
Parameter Identification in Partial Differential Equations Differentiation of data Not strictly a parameter identification problem, but good motivation. Appears often as a subproblem. Given noisy observation
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 informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationCoordinate Update Algorithm Short Course Operator Splitting
Coordinate Update Algorithm Short Course Operator Splitting Instructor: Wotao Yin (UCLA Math) Summer 2016 1 / 25 Operator splitting pipeline 1. Formulate a problem as 0 A(x) + B(x) with monotone operators
More informationContinuous Primal-Dual Methods in Image Processing
Continuous Primal-Dual Methods in Image Processing Michael Goldman CMAP, Polytechnique August 2012 Introduction Maximal monotone operators Application to the initial problem Numerical illustration Introduction
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 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 informationA generalized MBO diffusion generated motion for constrained harmonic maps
A generalized MBO diffusion generated motion for constrained harmonic maps Dong Wang Department of Mathematics, University of Utah Joint work with Braxton Osting (U. Utah) Workshop on Modeling and Simulation
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 informationFinite difference method for heat equation
Finite difference method for heat equation 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 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 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 informationContraction Methods for Convex Optimization and Monotone Variational Inequalities No.16
XVI - 1 Contraction Methods for Convex Optimization and Monotone Variational Inequalities No.16 A slightly changed ADMM for convex optimization with three separable operators Bingsheng He Department of
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 informationAccelerated primal-dual methods for linearly constrained convex problems
Accelerated primal-dual methods for linearly constrained convex problems Yangyang Xu SIAM Conference on Optimization May 24, 2017 1 / 23 Accelerated proximal gradient For convex composite problem: minimize
More informationApproximation of self-avoiding inextensible curves
Approximation of self-avoiding inextensible curves Sören Bartels Department of Applied Mathematics University of Freiburg, Germany Joint with Philipp Reiter (U Duisburg-Essen) & Johannes Riege (U Freiburg)
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 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 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 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 informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 155 A posteriori error estimates for stationary slow flows of power-law fluids Michael Bildhauer,
More informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) LMS
More informationOPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL
OPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL EDUARDO CASAS, ROLAND HERZOG, AND GERD WACHSMUTH Abstract. Semilinear elliptic optimal control
More informationOn the Optimal Insulation of Conductors 1
JOURNAL OF OPTIMIZATION THEORY AND APPLICATIONS: Vol. 100, No. 2, pp. 253-263, FEBRUARY 1999 On the Optimal Insulation of Conductors 1 S. J. COX, 2 B. KAWOHL, 3 AND P. X. UHLIG 4 Communicated by K. A.
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationOn a multiscale representation of images as hierarchy of edges. Eitan Tadmor. University of Maryland
On a multiscale representation of images as hierarchy of edges Eitan Tadmor Center for Scientific Computation and Mathematical Modeling (CSCAMM) Department of Mathematics and Institute for Physical Science
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 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 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 informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
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 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 informationI P IANO : I NERTIAL P ROXIMAL A LGORITHM FOR N ON -C ONVEX O PTIMIZATION
I P IANO : I NERTIAL P ROXIMAL A LGORITHM FOR N ON -C ONVEX O PTIMIZATION Peter Ochs University of Freiburg Germany 17.01.2017 joint work with: Thomas Brox and Thomas Pock c 2017 Peter Ochs ipiano c 1
More informationAMS 212A Applied Mathematical Methods I Appendices of Lecture 06 Copyright by Hongyun Wang, UCSC. ( ) cos2
AMS 22A Applied Mathematical Methods I Appendices of Lecture 06 Copyright by Hongyun Wang UCSC Appendix A: Proof of Lemma Lemma : Let (x ) be the solution of x ( r( x)+ q( x) )sin 2 + ( a) 0 < cos2 where
More informationAdaptive Primal Dual Optimization for Image Processing and Learning
Adaptive Primal Dual Optimization for Image Processing and Learning Tom Goldstein Rice University tag7@rice.edu Ernie Esser University of British Columbia eesser@eos.ubc.ca Richard Baraniuk Rice University
More informationRecent Developments of Alternating Direction Method of Multipliers with Multi-Block Variables
Recent Developments of Alternating Direction Method of Multipliers with Multi-Block Variables Department of Systems Engineering and Engineering Management The Chinese University of Hong Kong 2014 Workshop
More informationSingular Diffusion Equations With Nonuniform Driving Force. Y. Giga University of Tokyo (Joint work with M.-H. Giga) July 2009
Singular Diffusion Equations With Nonuniform Driving Force Y. Giga University of Tokyo (Joint work with M.-H. Giga) July 2009 1 Contents 0. Introduction 1. Typical Problems 2. Variational Characterization
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 informationPartial Differential Equations, 2nd Edition, L.C.Evans The Calculus of Variations
Partial Differential Equations, 2nd Edition, L.C.Evans Chapter 8 The Calculus of Variations Yung-Hsiang Huang 2018.03.25 Notation: denotes a bounded smooth, open subset of R n. All given functions are
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
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 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 informationConvergence rates of convex variational regularization
INSTITUTE OF PHYSICS PUBLISHING Inverse Problems 20 (2004) 1411 1421 INVERSE PROBLEMS PII: S0266-5611(04)77358-X Convergence rates of convex variational regularization Martin Burger 1 and Stanley Osher
More informationA posteriori error estimates in FEEC for the de Rham complex
A posteriori error estimates in FEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationSplitting Techniques in the Face of Huge Problem Sizes: Block-Coordinate and Block-Iterative Approaches
Splitting Techniques in the Face of Huge Problem Sizes: Block-Coordinate and Block-Iterative Approaches Patrick L. Combettes joint work with J.-C. Pesquet) Laboratoire Jacques-Louis Lions Faculté de Mathématiques
More informationZdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)
Navier-Stokes equations with constrained L 2 energy of the solution Zdzislaw Brzeźniak Department of Mathematics University of York joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York) Stochastic
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 informationCLASSICAL ITERATIVE METHODS
CLASSICAL ITERATIVE METHODS LONG CHEN In this notes we discuss classic iterative methods on solving the linear operator equation (1) Au = f, posed on a finite dimensional Hilbert space V = R N equipped
More informationConvergence of the MAC scheme for incompressible flows
Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba Calanque de Cortiou
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 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 informationObstacle Problems Involving The Fractional Laplacian
Obstacle Problems Involving The Fractional Laplacian Donatella Danielli and Sandro Salsa January 27, 2017 1 Introduction Obstacle problems involving a fractional power of the Laplace operator appear in
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR TERM
Internat. J. Math. & Math. Sci. Vol. 22, No. 3 (999 587 595 S 6-72 9922587-2 Electronic Publishing House ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationLecture No 2 Degenerate Diffusion Free boundary problems
Lecture No 2 Degenerate Diffusion Free boundary problems Columbia University IAS summer program June, 2009 Outline We will discuss non-linear parabolic equations of slow diffusion. Our model is the porous
More informationBeyond Heuristics: Applying Alternating Direction Method of Multipliers in Nonconvex Territory
Beyond Heuristics: Applying Alternating Direction Method of Multipliers in Nonconvex Territory Xin Liu(4Ð) State Key Laboratory of Scientific and Engineering Computing Institute of Computational Mathematics
More informationConvex Hodge Decomposition and Regularization of Image Flows
Convex Hodge Decomposition and Regularization of Image Flows Jing Yuan, Christoph Schnörr, Gabriele Steidl April 14, 2008 Abstract The total variation (TV) measure is a key concept in the field of variational
More informationBias-free Sparse Regression with Guaranteed Consistency
Bias-free Sparse Regression with Guaranteed Consistency Wotao Yin (UCLA Math) joint with: Stanley Osher, Ming Yan (UCLA) Feng Ruan, Jiechao Xiong, Yuan Yao (Peking U) UC Riverside, STATS Department March
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 informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationShape optimization problems for variational functionals under geometric constraints
Shape optimization problems for variational functionals under geometric constraints Ilaria Fragalà 2 nd Italian-Japanese Workshop Cortona, June 20-24, 2011 The variational functionals The first Dirichlet
More informationPreconditioned Alternating Direction Method of Multipliers for the Minimization of Quadratic plus Non-Smooth Convex Functionals
SpezialForschungsBereich F 32 Karl Franzens Universita t Graz Technische Universita t Graz Medizinische Universita t Graz Preconditioned Alternating Direction Method of Multipliers for the Minimization
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 informationStrauss conjecture for nontrapping obstacles
Chengbo Wang Joint work with: Hart Smith, Christopher Sogge Department of Mathematics Johns Hopkins University Baltimore, Maryland 21218 wangcbo@jhu.edu November 3, 2010 1 Problem and Background Problem
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. I - Basic concepts. Alberto Bressan. Department of Mathematics, Penn State University
Hyperbolic Systems of Conservation Laws in One Space Dimension I - Basic concepts Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 The Scalar Conservation
More informationLINEARIZED BREGMAN ITERATIONS FOR FRAME-BASED IMAGE DEBLURRING
LINEARIZED BREGMAN ITERATIONS FOR FRAME-BASED IMAGE DEBLURRING JIAN-FENG CAI, STANLEY OSHER, AND ZUOWEI SHEN Abstract. Real images usually have sparse approximations under some tight frame systems derived
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 informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationA fast nonstationary preconditioning strategy for ill-posed problems, with application to image deblurring
A fast nonstationary preconditioning strategy for ill-posed problems, with application to image deblurring Marco Donatelli Dept. of Science and High Tecnology U. Insubria (Italy) Joint work with M. Hanke
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik A finite element approximation to elliptic control problems in the presence of control and state constraints Klaus Deckelnick and Michael Hinze Nr. 2007-0
More informationRate optimal adaptive FEM with inexact solver for strongly monotone operators
ASC Report No. 19/2016 Rate optimal adaptive FEM with inexact solver for strongly monotone operators G. Gantner, A. Haberl, D. Praetorius, B. Stiftner Institute for Analysis and Scientific Computing Vienna
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 informationConvergence rate estimates for the gradient differential inclusion
Convergence rate estimates for the gradient differential inclusion Osman Güler November 23 Abstract Let f : H R { } be a proper, lower semi continuous, convex function in a Hilbert space H. The gradient
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
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 informationSparse Optimization Lecture: Dual Methods, Part I
Sparse Optimization Lecture: Dual Methods, Part I Instructor: Wotao Yin July 2013 online discussions on piazza.com Those who complete this lecture will know dual (sub)gradient iteration augmented l 1 iteration
More informationTime domain boundary elements for dynamic contact problems
Time domain boundary elements for dynamic contact problems Heiko Gimperlein (joint with F. Meyer 3, C. Özdemir 4, D. Stark, E. P. Stephan 4 ) : Heriot Watt University, Edinburgh, UK 2: Universität Paderborn,
More informationEXISTENCE RESULTS FOR QUASILINEAR HEMIVARIATIONAL INEQUALITIES AT RESONANCE. Leszek Gasiński
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 2007 pp. 409 418 EXISTENCE RESULTS FOR QUASILINEAR HEMIVARIATIONAL INEQUALITIES AT RESONANCE Leszek Gasiński Jagiellonian
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More information