# American Institute of Mathematical Sciences

• Previous Article
Numerical and mathematical analysis of blow-up problems for a stochastic differential equation
• DCDS-S Home
• This Issue
• Next Article
A method how to determine parameters arising in a smoldering evolution equation by image segmentation for experiment's movies
March  2021, 14(3): 893-907. doi: 10.3934/dcdss.2020390

## Numerical analysis of an ODE and a level set methods for evolving spirals by crystalline eikonal-curvature flow

 1 Department of Mathematical Sciences, Shibaura Institute of Technology, Fukasaku 309, Minuma-ku, Saitama 337-8570, Japan 2 Division of Pure and Applied Science, Faculty of Science and Technology, Gunma University, Aramaki-machi 4-2, Maebashi, 371-8510 Gunma, Japan

* Corresponding author: Takeshi Ohtsuka

Received  January 2019 Revised  February 2020 Published  June 2020

In this paper, the evolution of a polygonal spiral curve by the crystalline curvature flow with a pinned center is considered from two viewpoints; a discrete model consisting of an ODE system describing facet lengths and another using level set method. We investigate the difference of these models numerically by calculating the area of an interposed region by their spiral curves. The area difference is calculated by the normalized $L^1$ norm of the difference of step-like functions which are branches of $\arg (x)$ whose discontinuities are on the spirals. We find that the differences in the numerical results are small, even though the model equations around the center and the farthest facet are slightly different.

Citation: Tetsuya Ishiwata, Takeshi Ohtsuka. Numerical analysis of an ODE and a level set methods for evolving spirals by crystalline eikonal-curvature flow. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 893-907. doi: 10.3934/dcdss.2020390
##### References:

show all references

##### References:
A figure of two spirals (the solid and dashed lines) and the interposed region by them. The function $\mathcal{D} (t)$ indicates the area of the gray regions
Description of $\Gamma_D = \bigcup_{j = 0}^k L_j (t)$. Note that, for the simplicity, the variable $t$ of $L_j$ and $y_j$ is omitted in the above figure
Construction of $\theta_D (t, x)$; we construct a branch of $\arg (x)$ whose discontinuities are only on $\Gamma (t)$(the dashed line in (1)). For this purpose we first construct $\vartheta (x) = \arg (x)$ whose discontinuities are only on $\mathcal{L}_k (t)$ (the solid line in (2)). Then, we make go down the height of $\vartheta (x)$ on $R_{j} (t)$ (the gray region in (3) or (4)) with the jump-height $2 \pi$ from $j = k-1$ to $j = 0$ inductively to remove illegal discontinuities. The solid line in figure (3) or (4) denotes the discontinuity of $\Theta_{k,k-1}$ or $\Theta_{k,k-2}$, respectively
Profiles of the square spiral at $t = 1$. The level set method is calculated using $\rho = 0.02$ and $\Delta x = 0.0050$
Graphs of functions $\mathcal{D} (t)$ for the square spiral with a fixed center radius $\rho = 0.02$(left), and with a reduced center radius $\rho = 2 \Delta x$(right)
Profiles of the diagonal spiral at $t = 1$. The level set method is calculated using $\rho = 0.02$ and $\Delta x = 0.0050$
Graphs of $\mathcal{D}(t)$ for the diagonal spiral with a fixed center radius $\rho = 0.02$(left), and with a reduced center radius $\rho = 4 \Delta x$(right)
Profiles of the triangle spiral at $t = 0.8$. The level set method is calculated using $\rho = 0.02$ and $\Delta x = 0.0050$
Graphs of $\mathcal{D} (t)$ for the triangle spiral with a fixed center radius $\rho = 0.02$(left), and with a reduced center radius $\rho = 4\Delta x$(right)
 [1] Feng Luo. A combinatorial curvature flow for compact 3-manifolds with boundary. Electronic Research Announcements, 2005, 11: 12-20. [2] Guido De Philippis, Antonio De Rosa, Jonas Hirsch. The area blow up set for bounded mean curvature submanifolds with respect to elliptic surface energy functionals. Discrete & Continuous Dynamical Systems - A, 2019, 39 (12) : 7031-7056. doi: 10.3934/dcds.2019243 [3] Simone Cacace, Maurizio Falcone. A dynamic domain decomposition for the eikonal-diffusion equation. Discrete & Continuous Dynamical Systems - S, 2016, 9 (1) : 109-123. doi: 10.3934/dcdss.2016.9.109 [4] Arseny Egorov. Morse coding for a Fuchsian group of finite covolume. Journal of Modern Dynamics, 2009, 3 (4) : 637-646. doi: 10.3934/jmd.2009.3.637 [5] Ronald E. Mickens. Positivity preserving discrete model for the coupled ODE's modeling glycolysis. Conference Publications, 2003, 2003 (Special) : 623-629. doi: 10.3934/proc.2003.2003.623 [6] Shu-Yu Hsu. Existence and properties of ancient solutions of the Yamabe flow. Discrete & Continuous Dynamical Systems - A, 2018, 38 (1) : 91-129. doi: 10.3934/dcds.2018005 [7] Matthias Erbar, Jan Maas. Gradient flow structures for discrete porous medium equations. Discrete & Continuous Dynamical Systems - A, 2014, 34 (4) : 1355-1374. doi: 10.3934/dcds.2014.34.1355 [8] Bin Pei, Yong Xu, Yuzhen Bai. Convergence of p-th mean in an averaging principle for stochastic partial differential equations driven by fractional Brownian motion. Discrete & Continuous Dynamical Systems - B, 2020, 25 (3) : 1141-1158. doi: 10.3934/dcdsb.2019213 [9] Caifang Wang, Tie Zhou. The order of convergence for Landweber Scheme with $\alpha,\beta$-rule. Inverse Problems & Imaging, 2012, 6 (1) : 133-146. doi: 10.3934/ipi.2012.6.133 [10] Hakan Özadam, Ferruh Özbudak. A note on negacyclic and cyclic codes of length $p^s$ over a finite field of characteristic $p$. Advances in Mathematics of Communications, 2009, 3 (3) : 265-271. doi: 10.3934/amc.2009.3.265 [11] Armin Lechleiter, Tobias Rienmüller. Factorization method for the inverse Stokes problem. Inverse Problems & Imaging, 2013, 7 (4) : 1271-1293. doi: 10.3934/ipi.2013.7.1271 [12] Wolf-Jüergen Beyn, Janosch Rieger. The implicit Euler scheme for one-sided Lipschitz differential inclusions. Discrete & Continuous Dynamical Systems - B, 2010, 14 (2) : 409-428. doi: 10.3934/dcdsb.2010.14.409 [13] Alina Chertock, Alexander Kurganov, Mária Lukáčová-Medvi${\rm{\check{d}}}$ová, Șeyma Nur Özcan. An asymptotic preserving scheme for kinetic chemotaxis models in two space dimensions. Kinetic & Related Models, 2019, 12 (1) : 195-216. doi: 10.3934/krm.2019009 [14] Peter Benner, Jens Saak, M. Monir Uddin. Balancing based model reduction for structured index-2 unstable descriptor systems with application to flow control. Numerical Algebra, Control & Optimization, 2016, 6 (1) : 1-20. doi: 10.3934/naco.2016.6.1 [15] M. Grasselli, V. Pata. Asymptotic behavior of a parabolic-hyperbolic system. Communications on Pure & Applied Analysis, 2004, 3 (4) : 849-881. doi: 10.3934/cpaa.2004.3.849 [16] Elena Bonetti, Pierluigi Colli, Gianni Gilardi. Singular limit of an integrodifferential system related to the entropy balance. Discrete & Continuous Dynamical Systems - B, 2014, 19 (7) : 1935-1953. doi: 10.3934/dcdsb.2014.19.1935 [17] Dmitry Treschev. A locally integrable multi-dimensional billiard system. Discrete & Continuous Dynamical Systems - A, 2017, 37 (10) : 5271-5284. doi: 10.3934/dcds.2017228 [18] Qiang Guo, Dong Liang. An adaptive wavelet method and its analysis for parabolic equations. Numerical Algebra, Control & Optimization, 2013, 3 (2) : 327-345. doi: 10.3934/naco.2013.3.327 [19] Tomáš Roubíček. An energy-conserving time-discretisation scheme for poroelastic media with phase-field fracture emitting waves and heat. Discrete & Continuous Dynamical Systems - S, 2017, 10 (4) : 867-893. doi: 10.3934/dcdss.2017044 [20] Chih-Chiang Fang. Bayesian decision making in determining optimal leased term and preventive maintenance scheme for leased facilities. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2020127

2019 Impact Factor: 1.233