Hostname: page-component-586b7cd67f-t8hqh Total loading time: 0 Render date: 2024-11-26T04:06:04.163Z Has data issue: false hasContentIssue false
  • English
  • Français

Interior Penalty Discontinuous Galerkin Based Isogeometric Analysis for Allen-Cahn Equations on Surfaces

Part of: Parabolic equations and systems Numerical approximation and computational geometry Partial differential equations, initial value and time-dependent initial-boundary value problems

Published online by Cambridge University Press:  23 November 2015

Futao Zhang ,
Yan Xu ,
Falai Chen  and
Ruihan Guo
Futao Zhang*
Affiliation:
School of Mathematical Sciences, University of Science and Technology of China, Hefei, Anhui 230026, P.R. China.
Yan Xu
Affiliation:
School of Mathematical Sciences, University of Science and Technology of China, Hefei, Anhui 230026, P.R. China.
Falai Chen
Affiliation:
School of Mathematical Sciences, University of Science and Technology of China, Hefei, Anhui 230026, P.R. China.
Ruihan Guo
Affiliation:
School of Mathematical Sciences, University of Science and Technology of China, Hefei, Anhui 230026, P.R. China.
*
*Corresponding author. Email addresses:[email protected] (F. Zhang), [email protected] (Y. Xu), chenf [email protected] (F. Chen), [email protected] (R. Guo)
Get access

Abstract

We propose a method that combines Isogeometric Analysis (IGA) with the interior penalty discontinuous Galerkin (IPDG) method for solving the Allen-Cahn equation, arising from phase transition in materials science, on three-dimensional (3D) surfaces consisting of multiple patches. DG ideology is adopted at patch level, i.e., we employ the standard IGA within each patch, and employ the IPDG method across the patch interfaces. IGA is very suitable for solving Partial Differential Equations (PDEs) on (3D) surfaces and the IPDG method is used to glue the multiple patches together to get the right solution. Our method takes advantage of both IGA and the IPDG method, which allows us to design a superior semi-discrete (in time) IPDG scheme. First and most importantly, the time-consuming mesh generation process in traditional Finite Element Analysis (FEA) is no longer necessary and refinements, including h-refinement and p-refinement which both maintain the original geometry, can be easily performed at any level. Moreover, the flexibility of the IPDG method makes our method very easy to handle cases with non-conforming patches and different degrees across the patch interfaces. Additionally, the geometrical error is eliminated (for all conic sections) or significantly reduced at the beginning due to the geometric flexibility of IGA basis functions, especially the use of multiple patches. Finally, this method can be easily formulated and implemented. We present our semi-discrete IPDG scheme after generally describe the problem, and then briefly introduce the time marching method employed in this paper. Theoretical analysis is carried out to show that our method satisfies a discrete energy law, and achieves the optimal convergence rate with respect to the L2 norm. Furthermore, we propose an elliptic projection operator on (3D) surfaces and prove an approximation error estimate which are vital for us to obtain the error estimate in the L2 norm. Numerical tests are given to validate the theory and gauge the good performance of our method.

MSC classification

Secondary: 65M60: Finite elements, Rayleigh-Ritz and Galerkin methods, finite methods 65D17: Computer aided design (modeling of curves and surfaces) 68U07: Computer-aided design 35K93: Quasilinear parabolic equations with mean curvature operator
Type
Research Article
Information
Communications in Computational Physics , Volume 18 , Issue 5 , November 2015 , pp. 1380 - 1416
Copyright
Copyright © Global-Science Press 2015 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Alexander, R., Diagonally implicit Runge-Kutta methods for stiff O.D.E'.s, SIAM Journal on Numerical Analysis, 14 (1977), pp. 1006–1021.Google Scholar
[2]Allen, S. M. and Cahn, J. W., A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening, Acta Metallurgica, 27 (1979), pp. 1084–1095.Google Scholar
[3]Arnold, D. N., Brezzi, F., Cockburn, B., and Marini, L. D., Unified analysis of Discontinuous Galerkin methods for elliptic problems, SIAM Journal on Numerical Analysis, 39 (2002), pp. 1749–1779.Google Scholar
[4]Auricchio, F., Beirão de Veiga, L., Buffa, A., and Sangalli, G., A fully “locking-free” isogeometric approach for plane linear elasticity problems: A stream function formulation, Computer Methods in Applied Mechanics and Engineering, 197.1 (2007), pp. 160–172.Google Scholar
[5]Bartels, S. and Müller, R., Quasi-optimal and robust a posteriori error estimates in L∞(L2)for the approximation of Allen-Cahn equations past singularities, Mathematics of Computation, 80 (2011), pp. 761–780.Google Scholar
[6]Bazilevs, Y., Beirão de Veiga, L., Cottrell, J. A., Hughes, T. J. R., and Sangalli, G., Isogeometric analysis: approximation, stability and error estimates for h-refined meshes, Mathematical Models and Methods in Applied Sciences, 16.07 (2006), pp. 1031–1090.Google Scholar
[7]Bazilevs, Y., Calo, V. M., Cottrell, J. A., Hughes, T. J. R., Reali, A., and Scovazzi, G., Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows, Computer Methods in Applied Mechanics and Engineering, 197 (2007), pp. 173–201.Google Scholar
[8]Bazilevs, Y. and Hughes, T. J. R., NURBS-based isogeometric analysis for the computation of flows about rotating components, Computational Mechanics, 43 (1) (2008), pp. 143–150.Google Scholar
[9]Bazilevs, Y., Calo, V. M., Hughes, T. J. R., and Zhang, Y., Isogeometric fluid-structure interaction: Theory, algorithms, and computations, Computational Mechanics, 43 (2008), pp. 3–37.Google Scholar
[10]Brunero, F., Pavarino, L. F., and Pechstein, C., Discontinuous Galerkin Methods for Isogeometric Analysis, Diploma thesis, 2012.Google Scholar
[11]Chen, L. Q., Phase-field models for microstructure evolution, Annual Review of Materials Research, 32 (2002), pp. 113–140.Google Scholar
[12]Chen, X., Generation and propagation of interfaces for reaction-diffusion equations, Journal of Differential equations, 96 (1) (1992), pp. 116–141.Google Scholar
[13]Chen, X., Elliott, C. M., Gardiner, A., and Zhao, J. J., Convergence of numerical solutions to the Allen-Cahn equation, Applicable Analysis, 69 (1998), pp. 47–56.Google Scholar
[14]Choi, J. W., Lee, H. G., Jeong, D., and Kim, J., An unconditionally gradient stable numerical method for solving the Allen-Cahn equation, Physica A: Statistical Mechanics and its Applications, 388(2009), pp. 1791–1803.Google Scholar
[15]Cockburn, B., Karniadakis, G., and Shu, C. W., The development of discontinuous Galerkin methods. In Discontinuous Galerkin Methods: Theory, Computation and Applications, Cock-burn, B., Karniadakis, G. and Shu, C. W., editors, Lecture Notes in Computational Science and Engineering, 11, Springer, (2000), Part I: Overview, pp. 3–50.Google Scholar
[16]Cottrell, J. A., Reali, A., Bazilevs, Y., and Hughes, T. J. R., Isogeometric analysis of structural vibrations, Computer Methods in Applied Mechanics and Engineering, 195 (2006), pp. 5257–5296.CrossRefGoogle Scholar
[17]Cottrell, J. A., Hughes, T. J. R., Reali, A., Studies of refinement and continuity in isogeometric analysis, Computer Methods in Applied Mechanics and Engineering, 196 (2007), pp. 4160–4183.Google Scholar
[18]Cottrell, J. A., Hughes, T. J. R., and Bazilevs, Y., Isogeometric Analysis: Toward Integration of CAD and FEA, Chichester: Wiley; 2009.Google Scholar
[19]de Mottoni, P. and Schatzman, M., Geometrical evolution of developed interfaces, Transactions of the American Mathematical Society, 347(5) (1995), pp. 1533–1589.Google Scholar
[20]Dziuk, G. and Elliot, C. M., Finite element methods for surface PDEs, Acta Numerica, 22 (2013), pp. 289–396.Google Scholar
[21]Evans, L.C, Soner, H.M., and Souganidis, P.E., Phase transitions and generalized motion by mean curvature, Communications on Pure and Applied Mathematics, 45(9) (1992), pp. 1097–1123.Google Scholar
[22]Evans, J.A. and Hughes, T.J.R., Isogeometric divergence-conforming B-splinesfor the steady Navier-Stokes equations, Mathematical Models and Methods in Applied Sciences, 23(08) (2013), pp. 1421–1478.Google Scholar
[23]Farin, G., Curves and Surfaces for CAGD, 5th ed., Morgan Kaufmann Publishers, 2002.Google Scholar
[24]Feng, X. B. and Prohl, A., Numerical analysis of the Allen-Cahn equation and approximation for mean curvature flows, Numerische Mathematik, 94 (2003), pp. 33–65.Google Scholar
[25]Feng, X. B. and Wu, H. J., A posteriori error estimates and an adaptive finite element method for the Allen-Cahn equation and the mean curvature flow, Journal of Scientific Computing, 24 (2005), pp. 121–146.Google Scholar
[26]Feng, X. B. and Li, Y. K., Analysis of interior penalty discontinuous Galerkin methods for the Allen-Cahn equation and the mean curvature flow, IMA Journal of Numerical Analysis, 2014.Google Scholar
[27]Fife, P. C., Dynamics of Internal Layers and Diffusive Interfaces, Society for Industrial and Applied Mathematics, Philadelphia, PA, 1988.Google Scholar
[28]Guo, R., Ji, L. and Xu, Y., High order local discontinuous Galerkin methods for the Allen-Cahn equation: analysis and simulation, Journal of Computational Mathematics, to appear.Google Scholar
[29]Hughes, T. J. R., Cottrell, J. A., and Bazilevs, Y., Isogeometric analysis: CAD, finite elements, NURBS, exact geometry, and mesh refinement, Computer Methods in Applied Mechanics and Engineering, 194 (2005), pp. 4135–4195.Google Scholar
[30]Kessler, D., Nochetto, R. H., and Schmidt, A., A posteriori error control for the Allen-Cahn problem: Circumventing Gronwall's inequality, Mathematical Modelling and Numerical Analysis, 38 (2004), pp. 129–142.Google Scholar
[31]McFadden, G.B., Phase field models of solidification. In Recent Advances in Numerical Methods for Partial Differential Equations and Applications, Contemporary Mathematics, 306 (2002), American Mathematical Society, Providence, RI, pp. 107–145.Google Scholar
[32]Nguyen, V.P., Kerfriden, P., Brino, M., et al., Nitsche's method for two and three dimensional NURBS patch coupling, Computational Mechanics, 53 (2014), pp. 1163–1182.Google Scholar
[33]Özcan, A.I., A Nitsche-based mortar method for non-matching multi-patch discretizations in the finite cell method, Master's thesis, Technische Universität München, Germany, 2012.Google Scholar
[34]Ruess, M., Schillinger, D., et al., Weak coupling for isogeometric analysis of non-matching and trimmed multi-patch geometries, Computer Methods in Applied Mechanics and Engineering, 269 (2014), pp. 46–71.Google Scholar
[35]Wheeler, A.A., Boettinger, W.J., and McFadden, G.B., Phase-field model for isothermal phase transitions in binary alloys, Physical Review A, 45 (1992), pp. 7424–7439.Google Scholar
[36]Xia, Y., Xu, Y., and Shu, C.W., Efficient time discretization for local discontinuous Galerkin methods, Discrete and Continuous Dynamical Systems Series B, 8 (2007), pp. 677–693.Google Scholar
[37]Xu, G.L. and Zhang, Q., Geometric Partial Differential Equation Methods in Computational Geometry, Beijing: Science Press; 2013.Google Scholar
[38]Yang, X., Error analysis of stabilized semi-implicit method of Allen-Cahn equation, Discrete and Continuous Dynamical Systems Series B, 11 (2009), pp. 1057–1070.Google Scholar
[39]Zhang, F., Xu, Y., and Chen, F., Discontinuous Galerkin methods for isogeometric analysis for elliptic equations on surfaces, Communications in Mathematics and Statistics, 2 (2014), pp. 431–461.CrossRefGoogle Scholar
[40]Zhang, J. and Du, Q., Numerical studies of discrete approximations to the Allen-Cahn equation in the sharp interface limit, SIAM Journal on Scientific Computing, 31 (2009), pp. 3042–3063.Google Scholar