Figure 1.
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
-
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 |
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.
Keywords:
Crystalline eikonal-curvature flow,
evolution of a polygonal spiral,
ODE system for crystalline motion,
level set method,
finite difference scheme.
Mathematics Subject Classification:
Primary: 34A34, 53E10; Secondary: 53A04.
Citation:
Tetsuya Ishiwata, Takeshi Ohtsuka. Numerical analysis of an ODE and a level set methods for evolving spirals by crystalline eikonal-curvature flow. Discrete and Continuous Dynamical Systems - S,
2021, 14
(3)
: 893-907.
doi: 10.3934/dcdss.2020390
![]()
References:
| [1] |
F. Almgren and J. E. Taylor,
Flat flow is motion by crystalline curvature for curves with crystalline energies, J. Differential Geom., 42 (1995), 1-22.
doi: 10.4310/jdg/1214457030. |
| [2] |
F. Almgren, J. E. Taylor and L. Wang,
Curvature-driven flows: A variational approach, SIAM J. Control Optim., 31 (1993), 387-438.
doi: 10.1137/0331020. |
| [3] |
S. Angenent and M. E. Gurtin,
Multiphase thermomechanics with interfacial structure. Ⅱ. Evolution of an isothermal interface, Arch. Rational Mech. Anal., 108 (1989), 323-391.
doi: 10.1007/BF01041068. |
| [4] |
G. Bellettini and M. Paolini,
Anisotropic motion by mean curvature in the context of Finsler geometry, Hokkaido Math. J., 25 (1996), 537-566.
doi: 10.14492/hokmj/1351516749. |
| [5] |
A. Chambolle,
An algorithm for mean curvature motion, Interfaces Free Bound., 6 (2004), 195-218.
doi: 10.4171/IFB/97. |
| [6] |
A. Chambolle, M. Morini, M. Novaga and M. Ponsiglione,
Existence and uniqueness for anisotropic and crystalline mean curvature flows, J. Amer. Math. Soc., 32 (2019), 779-824.
doi: 10.1090/jams/919. |
| [7] |
A. Chambolle, M. Morini and M. Ponsiglione,
Existence and uniqueness for a crystalline mean curvature flow, Comm. Pure Appl. Math., 70 (2017), 1084-1114.
doi: 10.1002/cpa.21668. |
| [8] |
B. Engquist, A.-K. Tornberg and R. Tsai,
Discretization of Dirac delta functions in level set methods, J. Comput. Phys., 207 (2005), 28-51.
doi: 10.1016/j.jcp.2004.09.018. |
| [9] |
M.-H. Giga and Y. Giga,
Generalized motion by nonlocal curvature in the plane, Arch. Ration. Mech. Anal., 159 (2001), 295-333.
doi: 10.1007/s002050100154. |
| [10] |
Y. Giga, Surface Evolution Equations: A Level Set Approach, Monographs in Mathematics, 99. Birkhäuser Verlag, Basel, 2006. |
| [11] |
Y. Giga and N. Požár, A level set crystalline mean curvature flow of surfaces, Adv. Differential Equations, 21 (2016), 631–698, http://projecteuclid.org/euclid.ade/1462298654. |
| [12] |
S. Goto, M. Nakagawa and T. Ohtsuka,
Uniqueness and existence of generalized motion for spiral crystal growth, Indiana University Mathematics Journal, 57 (2008), 2571-2599.
doi: 10.1512/iumj.2008.57.3350. |
| [13] |
M. E. Gurtin, Thermomechanics of Evolving Phase Boundaries in the Plane, Oxford Mathematical Monographs, The Clarendon Press, Oxford University Press, New York, 1993. |
| [14] |
T. Ishiwata,
Crystalline motion of spiral-shaped polygonal curves with a tip motion, Discrete Contin. Dyn. Syst. Ser. S, 7 (2014), 53-62.
doi: 10.3934/dcdss.2014.7.53. |
| [15] |
T. Ishiwata and T. Ohtsuka,
Evolution of spiral-shaped polygonal curve by crystalline curvature flow with a pinned tip, Discrete Contin. Dyn. Syst. Ser. B, 24 (2019), 5261-5295.
doi: 10.3934/dcdsb.2019058. |
| [16] |
A. Oberman, S. Osher, R. Takei and R. Tsai,
Numerical methods for anisotropic mean curvature flow based on a discrete time variational formulation, Commun. Math. Sci., 9 (2011), 637-662.
doi: 10.4310/CMS.2011.v9.n3.a1. |
| [17] |
T. Ohtsuka, Y.-H. R. Tsai and Y. Giga,
A level set approach reflecting sheet structure with single auxiliary function for evolving spirals on crystal surfaces, Journal of Scientific Computing, 62 (2015), 831-874.
doi: 10.1007/s10915-014-9877-2. |
| [18] |
T. Ohtsuka,
A level set method for spiral crystal growth, Advances in Mathematical Sciences and Applications, 13 (2003), 225-248.
|
| [19] |
T. Ohtsuka,
Minimizing movement approach for spirals evolving by crystalline curvature using level set functions, Oberwolfach Reports, 14 (2017), 314-317.
|
| [20] |
T. Ohtsuka, Minimizing movement approach without using distance function for evolving spirals by the crystalline curvature with driving force, RIMS Kôkyûroku No.2121, 74–87, http://www.kurims.kyoto-u.ac.jp/ kyodo/kokyuroku/contents/pdf/2121-06.pdf. |
| [21] |
T. Ohtsuka, Y.-H. R. Tsai and Y. Giga,
Growth rate of crystal surfaces with several dislocation centers, Crystal Growth & Design, 18 (2018), 1917-1929.
doi: 10.1021/acs.cgd.7b00833. |
| [22] |
R. T. Rockafellar, Convex Analysis, Princeton Mathematical Series, No. 28, Princeton University Press, Princeton, N.J., 1970. |
| [23] |
J. E. Taylor,
Constructions and conjectures in crystalline nondifferential geometry, Differential Geometry, Pitman Monogr. Surveys Pure Appl. Math., Longman Sci. Tech., Harlow, 52 (1991), 321-336.
|
show all references
References:
| [1] |
F. Almgren and J. E. Taylor,
Flat flow is motion by crystalline curvature for curves with crystalline energies, J. Differential Geom., 42 (1995), 1-22.
doi: 10.4310/jdg/1214457030. |
| [2] |
F. Almgren, J. E. Taylor and L. Wang,
Curvature-driven flows: A variational approach, SIAM J. Control Optim., 31 (1993), 387-438.
doi: 10.1137/0331020. |
| [3] |
S. Angenent and M. E. Gurtin,
Multiphase thermomechanics with interfacial structure. Ⅱ. Evolution of an isothermal interface, Arch. Rational Mech. Anal., 108 (1989), 323-391.
doi: 10.1007/BF01041068. |
| [4] |
G. Bellettini and M. Paolini,
Anisotropic motion by mean curvature in the context of Finsler geometry, Hokkaido Math. J., 25 (1996), 537-566.
doi: 10.14492/hokmj/1351516749. |
| [5] |
A. Chambolle,
An algorithm for mean curvature motion, Interfaces Free Bound., 6 (2004), 195-218.
doi: 10.4171/IFB/97. |
| [6] |
A. Chambolle, M. Morini, M. Novaga and M. Ponsiglione,
Existence and uniqueness for anisotropic and crystalline mean curvature flows, J. Amer. Math. Soc., 32 (2019), 779-824.
doi: 10.1090/jams/919. |
| [7] |
A. Chambolle, M. Morini and M. Ponsiglione,
Existence and uniqueness for a crystalline mean curvature flow, Comm. Pure Appl. Math., 70 (2017), 1084-1114.
doi: 10.1002/cpa.21668. |
| [8] |
B. Engquist, A.-K. Tornberg and R. Tsai,
Discretization of Dirac delta functions in level set methods, J. Comput. Phys., 207 (2005), 28-51.
doi: 10.1016/j.jcp.2004.09.018. |
| [9] |
M.-H. Giga and Y. Giga,
Generalized motion by nonlocal curvature in the plane, Arch. Ration. Mech. Anal., 159 (2001), 295-333.
doi: 10.1007/s002050100154. |
| [10] |
Y. Giga, Surface Evolution Equations: A Level Set Approach, Monographs in Mathematics, 99. Birkhäuser Verlag, Basel, 2006. |
| [11] |
Y. Giga and N. Požár, A level set crystalline mean curvature flow of surfaces, Adv. Differential Equations, 21 (2016), 631–698, http://projecteuclid.org/euclid.ade/1462298654. |
| [12] |
S. Goto, M. Nakagawa and T. Ohtsuka,
Uniqueness and existence of generalized motion for spiral crystal growth, Indiana University Mathematics Journal, 57 (2008), 2571-2599.
doi: 10.1512/iumj.2008.57.3350. |
| [13] |
M. E. Gurtin, Thermomechanics of Evolving Phase Boundaries in the Plane, Oxford Mathematical Monographs, The Clarendon Press, Oxford University Press, New York, 1993. |
| [14] |
T. Ishiwata,
Crystalline motion of spiral-shaped polygonal curves with a tip motion, Discrete Contin. Dyn. Syst. Ser. S, 7 (2014), 53-62.
doi: 10.3934/dcdss.2014.7.53. |
| [15] |
T. Ishiwata and T. Ohtsuka,
Evolution of spiral-shaped polygonal curve by crystalline curvature flow with a pinned tip, Discrete Contin. Dyn. Syst. Ser. B, 24 (2019), 5261-5295.
doi: 10.3934/dcdsb.2019058. |
| [16] |
A. Oberman, S. Osher, R. Takei and R. Tsai,
Numerical methods for anisotropic mean curvature flow based on a discrete time variational formulation, Commun. Math. Sci., 9 (2011), 637-662.
doi: 10.4310/CMS.2011.v9.n3.a1. |
| [17] |
T. Ohtsuka, Y.-H. R. Tsai and Y. Giga,
A level set approach reflecting sheet structure with single auxiliary function for evolving spirals on crystal surfaces, Journal of Scientific Computing, 62 (2015), 831-874.
doi: 10.1007/s10915-014-9877-2. |
| [18] |
T. Ohtsuka,
A level set method for spiral crystal growth, Advances in Mathematical Sciences and Applications, 13 (2003), 225-248.
|
| [19] |
T. Ohtsuka,
Minimizing movement approach for spirals evolving by crystalline curvature using level set functions, Oberwolfach Reports, 14 (2017), 314-317.
|
| [20] |
T. Ohtsuka, Minimizing movement approach without using distance function for evolving spirals by the crystalline curvature with driving force, RIMS Kôkyûroku No.2121, 74–87, http://www.kurims.kyoto-u.ac.jp/ kyodo/kokyuroku/contents/pdf/2121-06.pdf. |
| [21] |
T. Ohtsuka, Y.-H. R. Tsai and Y. Giga,
Growth rate of crystal surfaces with several dislocation centers, Crystal Growth & Design, 18 (2018), 1917-1929.
doi: 10.1021/acs.cgd.7b00833. |
| [22] |
R. T. Rockafellar, Convex Analysis, Princeton Mathematical Series, No. 28, Princeton University Press, Princeton, N.J., 1970. |
| [23] |
J. E. Taylor,
Constructions and conjectures in crystalline nondifferential geometry, Differential Geometry, Pitman Monogr. Surveys Pure Appl. Math., Longman Sci. Tech., Harlow, 52 (1991), 321-336.
|
Figure 2.
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
Figure 3.
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
Figure 4.
Profiles of the square spiral at $ t = 1 $ . The level set method is calculated using $ \rho = 0.02 $ and $ \Delta x = 0.0050 $
Figure 5.
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)
Figure 6.
Profiles of the diagonal spiral at $ t = 1 $ . The level set method is calculated using $ \rho = 0.02 $ and $ \Delta x = 0.0050 $
Figure 7.
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)
Figure 8.
Profiles of the triangle spiral at $ t = 0.8 $ . The level set method is calculated using $ \rho = 0.02 $ and $ \Delta x = 0.0050 $
Figure 9.
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] |
Tetsuya Ishiwata, Takeshi Ohtsuka. Evolution of a spiral-shaped polygonal curve by the crystalline curvature flow with a pinned tip. Discrete and Continuous Dynamical Systems - B, 2019, 24 (10) : 5261-5295. doi: 10.3934/dcdsb.2019058 |
| [2] |
Tetsuya Ishiwata. On the motion of polygonal curves with asymptotic lines by crystalline curvature flow with bulk effect. Discrete and Continuous Dynamical Systems - S, 2011, 4 (4) : 865-873. doi: 10.3934/dcdss.2011.4.865 |
| [3] |
Tetsuya Ishiwata. Crystalline motion of spiral-shaped polygonal curves with a tip motion. Discrete and Continuous Dynamical Systems - S, 2014, 7 (1) : 53-62. doi: 10.3934/dcdss.2014.7.53 |
| [4] |
Tetsuya Ishiwata. On spiral solutions to generalized crystalline motion with a rotating tip motion. Discrete and Continuous Dynamical Systems - S, 2015, 8 (5) : 881-888. doi: 10.3934/dcdss.2015.8.881 |
| [5] |
Tetsuya Ishiwata. Motion of polygonal curved fronts by crystalline motion: v-shaped solutions and eventual monotonicity. Conference Publications, 2011, 2011 (Special) : 717-726. doi: 10.3934/proc.2011.2011.717 |
| [6] |
Peter Frolkovič, Viera Kleinová. A new numerical method for level set motion in normal direction used in optical flow estimation. Discrete and Continuous Dynamical Systems - S, 2021, 14 (3) : 851-863. doi: 10.3934/dcdss.2020347 |
| [7] |
Mi-Ho Giga, Yoshikazu Giga, Ryo Kuroda, Yusuke Ochiai. Crystalline flow starting from a general polygon. Discrete and Continuous Dynamical Systems, 2022, 42 (4) : 2027-2051. doi: 10.3934/dcds.2021182 |
| [8] |
Mi-Ho Giga, Yoshikazu Giga. A subdifferential interpretation of crystalline motion under nonuniform driving force. Conference Publications, 1998, 1998 (Special) : 276-287. doi: 10.3934/proc.1998.1998.276 |
| [9] |
Yoshikazu Giga, Hiroyoshi Mitake, Hung V. Tran. Remarks on large time behavior of level-set mean curvature flow equations with driving and source terms. Discrete and Continuous Dynamical Systems - B, 2020, 25 (10) : 3983-3999. doi: 10.3934/dcdsb.2019228 |
| [10] |
Ronald E. Mickens. A nonstandard finite difference scheme for the drift-diffusion system. Conference Publications, 2009, 2009 (Special) : 558-563. doi: 10.3934/proc.2009.2009.558 |
| [11] |
Tetsuya Ishiwata, Shigetoshi Yazaki. A fast blow-up solution and degenerate pinching arising in an anisotropic crystalline motion. Discrete and Continuous Dynamical Systems, 2014, 34 (5) : 2069-2090. doi: 10.3934/dcds.2014.34.2069 |
| [12] |
Bendong Lou. Spiral rotating waves of a geodesic curvature flow on the unit sphere. Discrete and Continuous Dynamical Systems - B, 2012, 17 (3) : 933-942. doi: 10.3934/dcdsb.2012.17.933 |
| [13] |
Zhenlin Guo, Ping Lin, Guangrong Ji, Yangfan Wang. Retinal vessel segmentation using a finite element based binary level set method. Inverse Problems and Imaging, 2014, 8 (2) : 459-473. doi: 10.3934/ipi.2014.8.459 |
| [14] |
Yones Esmaeelzade Aghdam, Hamid Safdari, Yaqub Azari, Hossein Jafari, Dumitru Baleanu. Numerical investigation of space fractional order diffusion equation by the Chebyshev collocation method of the fourth kind and compact finite difference scheme. Discrete and Continuous Dynamical Systems - S, 2021, 14 (7) : 2025-2039. doi: 10.3934/dcdss.2020402 |
| [15] |
Maria Laura Delle Monache, Paola Goatin. A front tracking method for a strongly coupled PDE-ODE system with moving density constraints in traffic flow. Discrete and Continuous Dynamical Systems - S, 2014, 7 (3) : 435-447. doi: 10.3934/dcdss.2014.7.435 |
| [16] |
Jamie Cruz, Miguel Gutiérrez. Spiral motion in classical mechanics. Conference Publications, 2009, 2009 (Special) : 191-197. doi: 10.3934/proc.2009.2009.191 |
| [17] |
Francis Michael Russell, J. C. Eilbeck. Persistent mobile lattice excitations in a crystalline insulator. Discrete and Continuous Dynamical Systems - S, 2011, 4 (5) : 1267-1285. doi: 10.3934/dcdss.2011.4.1267 |
| [18] |
Annalisa Malusa, Matteo Novaga. Crystalline evolutions in chessboard-like microstructures. Networks and Heterogeneous Media, 2018, 13 (3) : 493-513. doi: 10.3934/nhm.2018022 |
| [19] |
Weizhu Bao, Chunmei Su. Uniform error estimates of a finite difference method for the Klein-Gordon-Schrödinger system in the nonrelativistic and massless limit regimes. Kinetic and Related Models, 2018, 11 (4) : 1037-1062. doi: 10.3934/krm.2018040 |
| [20] |
Jiangfeng Huang, Zhiliang Deng, Liwei Xu. A Bayesian level set method for an inverse medium scattering problem in acoustics. Inverse Problems and Imaging, 2021, 15 (5) : 1077-1097. doi: 10.3934/ipi.2021029 |
2020 Impact Factor: 2.425
Tools
Metrics
Other articles
by authors
[Back to Top]