scholarly journals Motion by anisotropic mean curvature as sharp interface limit of an inhomogeneous and anisotropic Allen–Cahn equation

Author(s):  
Matthieu Alfaro ◽  
Harald Garcke ◽  
Danielle Hilhorst ◽  
Hiroshi Matano ◽  
Reiner Schätzle

We consider the spatially inhomogeneous and anisotropic reaction–diffusion equation ut = m(x)−1 div[m(x)ap(x,∇u)] + ε−2f(u), involving a small parameter ε > 0 and a bistable nonlinear term whose stable equilibria are 0 and 1. We use a Finsler metric related to the anisotropic diffusion term and work in relative geometry. We prove a weak comparison principle and perform an analysis of both the generation and the motion of interfaces. More precisely, we show that, within the time-scale of order ε2|ln ε|, the unique weak solution uε develops a steep transition layer that separates the regions {uε ≈ 0} and {uε | 1}. Then, on a much slower time-scale, the layer starts to propagate. Consequently, as ε → 0, the solution uε converges almost everywhere (a.e.) to 0 in Ω−t and 1 in Ω+t , where Ω−t and Ω+t are sub-domains of Ω separated by an interface Гt, whose motion is driven by its anisotropic mean curvature. We also prove that the thickness of the transition layer is of order ε.

1998 ◽  
Vol 08 (04) ◽  
pp. 573-601 ◽  
Author(s):  
FRANCESCA FIERRO ◽  
ROBERTA GOGLIONE ◽  
MAURIZIO PAOLINI

In this paper we present and discuss the results of some numerical simulations in order to investigate the mean curvature flow problem in the presence of a nonconvex anisotropy. Mathematically, nonconvexity of the anisotropy leads to the ill-posedness of the evolution problem, which becomes forward–backward parabolic. Simulations presented here refer to two different settings: curvature driven vertical motion of graphs (nonparametric setting) and motion in the normal direction by anisotropic mean curvature of surfaces (parametric setting). In the latter we first relax the problem via an Allen–Cahn type reaction-diffusion equation, in the context of Finsler geometry (diffused interface approximation). Our results suggest three main points. A nonconvex anisotropy and its convexification give rise, for both settings and the discretizations considered, to different evolutions. Wrinkled regions seem to appear only in correspondence to locally concave parts of the anisotropy. Moreover, locally convex regions (interior to the convexification of the so-called Frank diagram) seem to play an important role.


2018 ◽  
Vol 41 (18) ◽  
pp. 9203-9217 ◽  
Author(s):  
Natalia T. Levashova ◽  
Nikolay N. Nefedov ◽  
Olga A. Nikolaeva ◽  
Andrey O. Orlov ◽  
Alexander A. Panin

Mathematics ◽  
2020 ◽  
Vol 8 (12) ◽  
pp. 2111
Author(s):  
Mitică CRAUS ◽  
Silviu-Dumitru PAVĂL

In this paper we propose and compare two methods to optimize the numerical computations for the diffusion term in a nonlocal formulation for a reaction-diffusion equation. The diffusion term is particularly computationally intensive due to the integral formulation, and thus finding a better way of computing its numerical approximation could be of interest, given that the numerical analysis usually takes place on large input domains having more than one dimension. After introducing the general reaction-diffusion model, we discuss a numerical approximation scheme for the diffusion term, based on a finite difference method. In the next sections we propose two algorithms to solve the numerical approximation scheme, focusing on finding a way to improve the time performance. While the first algorithm (sequential) is used as a baseline for performance measurement, the second algorithm (parallel) is implemented using two different memory-sharing parallelization technologies: Open Multi-Processing (OpenMP) and CUDA. All the results were obtained by using the model in image processing applications such as image restoration and segmentation.


2012 ◽  
Vol 22 (11) ◽  
pp. 1250285
Author(s):  
CHANGBING TANG ◽  
FANGYUE CHEN ◽  
JIANBO WANG ◽  
XIANG LI

The fundamental local passivity theory asserts that a wide spectrum of complex behaviors may exist if the cells in the reaction–diffusion are not locally passive. This local passivity principle has provided a powerful tool for studying the complexity in a homogeneous lattice formed by coupled cells. In this paper, the complexity matrix YQ(s), which is the tool for testing the local passivity theory, is modified based on the characteristic polynomial AQ(λ). Then, the local passivity theory is applied to the study of the Oregonator CNN to judge if the cell parameters of a CNN are chosen at the edge of chaos. Analysis of the bifurcation and the numerical simulations show that nonzero diffusion term in Oregonator CNN may cause a reaction–diffusion equation oscillating under the appropriate choice of diffusion coefficient if the local passivity theory is not satisfied. That is, if the cell parameters of a CNN are chosen at the edge of chaos, the system is potentially unstable.


1996 ◽  
Vol 7 (3) ◽  
pp. 287-301 ◽  
Author(s):  
J. W. Cahn ◽  
C. M. Elliott ◽  
A. Novick-Cohen

We show by using formal asymptotics that the zero level set of the solution to the Cahn–Hilliard equation with a concentration dependent mobility approximates to lowest order in ɛ. an interface evolving according to the geometric motion,(where V is the normal velocity, Δ8 is the surface Laplacian and κ is the mean curvature of the interface), both in the deep quench limit and when the temperature θ is where є2 is the coefficient of gradient energy. Equation (0.1) may be viewed as motion by surface diffusion, and as a higher-order analogue of motion by mean curvature predicted by the bistable reaction-diffusion equation.


Sign in / Sign up

Export Citation Format

Share Document