A generalized MBO diffusion generated motion for constrained harmonic maps
|
|
- Evan Tyler
- 5 years ago
- Views:
Transcription
1 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 of Interface-related Problems IMS, NUS, Singapore, May. 3, 2018
2 Motion by mean curvature 1 Mean curvature flow arises in a variety of physical applications: Related to surface tension A model for the formation of grain boundaries in crystal growth Some ideas for numerical computation: We could parameterize the surface and compute H = 1 2 ˆn If the surface is implicitly defined by the equation F(x, y, z) = 0, then mean curvature can be computed H = 1 ( ) F 2 F
3 MBO diffusion generated motion 2 In 1989, Merriman, Bence, and Osher (MBO) developed an iterative method for evolving an interface by mean curvature. Repeat until convergence: Step 1. Solve the Cauchy problem for the diffusion equation (heat equation) u t = u u(x, t = 0) = χ D, with initial condition given by the indicator function χ D of a domain D until time τ to obtain the solution u(x, τ). Step 2. Obtain a domain D new by thresholding: D new = { x R d : u(x, τ) 1 2 }.
4 How to understand the MBO method? 3 Intuitively, from pictures, one can easily see: diffusion quickly blunts sharp points on the boundary and diffusion has little effect on the flatter parts of the boundary. Formally, consider a point P D. In local polar coordinates, the diffusion equation is given by u t = 1 r u r + 2 u r u 2 r 2 θ. 2 Considering local symmetry, we have u t = 1 r u r + 2 u r 2 = H u r + 2 u r 2. The 1 level set will move in the normal 2 direction with velocity given by the mean curvature, H. Initial t = t = t = 0.01 t = 0.02
5 A variational point of view: Modica+Mortola 4 Define the energy J ε (u) = Ω 1 2 u W (u) du ε2 ( where W (u) = 1 4 u 2 u ) 2 is a double well potential. Theorem (Modica+ Mortola, 1977) A minimizing sequence (u ε ) converges (along a subsequence) to χ D in L 1 for some D Ω. Furthermore, εj ε (u ε ) Hd 1 ( D) as ε 0. Gradient flow. The L 2 gradient flow of J ε gives the Allen-Cahn equation: u t = u 1 ε 2 W (u) in Ω.
6 A variational point of view: 5 Operator/energy splitting. Repeat the following two steps until convergence: Step 1. Solve the diffusion equation until time τ with initial condition u(x, t = 0) = χ D t u = u Step 2*. Solve the (pointwise defined!) equation until time τ: φ t = W (φ)/ε 2, φ(x, 0) = u(x, τ), in Ω. Step 2. Rescaling t = ε 2 t, we have as ε 0, ε 2 τ. So, Step 2* is equivalent to thresholding: { 1 if φ(x, 0) > 1/2 φ(x, ) = 0 if φ(x, 0) < 1/2.
7 Analysis, extensions, applications, connections, and computation 6 Proof of convergence of the MBO method to mean curvature flow [Evans1993, Barles and Georgelin 1995, Chambolle and Novaga 2006, Laux and Swartz 2017, Swartz and Yip 2017, Laux and Yip 2018]. Multi-phase problems with arbitrary surface tensions [Esedoglu and Otto 2015, Laux and Otto 2016] Numerical algorithms [Ruuth 1996, Ruuth 1998, Jiang et al. 2017] Area or volume preserving interface motion [Ruuth and Weston 2003] Image processing [Esedoglu et al. 2006, Merkurjev et al. 2013, Wang et al. 2017, Jacobs 2017] Problems of anisotropic interface motion [Merriman et al. 2000, Ruuth et al. 2001, Bonnetier et al. 2010, Elsey et al. 2016] Diffusion generated motion using signed distance function [Esedoglu et al. 2009] High order geometric motion [Esedoglu et al. 2008] Harmonic map of heat flow [E and Wang 2000, Wang et al. 2001, Ruuth et al. 2002] Nonlocal threshold dynamics method [Caffarelli and Souganidis 2010] Wetting problem on solid surfaces [Xu et al. 2017] Graph partitioning and data clustering [Van Gennip et al. 2013] Auction dynamics [Jacobs et al. 2017] Centroidal Voronoi Tessellation [Du et al. 1999] Quad meshing [Viertel and Osting 2017]
8 Generalized energies 7 Let Ω R d be a bounded domain with smooth boundary. Let T R k be "target set" and f : R k R + be a smooth function such that T = f 1 (0). T is the set of global minimizers of f. Roughly, we want f (x) dist 2 (x; T ). Consider a generalized variational problem, inf E(u) where E(u) = 1 u 2 dx u:ω T 2 Ω Relax the energy to obtain: 1 min E ε(u) where E ε(u) = u H 1 (Ω,R k ) Ω 2 u f (u(x)) dx. ε2 Examples. k T f(x) comment 1 1 {±1} 4 (x 2 1) 2 Allen-Cahn 2 S ( x 2 1) 2 Ginzburg-Landau k S k ( x 2 1) 2 n 2 1 O(n) 4 x t x I n 2 F Orthogonal matrix valued fields 1 k coordinate axes, Σ k 4 i j x i 2 xj 2 Dirichlet partitions RP 2 Landau-de Gennes model.
9 Orthogonal matrix valued fields 8 Let O n M n = R n n be the group of orthogonal matrices. inf E(A), where E(A) := 1 A 2 F dx. A:Ω O n 2 Relaxation: min E ε(a), where E ε(a) := A H 1 (Ω,M n) Ω Ω 1 2 A 2 F + 1 4ε 2 At A I n 2 F dx. The penalty term can be written: 1 4ε 2 At A I n 2 F = 1 n W (σ ε 2 i (A)), where W (x) = 1 4 i=1 Gradient Flow. The gradient flow of E ε is Special cases. ta = A E 2,ε (A) = A ε 2 A(A t A I n). For n = 1, we recover Allen-Cahn equation. For n = 2, if the initial condition is taken to be in SO(2) = S 1, we recover the complex Ginzburg-Landau equation. ( x 2 1) 2.
10 Diffusion generated motion for O n valued fields. 9 E ε(a) := Ω 1 2 A 2 F + 1 4ε 2 At A I n 2 F dx. k T f(x) comment n 2 O(n) 1 x I 4 n 2 F Orthogonal matrix valued fields Lemma. The nearest-point projection map, Π T : R n n T, for T = O n is given by Π T A = A(A t A) 1 2 = UV t, where A has the singular value decomposition, A = UΣV t. Diffusion generated method. For i = 1, 2,..., Step 1. Solve the diffusion equation until time τ with initial value given by A i (x): tu = u, u(x, t = 0) = A i (x). Step 2. Point-wise, apply the nearest-point projection map: A i+1 (x) = Π T u(x, τ).
11 Computational Example I: flat torus, n = 2 10 closed line defect vol. const. closed line defect parallel lines defect parallel lines defect O(n) = SO(n) SO (n), SO(2) = S 1 x is yellow det(a(x)) = 1 A(x) SO(n) x is blue det(a(x)) = 1 A(x) SO (n).
12 Computational Example II: sphere, n = 3 11 shrinking on sphere vol. const. on sphere
13 Computational Example III: peanut, n = 3 12 x(t, θ) = (3t t 3 ), y(t, θ) = 1 (1 + x 2 2 )(4 x 2 ) cos(θ), z(t, θ) = 1 (1 + x 2 2 )(4 x 2 ) sin(θ) peanut with closed geodesic
14 Lyapunov functional for MBO iterates 13 Motivated by [Esedoglu + Otto, 2015], we define the functional E τ : H 1 (Ω, M n ) R, given by E τ (A) := 1 n A, e τ A F dx τ Ω Here, e τ A denotes the solution to the heat equation at time τ with initial condition at time t = 0 given by A = A(x). Denoting the spectral norm by A 2 = σ max (A), the convex hull of O n is K n = conv O n = {A M n : A 2 1}.
15 Lyapunov functional for MBO iterates 14 Lemma. The functional E τ has the following elementary properties. (i) For A L 2 (Ω, O n ), (ii) E τ (A) is concave. (iii) We have E τ (A) = E(A) + O(τ). min E τ (A) = min E τ (A). A L 2 (Ω,O n) A L 2 (Ω,K n) (iv) E τ (A) is Fréchet differentiable with derivative L τ A : L (Ω, M n ) R at A in the direction B given by L τ A(B) = 2 e τ A, B F dx. τ Ω
16 Stability 15 The sequential linear programming approach to minimizing E τ (A) subject to A L (Ω, K n ) is to consider a sequence of functions {A s } s=0 which satisfies A s+1 = arg min A L (Ω,K Lτ A s (A), n) A 0 L (Ω, O n ) given. Lemma. If e Asτ = UΣV t, the solution to the linear optimization problem, min A L (Ω,K n) Lτ A s (A). is attained by the function A = UV t L (Ω, O n ). Thus, A s L (Ω, O n ) for all s 0 and these are precisely the iterations generated by the generalized MBO diffusion generated motion! Theorem (Stability). [Osting + W., 2017] The functional E τ is non-increasing on the iterates {A s } s=1, i.e., E τ (A s+1 ) E τ (A s ).
17 Convergence 16 We consider a discrete grid Ω = {x i } Ω i=1 Ω and a standard finite difference approximation of the Laplacian,, on Ω. For A: Ω O n, define the discrete functional Ẽ τ (A) = 1 τ 1 A i, (e τ A) i F x i Ω and its linearization by L τ A(B) = 2 τ B i, (e τ A) i F. x i Ω
18 Convergence 17 Theorem (Convergence for n = 1.) [Osting + W., 2017] Let n = 1. Non-stationary iterations of the generalized MBO diffusion generated motion strictly decrease the value of Ẽ τ and since the state space is finite, {±1} Ω, the algorithm converges in a finite number of iterations. Furthermore, for m := e τ, each iteration reduces the value of J by at least 2m, so the total number of iterations is less than Ẽ τ (A 0 )/2m. Theorem (Convergence for n 2.) [Osting + W., 2017] Let n 2. The non-stationary iterations of the generalized MBO diffusion generated motion strictly decrease the value of Ẽ τ. For a given initial condition A 0 : Ω O n, there exists a partition Ω = Ω + Ω and an S N such that for s S, { +1 x i det A s (x i ) = Ω +. 1 x i Ω Lemma. dist (SO(n), SO (n)) = 2.
19 Volume constrained problem 18 Definitions: D n : A diagonal matrix with diagonal entries 1 everywhere except in the n th position, where it is 1. T + (A) = { UV t, if det A > 0, UD n V t, if det A < 0. T (A) = { UDn V t, if det A > 0, UV t, if det A < 0. E(x) = T + (A(τ, x)) T (A(τ, x)), A(τ, x) F, Ω λ + = {x : E(x) λ}, Ω λ = Ω \ Ω λ +.
20 Volume constrained problem 19 Similar as that in [Ruuth+Weston, 2003], we treat the volume of Ω λ + as a function of λ and identify the value λ such that f (λ) = V. For the matrix case, the following lemma leads us to the optimal choice: Lemma. Assume λ 0 satisfies f (λ 0 ) = V. Then, { B T + (A(τ, x)) if E(x) λ 0 =. T (A(τ, x)) if E(x) < λ 0 attains the minimum in min L τ A(B) B L (Ω,O n) s.t. vol({x : B(x) SO(n)}) = V. Theorem (Stability). [Osting+W., 2017] For any τ > 0, the functional E τ, is non-increasing on the volume-preserving iterates {A s} s=1, i.e., E τ (A s+1 ) E τ (A s).
21 Dirichlet partitions 20 Let U R d be either an open bounded domain with Lipschitz boundary. Dirichlet k-partition: A collection of k disjoint open sets, U 1, U 2,..., U k U, attains inf U l U U l U m= k l=1 λ 1 (U l ) where λ 1 (U) := min E(u). u H0 1 (U) u L 2 (U) =1 E(u) = U u 2 dx is the Dirichlet energy and u L 2 (U) := ( U u2 (x)dx) 1 2. λ 1 (U) is the first Dirichlet eigenvalue of the Laplacian,. Monotonicity of eigenvalues = U = k l=1u l.
22 A mapping formulation of Dirichlet partitions [Cafferelli and Lin 2007] 21 Consider vector valued functions u = (u 1, u 2,..., u k ), that takes values in the singular space, Σ k, given by the coordinate axes, k Σ k := x Rk : xi 2 xj 2 = 0. i j The Dirichlet partition problem for U is equivalent to the mapping problem { } min E(u): u = (u 1,..., u k ) H0 1 (U; Σ k ), ul 2 (x) dx = 1 l [k], U where E(u) := k l=1 U u l 2 dx is the Dirichlet energy of u and H0 1(U; Σ k ) = {u H0 1(U, Rk ) : u Σ k a.e.}. Refer to minimizers u as ground states, and WLOG take u > 0 and quasi-continuous. U = l u is a ground state U l with U l = u 1 l ((0, )) for l [k] is a Dirichlet partition.
23 Diffusion generated method for computing Dirichlet partitions 22 Relaxed problem: 1 E ε(u) := Ω 2 u 2 F + 1 u 2 4ε 2 i (x)u 2 j (x) dx. i j k T f(x) comment 1 k coordinate axes, Σ k 4 i j x 2 i x 2 j Dirichlet partitions min E ε(u) s.t. u j L 2 (Ω) = 1. u H1(Ω;R k ) The nearest-point projection map, Π T : R k T, for T = Σ k is given by { x (Π T x) i = i x i = max j x j 0 otherwise Diffusion generated method. For i = 1, 2,..., Step 1. Solve the diffusion equation until time τ t φ = φ, φ(x, t = 0) = u i (x). Step 2. Point-wise, apply the nearest-point projection map: ũ i+1 (x) = Π T φ(x, τ). Step 3. Nomarlize: u i+1 (x) = ũi+1(x) ũ i+1 L 2 (Ω).
24 Results for 2D flat tori, k = 3 9, 11, 12, 15, 16, and 20 23
25 Results for 3D flat tori, k = 2 24
26 Results for 3D flat tori, k = 4, tessellation by rhombic dodecahedra 25
27 Results for 3D flat tori, k = 8, Weaire-Phelan structure 26
28 Results for 3D flat tori, k = 12, Kelvin s structure composed of truncated octahedra 27
29 Results for 4D flat tori, k = 8, 24-cell honeycomb 28 Dong Wang Generalized MBO May. 3, 2018 IMS, NUS, Singapore
30 Discussion and future directions for generalized MBO methods 29 We only considered a single matrix valued field that has two "phases given by when the determinant is positive or negative. It would be very interesting to extend this work to the multi-phase problem as was accomplished for n = 1 in [Esedoglu+Otto, 2015]. For O(n) valued fields with n = 2, the motion law for the interface is unknown. For n = 2 on a two-dimensional flat torus, further analysis regarding the winding number of the field is required. Is it possible to determine the final solution based on the winding number of the initial field? For problems with a non-trivial boundary condition, it not obvious how to adapt the Lyapunov functional. Thanks! Questions? dwang@math.utah.edu B. Osting and D. Wang, A generalized MBO diffusion generated motion for orthogonal matrix valued fields, submitted, arxiv: (2017). D. Wang and B. Osting, A diffusion generated method for computing Dirichlet partitions, submitted, arxiv: (2018).
An improved threshold dynamics method for wetting dynamics
An improved threshold dynamics method for wetting dynamics Dong Wang a, Xiao-Ping Wang b,, Xianmin Xu c,d a Department of Mathematics, University of Utah, Salt Lake City, UT, 8, USA. b Department of Mathematics,
More informationarxiv: v1 [math.oc] 22 Dec 2018
AN EFFICIENT THRESHOLD DYNAMICS METHOD FOR TOPOLOGY OPTIMIZATION FOR FLUIDS HUANGXIN CHEN, HAITAO LENG, DONG WANG, AND XIAO-PING WANG arxiv:1812.09437v1 [math.oc] 22 Dec 2018 Abstract. We propose an efficient
More informationDiffuse interface methods on graphs: Data clustering and Gamma-limits
Diffuse interface methods on graphs: Data clustering and Gamma-limits Yves van Gennip joint work with Andrea Bertozzi, Jeff Brantingham, Blake Hunter Department of Mathematics, UCLA Research made possible
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 informationBranched transport limit of the Ginzburg-Landau functional
Branched transport limit of the Ginzburg-Landau functional Michael Goldman CNRS, LJLL, Paris 7 Joint work with S. Conti, F. Otto and S. Serfaty Introduction Superconductivity was first observed by Onnes
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 informationESTIMATES FOR THE MONGE-AMPERE EQUATION
GLOBAL W 2,p ESTIMATES FOR THE MONGE-AMPERE EQUATION O. SAVIN Abstract. We use a localization property of boundary sections for solutions to the Monge-Ampere equation obtain global W 2,p estimates under
More informationSmall energy regularity for a fractional Ginzburg-Landau system
Small energy regularity for a fractional Ginzburg-Landau system Yannick Sire University Aix-Marseille Work in progress with Vincent Millot (Univ. Paris 7) The fractional Ginzburg-Landau system We are interest
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 information(Non)local phase transitions and minimal surfaces
Dipartimento di Matematica valdinoci@mat.uniroma2.it Italian-Japanese workshop Cortona, June 2011 Outline of the talk Outline of the talk ( ) s u = F 1( ξ 2s (Fu) ), where s (0, 1) and F is the Fourier
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 informationDeforming conformal metrics with negative Bakry-Émery Ricci Tensor on manifolds with boundary
Deforming conformal metrics with negative Bakry-Émery Ricci Tensor on manifolds with boundary Weimin Sheng (Joint with Li-Xia Yuan) Zhejiang University IMS, NUS, 8-12 Dec 2014 1 / 50 Outline 1 Prescribing
More informationALGORITHMS FOR THE EVOLUTION OF SURFACES
ALGORITHMS FOR THE EVOLUTION OF SURFACES MATTHEW STONE Abstract. We develop algorithms for simulating various curvature dependent motions of interfaces. These motions are described by a class of partial
More informationVISCOSITY SOLUTIONS. We follow Han and Lin, Elliptic Partial Differential Equations, 5.
VISCOSITY SOLUTIONS PETER HINTZ We follow Han and Lin, Elliptic Partial Differential Equations, 5. 1. Motivation Throughout, we will assume that Ω R n is a bounded and connected domain and that a ij C(Ω)
More 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 informationPartial regularity for fully nonlinear PDE
Partial regularity for fully nonlinear PDE Luis Silvestre University of Chicago Joint work with Scott Armstrong and Charles Smart Outline Introduction Intro Review of fully nonlinear elliptic PDE Our result
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More informationarxiv: v1 [math.ap] 18 Jan 2019
manuscripta mathematica manuscript No. (will be inserted by the editor) Yongpan Huang Dongsheng Li Kai Zhang Pointwise Boundary Differentiability of Solutions of Elliptic Equations Received: date / Revised
More informationDegenerate Monge-Ampère equations and the smoothness of the eigenfunction
Degenerate Monge-Ampère equations and the smoothness of the eigenfunction Ovidiu Savin Columbia University November 2nd, 2015 Ovidiu Savin (Columbia University) Degenerate Monge-Ampère equations November
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 information(a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? Solution: dim N(A) 1, since rank(a) 3. Ax =
. (5 points) (a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? dim N(A), since rank(a) 3. (b) If we also know that Ax = has no solution, what do we know about the rank of A? C(A)
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 informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationThe Convergence of Mimetic Discretization
The Convergence of Mimetic Discretization for Rough Grids James M. Hyman Los Alamos National Laboratory T-7, MS-B84 Los Alamos NM 87545 and Stanly Steinberg Department of Mathematics and Statistics University
More informationA sharp diffuse interface tracking method for approximating evolving interfaces
A sharp diffuse interface tracking method for approximating evolving interfaces Vanessa Styles and Charlie Elliott University of Sussex Overview Introduction Phase field models Double well and double obstacle
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 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 informationFoliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary
Foliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary David Chopp and John A. Velling December 1, 2003 Abstract Let γ be a Jordan curve in S 2, considered as the ideal
More informationPinning and depinning of interfaces in random media
Pinning and depinning of interfaces in random media Patrick Dondl joint work with Nicolas Dirr and Michael Scheutzow March 17, 2011 at Université d Orléans signal under crossed polarizers is measured and
More informationNumerical Solutions of Geometric Partial Differential Equations. Adam Oberman McGill University
Numerical Solutions of Geometric Partial Differential Equations Adam Oberman McGill University Sample Equations and Schemes Fully Nonlinear Pucci Equation! "#+ "#* "#) "#( "#$ "#' "#& "#% "#! "!!!"#$ "
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
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 information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationNon-trivial pinning threshold for an evolution equation involving long range interactions
Non-trivial pinning threshold for an evolution equation involving long range interactions Martin Jesenko joint work with Patrick Dondl Workshop: New trends in the variational modeling of failure phenomena
More informationDef. The euclidian distance between two points x = (x 1,...,x p ) t and y = (y 1,...,y p ) t in the p-dimensional space R p is defined as
MAHALANOBIS DISTANCE Def. The euclidian distance between two points x = (x 1,...,x p ) t and y = (y 1,...,y p ) t in the p-dimensional space R p is defined as d E (x, y) = (x 1 y 1 ) 2 + +(x p y p ) 2
More informationDiscrete Euclidean Curvature Flows
Discrete Euclidean Curvature Flows 1 1 Department of Computer Science SUNY at Stony Brook Tsinghua University 2010 Isothermal Coordinates Relation between conformal structure and Riemannian metric Isothermal
More informationMinimal k-partition for the p-norm of the eigenvalues
Minimal k-partition for the p-norm of the eigenvalues V. Bonnaillie-Noël DMA, CNRS, ENS Paris joint work with B. Bogosel, B. Helffer, C. Léna, G. Vial Calculus of variations, optimal transportation, and
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 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 informationMeasuring Ellipsoids 1
Measuring Ellipsoids 1 Igor Rivin Temple University 2 What is an ellipsoid? E = {x E n Ax = 1}, where A is a non-singular linear transformation of E n. Remark that Ax = Ax, Ax = x, A t Ax. The matrix Q
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 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 informationCoupled second order singular perturbations for phase transitions
Coupled second order singular perturbations for phase transitions CMU 06/09/11 Ana Cristina Barroso, Margarida Baía, Milena Chermisi, JM Introduction Let Ω R d with Lipschitz boundary ( container ) and
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
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 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 informationA Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term
A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term Peter Sternberg In collaboration with Dmitry Golovaty (Akron) and Raghav Venkatraman (Indiana) Department of Mathematics
More informationLinear algebra for computational statistics
University of Seoul May 3, 2018 Vector and Matrix Notation Denote 2-dimensional data array (n p matrix) by X. Denote the element in the ith row and the jth column of X by x ij or (X) ij. Denote by X j
More information4 Divergence theorem and its consequences
Tel Aviv University, 205/6 Analysis-IV 65 4 Divergence theorem and its consequences 4a Divergence and flux................. 65 4b Piecewise smooth case............... 67 4c Divergence of gradient: Laplacian........
More informationOptimization and Optimal Control in Banach Spaces
Optimization and Optimal Control in Banach Spaces Bernhard Schmitzer October 19, 2017 1 Convex non-smooth optimization with proximal operators Remark 1.1 (Motivation). Convex optimization: easier to solve,
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 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 informationIterative Reweighted Minimization Methods for l p Regularized Unconstrained Nonlinear Programming
Iterative Reweighted Minimization Methods for l p Regularized Unconstrained Nonlinear Programming Zhaosong Lu October 5, 2012 (Revised: June 3, 2013; September 17, 2013) Abstract In this paper we study
More informationMeasure and Integration: Solutions of CW2
Measure and Integration: s of CW2 Fall 206 [G. Holzegel] December 9, 206 Problem of Sheet 5 a) Left (f n ) and (g n ) be sequences of integrable functions with f n (x) f (x) and g n (x) g (x) for almost
More informationDECOUPLING LECTURE 6
18.118 DECOUPLING LECTURE 6 INSTRUCTOR: LARRY GUTH TRANSCRIBED BY DONGHAO WANG We begin by recalling basic settings of multi-linear restriction problem. Suppose Σ i,, Σ n are some C 2 hyper-surfaces in
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 informationSymmetric Matrices and Eigendecomposition
Symmetric Matrices and Eigendecomposition Robert M. Freund January, 2014 c 2014 Massachusetts Institute of Technology. All rights reserved. 1 2 1 Symmetric Matrices and Convexity of Quadratic Functions
More informationRandomized Coordinate Descent Methods on Optimization Problems with Linearly Coupled Constraints
Randomized Coordinate Descent Methods on Optimization Problems with Linearly Coupled Constraints By I. Necoara, Y. Nesterov, and F. Glineur Lijun Xu Optimization Group Meeting November 27, 2012 Outline
More informationRichard F. Bass Krzysztof Burdzy University of Washington
ON DOMAIN MONOTONICITY OF THE NEUMANN HEAT KERNEL Richard F. Bass Krzysztof Burdzy University of Washington Abstract. Some examples are given of convex domains for which domain monotonicity of the Neumann
More informationMEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW
MEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW GREGORY DRUGAN AND XUAN HIEN NGUYEN Abstract. We present two initial graphs over the entire R n, n 2 for which the mean curvature flow
More informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
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 informationAlgorithms for Nonsmooth Optimization
Algorithms for Nonsmooth Optimization Frank E. Curtis, Lehigh University presented at Center for Optimization and Statistical Learning, Northwestern University 2 March 2018 Algorithms for Nonsmooth Optimization
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
More informationAlgorithms for Multiphase Partitioning
Algorithms for Multiphase Partitioning by Matthew Jacobs A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Mathematics) in the University of Michigan
More informationBackground Mathematics (2/2) 1. David Barber
Background Mathematics (2/2) 1 David Barber University College London Modified by Samson Cheung (sccheung@ieee.org) 1 These slides accompany the book Bayesian Reasoning and Machine Learning. The book and
More informationLearning gradients: prescriptive models
Department of Statistical Science Institute for Genome Sciences & Policy Department of Computer Science Duke University May 11, 2007 Relevant papers Learning Coordinate Covariances via Gradients. Sayan
More informationCentre for Mathematics and Its Applications The Australian National University Canberra, ACT 0200 Australia. 1. Introduction
ON LOCALLY CONVEX HYPERSURFACES WITH BOUNDARY Neil S. Trudinger Xu-Jia Wang Centre for Mathematics and Its Applications The Australian National University Canberra, ACT 0200 Australia Abstract. In this
More informationReconstructing inclusions from Electrostatic Data
Reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: W. Rundell Purdue
More informationConvexity of level sets for solutions to nonlinear elliptic problems in convex rings. Paola Cuoghi and Paolo Salani
Convexity of level sets for solutions to nonlinear elliptic problems in convex rings Paola Cuoghi and Paolo Salani Dip.to di Matematica U. Dini - Firenze - Italy 1 Let u be a solution of a Dirichlet problem
More informationLecture 2: Isoperimetric methods for the curve-shortening flow and for the Ricci flow on surfaces
Lecture 2: Isoperimetric methods for the curve-shortening flow and for the Ricci flow on surfaces Ben Andrews Mathematical Sciences Institute, Australian National University Winter School of Geometric
More information2. Matrix Algebra and Random Vectors
2. Matrix Algebra and Random Vectors 2.1 Introduction Multivariate data can be conveniently display as array of numbers. In general, a rectangular array of numbers with, for instance, n rows and p columns
More informationDemixing in viscous fluids: Connection with OT. Felix Otto. joint work with: R. V. Kohn, Y. Brenier, C. Seis, D. Slepcev
Demixing in viscous fluids: Connection with OT Felix Otto Max Planck Institute for Mathematics in the Sciences Leipzig, Germany joint work with: R. V. Kohn, Y. Brenier, C. Seis, D. Slepcev Geometric evolution
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationConvexity and Non-Convexity in Partitioning and Interpolation Problems
Convexity and Non-Convexity in Partitioning and Interpolation Problems Jan Lellmann CIA/DAMTP, University of Cambridge Joint work with: C. Schnörr, F. Lenzen (IPA/HCI, University of Heidelberg), F. Widmann,
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 informationPotential Analysis meets Geometric Measure Theory
Potential Analysis meets Geometric Measure Theory T. Toro Abstract A central question in Potential Theory is the extend to which the geometry of a domain influences the boundary regularity of the solution
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
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 informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More information(x, y) = d(x, y) = x y.
1 Euclidean geometry 1.1 Euclidean space Our story begins with a geometry which will be familiar to all readers, namely the geometry of Euclidean space. In this first chapter we study the Euclidean distance
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 information1.The anisotropic plate model
Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 6 OPTIMAL CONTROL OF PIEZOELECTRIC ANISOTROPIC PLATES ISABEL NARRA FIGUEIREDO AND GEORG STADLER ABSTRACT: This paper
More informationWave equation on manifolds and finite speed of propagation
Wave equation on manifolds and finite speed of propagation Ethan Y. Jaffe Let M be a Riemannian manifold (without boundary), and let be the (negative of) the Laplace-Beltrami operator. In this note, we
More informationRegularity estimates for the solution and the free boundary of the obstacle problem for the fractional Laplacian
Regularity estimates for the solution and the free boundary of the obstacle problem for the fractional Laplacian Luis Caffarelli, Sandro Salsa and Luis Silvestre October 15, 2007 Abstract We use a characterization
More informationarxiv: v1 [math.ap] 10 Apr 2013
QUASI-STATIC EVOLUTION AND CONGESTED CROWD TRANSPORT DAMON ALEXANDER, INWON KIM, AND YAO YAO arxiv:1304.3072v1 [math.ap] 10 Apr 2013 Abstract. We consider the relationship between Hele-Shaw evolution with
More informationConvex solutions to the mean curvature flow
Annals of Mathematics 173 (2011), 1185 1239 doi: 10.4007/annals.2011.173.3.1 Convex solutions to the mean curvature flow By Xu-Jia Wang Abstract In this paper we study the classification of ancient convex
More informationOct : Lecture 15: Surface Integrals and Some Related Theorems
90 MIT 3.016 Fall 2005 c W.C Carter Lecture 15 Oct. 21 2005: Lecture 15: Surface Integrals and Some Related Theorems Reading: Kreyszig Sections: 9.4 pp:485 90), 9.5 pp:491 495) 9.6 pp:496 505) 9.7 pp:505
More informationAdaptive evolution : a population approach Benoît Perthame
Adaptive evolution : a population approach Benoît Perthame Adaptive dynamic : selection principle d dt n(x, t) = n(x, t)r( x, ϱ(t) ), ϱ(t) = R d n(x, t)dx. given x, n(x) = ϱ δ(x x), R( x, ϱ) = 0, ϱ( x).
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 informationSpectral Clustering. Spectral Clustering? Two Moons Data. Spectral Clustering Algorithm: Bipartioning. Spectral methods
Spectral Clustering Seungjin Choi Department of Computer Science POSTECH, Korea seungjin@postech.ac.kr 1 Spectral methods Spectral Clustering? Methods using eigenvectors of some matrices Involve eigen-decomposition
More informationEE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 17
EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 17 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory May 29, 2012 Andre Tkacenko
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 informationTHE RADON NIKODYM PROPERTY AND LIPSCHITZ EMBEDDINGS
THE RADON NIKODYM PROPERTY AND LIPSCHITZ EMBEDDINGS Motivation The idea here is simple. Suppose we have a Lipschitz homeomorphism f : X Y where X and Y are Banach spaces, namely c 1 x y f (x) f (y) c 2
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 informationThe Steepest Descent Algorithm for Unconstrained Optimization
The Steepest Descent Algorithm for Unconstrained Optimization Robert M. Freund February, 2014 c 2014 Massachusetts Institute of Technology. All rights reserved. 1 1 Steepest Descent Algorithm The problem
More information