538 research outputs found
Multiclass Data Segmentation using Diffuse Interface Methods on Graphs
We present two graph-based algorithms for multiclass segmentation of
high-dimensional data. The algorithms use a diffuse interface model based on
the Ginzburg-Landau functional, related to total variation compressed sensing
and image processing. A multiclass extension is introduced using the Gibbs
simplex, with the functional's double-well potential modified to handle the
multiclass case. The first algorithm minimizes the functional using a convex
splitting numerical scheme. The second algorithm is a uses a graph adaptation
of the classical numerical Merriman-Bence-Osher (MBO) scheme, which alternates
between diffusion and thresholding. We demonstrate the performance of both
algorithms experimentally on synthetic data, grayscale and color images, and
several benchmark data sets such as MNIST, COIL and WebKB. We also make use of
fast numerical solvers for finding the eigenvectors and eigenvalues of the
graph Laplacian, and take advantage of the sparsity of the matrix. Experiments
indicate that the results are competitive with or better than the current
state-of-the-art multiclass segmentation algorithms.Comment: 14 page
Extended Smoothed Boundary Method for Solving Partial Differential Equations with General Boundary Conditions on Complex Boundaries
In this article, we describe an approach for solving partial differential
equations with general boundary conditions imposed on arbitrarily shaped
boundaries. A continuous function, the domain parameter, is used to modify the
original differential equations such that the equations are solved in the
region where a domain parameter takes a specified value while boundary
conditions are imposed on the region where the value of the domain parameter
varies smoothly across a short distance. The mathematical derivations are
straightforward and generically applicable to a wide variety of partial
differential equations. To demonstrate the general applicability of the
approach, we provide four examples herein: (1) the diffusion equation with both
Neumann and Dirichlet boundary conditions; (2) the diffusion equation with both
surface diffusion and reaction; (3) the mechanical equilibrium equation; and
(4) the equation for phase transformation with the presence of additional
boundaries. The solutions for several of these cases are validated against
corresponding analytical and semi-analytical solutions. The potential of the
approach is demonstrated with five applications: surface-reaction-diffusion
kinetics with a complex geometry, Kirkendall-effect-induced deformation,
thermal stress in a complex geometry, phase transformations affected by
substrate surfaces, and a self-propelled droplet.Comment: This document is the revised version of arXiv:0912.1288v
Enhancement of Krylov Subspace Spectral Methods Through the Use of the Residual
Depending on the type of equation, finding the solution of a time-dependent partial differential equation can be quite challenging. Although modern time-stepping methods for solving these equations have become more accurate for a small number of grid points, in a lot of cases the scalability of those methods leaves much to be desired. That is, unless the timestep is chosen to be sufficiently small, the computed solutions might exhibit unreasonable behavior with large input sizes. Therefore, to improve accuracy as the number of grid points increases, the time-steps must be chosen to be even smaller to reach a reasonable solution.
Krylov subspace spectral (KSS) methods are componentwise, scalable, methods used to solve time-dependent, variable coefficient partial differential equations. The main idea behind KSS methods is to use an interpolating polynomial with frequency dependent interpolation points to approximate a solution operator for each Fourier coefficient.
This dissertation will discuss two techniques that were developed to eliminate error in the low frequency components of the solution computed using KSS methods. These two methods are a multigrid inspired technique (coarse grid residual correction) and developing a step size controller using the residual as an error approximation (adaptive time stepping)
Characterizing the stabilization size for semi-implicit Fourier-spectral method to phase field equations
Recent results in the literature provide computational evidence that
stabilized semi-implicit time-stepping method can efficiently simulate phase
field problems involving fourth-order nonlinear dif- fusion, with typical
examples like the Cahn-Hilliard equation and the thin film type equation. The
up-to-date theoretical explanation of the numerical stability relies on the
assumption that the deriva- tive of the nonlinear potential function satisfies
a Lipschitz type condition, which in a rigorous sense, implies the boundedness
of the numerical solution. In this work we remove the Lipschitz assumption on
the nonlinearity and prove unconditional energy stability for the stabilized
semi-implicit time-stepping methods. It is shown that the size of stabilization
term depends on the initial energy and the perturba- tion parameter but is
independent of the time step. The corresponding error analysis is also
established under minimal nonlinearity and regularity assumptions
Fast and Stable Schemes for Phase Fields Models
We propose and analyse new stabilized time marching schemes for Phase Fields model such as Allen-Cahn and Cahn-Hillard equations, when discretized in space with high order finite differences compact schemes. The stabilization applies to semi-implicit schemes for which the linear part is simplified using sparse pre-conditioners. The new methods allow to significant obtain a gain of CPU time. The numerical illustrations we give concern applications on pattern dynamics and on image processing (inpainting, segmentation) in two and three dimension cases
Krylov Subspace Spectral Method with Multigrid for a Time-Dependent, Variable-Coefficient Partial Differential Equation
Krylov Subspace Spectral (KSS) methods are traditionally used to solve time-dependent, variable-coefficient PDEs. They are high-order accurate, component-wise methods that are efficient with variable input sizes.
This thesis will demonstrate how one can make KSS methods even more efficient by using a Multigrid-like approach for low-frequency components. The essential ingredients of Multigrid, such as restriction, residual correction, and prolongation, are adapted to the timedependent case. Then a comparison of KSS, KSS with Multigrid, KSS-EPI and standard Krylov projection methods will be demonstrated
Primal-dual active set methods for Allen-Cahn variational inequalities
This thesis aims to introduce and analyse a primal-dual active set strategy for solving Allen-Cahn variational inequalities. We consider the standard Allen-Cahn
equation with non-local constraints and a vector-valued Allen-Cahn equation with and without non-local constraints. Existence and uniqueness results are derived
in a formulation involving Lagrange multipliers for local and non-local constraints. Local Convergence is shown by interpreting the primal-dual active set approach as
a semi-smooth Newton method. Properties of the method are discussed and several numerical simulations in two and three space dimensions demonstrate its efficiency.
In the second part of the thesis various applications of the Allen-Cahn equation are discussed. The non-local Allen-Cahn equation can be coupled with an elasticity
equation to solve problems in structural topology optimisation. The model can be extended to handle multiple structures by using the vector-valued Allen-Cahn
variational inequality with non-local constraints. Since many applications of the Allen-Cahn equation involve evolution of interfaces in materials an important extension of the standard Allen-Cahn model is to allow materials to exhibit anisotropic behaviour. We introduce an anisotropic version of the Allen-Cahn variational inequality and we show that it is possible to apply the primal-dual active set strategy efficiently to this model. Finally, the Allen-Cahn model is applied to problems in image processing, such as segmentation, denoising and inpainting.
The primal-dual active set method proves exible and reliable for all the applications considered in this thesis
- …