• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.1
• /
• pp.27-41
• /
• 2014

In this paper, we present an efficient numerical method for multiphase image segmentation using a multiphase-field model. The method combines the vector-valued Allen-Cahn phase-field equation with initial data fitting terms containing prescribed interface width and fidelity constants. An efficient numerical solution is achieved using the recently developed hybrid operator splitting method for the vector-valued Allen-Cahn phase-field equation. We split the modified vector-valued Allen-Cahn equation into a nonlinear equation and a linear diffusion equation with a source term. The linear diffusion equation is discretized using an implicit scheme and the resulting implicit discrete system of equations is solved by a multigrid method. The nonlinear equation is solved semi-analytically using a closed-form solution. And by treating the source term of the linear diffusion equation explicitly, we solve the modified vector-valued Allen-Cahn equation in a decoupled way. By decoupling the governing equation, we can speed up the segmentation process with multiple phases. We perform some characteristic numerical experiments for multiphase image segmentation.

• Journal of applied mathematics & informatics
• /
• v.29 no.5_6
• /
• pp.1477-1487
• /
• 2011

We propose a numerical method for solving Allen-Cahn equation, in both one-dimensional and two-dimensional cases. The new scheme that is explicit, stable, and easy to compute is obtained and the proposed method provides a straightforward and effective way for nonlinear evolution equations.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.1
• /
• pp.1-16
• /
• 2017

The Allen-Cahn equation is solved numerically by operator splitting Fourier spectral methods. The basic idea of the operator splitting method is to decompose the original problem into sub-equations and compose the approximate solution of the original equation using the solutions of the subproblems. The purpose of this paper is to characterize higher order operator splitting schemes and propose several higher order methods. Unlike the first and the second order methods, each of the heat and the free-energy evolution operators has at least one backward evaluation in higher order methods. We investigate the effect of negative time steps on a general form of third order schemes and suggest three third order methods for better stability and accuracy. Two fourth order methods are also presented. The traveling wave solution and a spinodal decomposition problem are used to demonstrate numerical properties and the order of convergence of the proposed methods.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.14 no.4
• /
• pp.201-210
• /
• 2010

We propose a new robust and accurate method for the numerical solution of medical image segmentation. The modified Allen-Cahn equation is used to model the boundaries of the image regions. Its numerical algorithm is based on operator splitting techniques. In the first step of the splitting scheme, we implicitly solve the heat equation with the variable diffusive coefficient and a source term. Then, in the second step, using a closed-form solution for the nonlinear equation, we get an analytic solution. We overcome the time step constraint associated with most numerical implementations of geometric active contours. We demonstrate performance of the proposed image segmentation algorithm on several artificial as well as real image examples.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.16 no.1
• /
• pp.15-29
• /
• 2012

In this paper we introduce 6-fold symmetry crystal growth using new phase-field models based on the modified Allen-Cahn equation. The proposed method is a hybrid method which uses both analytic and numerical solutions. We then show this method can be extended to $k$-fold case. The Wulff construction procedure is provided to understand and predict the shape of crystals. We also present a detailed mathematical proof of the validity of the Wulff construction. For computational results, we verify the accuracy and efficiency of the method for snow crystal growth.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2014.10a
• /
• pp.555-556
• /
• 2014

In this paper, topology optimization of two-dimensional acoustic lenses is presented by using the phase field method. The objective of the optimization is to maximize the acoustic pressure at a specified domain inside the acoustic domain for a given frequency, and the constraint is imposed on the amount of the material of the acoustic lens. Topology optimization of two-dimensional acoustic lenses are obtained as the steady state of the phase transition described by the Allen-Cahn equation. The Helmholtz equation modeling the wave propagation is solved by using a finite element method. The effectiveness of the proposed method is verified by applying it for several two-dimensional acoustic lens system design problems.

• Korean Journal of Mathematics
• /
• v.27 no.1
• /
• pp.193-219
• /
• 2019

The Navier-Stokes equations with variable density are challenging problems in numerical analysis community. We recently built the 2nd order stabilized Gauge-Uzawa method [SGUM] to solve the Navier-Stokes equations with constant density and have estimated theoretically optimal accuracy. Also we proved that SGUM is unconditionally stable. In this paper, we apply SGUM to the Navier-Stokes equations with nonconstant variable density and find out the stability condition of the algorithms. Because the condition is rather strong to apply to real problems, we consider Allen-Cahn scheme to construct unconditionally stable scheme.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.26 no.4
• /
• pp.333-342
• /
• 2022

In this paper, we perform a direct comparison study of real experimental data for domain rearrangement and the Cahn-Hilliard (CH) equation on the dynamics of morphological evolution. To validate a mathematical model for physical phenomena, we take initial conditions from experimental images by using an image segmentation technique. The image segmentation algorithm is based on the Mumford-Shah functional and the Allen-Cahn (AC) equation. The segmented phase-field profile is similar to the solution of the CH equation, that is, it has hyperbolic tangent profile across interfacial transition region. We use unconditionally stable schemes to solve the governing equations. As a test problem, we take domain rearrangement of lipid bilayers. Numerical results demonstrate that comparison of the evolutions with experimental data is a good benchmark test for validating a mathematical model.

• Nuclear Engineering and Technology
• /
• v.54 no.1
• /
• pp.226-233
• /
• 2022

In the present work, a phase-field model was developed to investigate the influence of recrystallization on bubble evolution during irradiation. Considering the interaction between bubbles and grain boundary (GB), a set of modified Cahn-Hilliard and Allen-Cahn equations, with field variables and order parameters evolving in space and time, was used in this model. Both the kinetics of recrystallization characterized in experiments and point defects generated during cascade were incorporated in the model. The bubble evolution in recrystallized polycrystalline of U-Mo alloy was also investigated. The simulation results showed that GB with a large area fraction generated by recrystallization accelerates the formation and growth of bubbles. With the formation of new grains, gas atoms are swept and collected by GBs. The simulation results of bubble size and distribution are consistent with the experimental results.