Provably unconditionally stable, second-order time-accurate, mixed variational methods for phase-field models

Journal of Computational Physics

Volumn 230, Issue 13, 2011, Pages 5310-5327

  Gomez, Hector ^a   Hughes, Thomas J R ^b  

^a UNIVERSITY OF A CORUÑA   (Spain)

^b UNIVERSITY OF TEXAS AT AUSTIN   (United States)

Author keywords

Cahn Hilliard; Mixed finite element; Nonlinear stability; Phase field; Time integration

Indexed keywords

FINITE ELEMENT METHOD; INTEGRATION; NUMERICAL METHODS; PHASE INTERFACES;

CAHN-HILLIARD; MIXED FINITE ELEMENTS; MIXED VARIATIONAL METHODS; NON-LINEAR STABILITIES; PHASE FIELD MODELS; PHASE FIELDS; SECOND ORDERS; TIME-ACCURATE; TIME-INTEGRATION; UNCONDITIONALLY STABLE;

FREE ENERGY;

EID: 79955469960     PISSN: 00219991     EISSN: 10902716     Source Type: Journal    
DOI: 10.1016/j.jcp.2011.03.033     Document Type: Article

References (76)

1. Akkerman I., Bazilevs Y., Calo V.M., Hughes T.J.R., Hulshoff S. The role of continuity in residual-based variational multiscale modeling of turbulence. Computational Mechanics 2007, 41:371-378.
2. Anderson D.M., McFadden G.B., Wheeler A.A. Diffuse-interface methods in fluid mechanics. Annual Review of Fluid Mechanics 1998, 30:139-165.
3. Armero F., Zambrana-Rojas C. Volume-preserving energy-momentum schemes for isochoric multiplicative plasticity. Computer Methods in Applied Mechanics and Engineering 2007, 196:4130-4159.
4. Auricchio F., Beirao da Veiga L., Hughes T.J.R., Reali A., Sangalli G. Isogeometric collocation methods. Mathematical Models and Methods in Applied Sciences 2010, 2:2075-2107.
5. Barenblatt G.I. Scaling, Self-similarity, and Intermediate Asymptotics 1996, Cambridge University Press.
6. Bazilevs Y., Calo V.M., Cottrell J.A., Evans J.A., Hughes T.J.R., Lipton S., Scott M.A., Sederberg T.W. Isogeometric Analysis using T-splines. Computer Methods in Applied Mechanics and Engineering 2010, 199:229-263.
7. Bazilevs Y., Calo V.M., Cottrell J.A., Hughes T.J.R., Reali A., Scovazzi G. Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows. Computer Methods in Applied Mechanics and Engineering 2007, 197:173-201.
8. Bazilevs Y., Hughes T.J.R. NURBS-based isogeometric analysis for the computation of flows about rotating components. Computational Mechanics 2008, 43:143-150.
9. Becker J., Grün G., Seemann R., Mantz H., Jacobs K., Mecke K.R., Blossey R. Complex dewetting scenarios captured by thin-film models. Nature Materials 2003, 2:59-63.
10. Benderskaya G., Clemens M., De Gersem H., Weiland T. Embedded Runge-Kutta methods for field-circuit coupled problems with switching elements. IEEE Transactions on Magnetics 2005, 41:1612-1615.
11. Buffa A., Sangalli G., Vázquez R. Isogeometric analysis in electromagnetics: B-splines approximation. Computer Methods in Applied Mechanics and Engineering 2010, 199:1143-1152.
12. Cahn J.W., Hilliard J.E. Free energy of a non-uniform system. I. Interfacial free energy. The Journal of Chemical Physics 1958, 28:258-267.
13. Cahn J.W., Hilliard J.E. Free energy of a non-uniform system. III. Nucleation in a two-component incompressible fluid. The Journal of Chemical Physics 1959, 31:688-699.
14. Chen L.Q. Phase-field models for microstructural evolution. Annual Review Materials Research 2002, 32:113-140.
15. Choksi R., Peletier M.A., Williams J.F. On the phase diagram for microphase separation of diblock copolymers: an approach via a nonlocal Cahn-Hilliard functional. SIAM Journal of Applied Mathematics 2009, 69:1712-1738.
16. Chung J., Hulbert G.M. A time integration algorithm for structural dynamics with improved numerical dissipation: the generalized-α method. Journal of Applied Mechanics 1993, 60:371-375.
17. Copetti M.I.M., Elliott C.M. Numerical analysis of the Cahn-Hilliard equation with a logarithmic free energy. Numerische Mathematik 1992, 63:39-65.
18. Cottrell J.A., Hughes T.J.R., Bazilevs Y. Isogeometric Analysis: Toward integration of CAD and FEA 2009, Wiley.
19. Cottrell J.A., Hughes T.J.R., Reali A. Studies of refinement and continuity in isogeometric structural analysis. Computer Methods in Applied Mechanics and Engineering 2007, 196:4160-4183.
20. Cristini V., Li X., Lowengrub J.S., Wise S.M. Nonlinear simulations of solid tumor growth using a mixture model: invasion and branching. Journal of Mathematical Biology 2009, 58:723-763.
21. Cueto-Felgueroso L., Juanes R. Nonlocal interface dynamics and pattern formation in gravity-driven unsaturated flow through porous media. Physical Review Letters 2008, 101:244504.
22. Cueto-Felgueroso L., Juanes R. A phase-field model of unsaturated flow. Water Resources Research 2009, 45:W10409.
23. Cueto-Felgueroso L., Peraire J. A time-adaptive finite volume method for the Cahn-Hilliard and Kuramoto-Sivashinsky equations. Journal of Computational Physics 2008, 227:9985-10017.
24. Du Q., Nicolaides R.A. Numerical analysis of a continuum model of phase transition. SIAM Journal of Numerical Analysis 1991, 28:1310-1322.
25. Elguedj T., Bazilevs Y., Calo V.M., Hughes T.J.R. B̄ and F̄ projection methods for nearly incompressible linear and non-linear elasticity and plasticity using higher-order NURBS elements. Computer Methods in Applied Mechanics and Engineering 2008, 197:2732-2762.
26. Elliott C.M., Garcke H. On the Cahn-Hilliard equation with degenerate mobility. SIAM Journal of Mathematical Analysis 1996, 27:404-423.
27. Emmerich H. The Diffuse Interface Approach in Materials Science 2003, Springer.
28. D.J. Eyre, An unconditionally stable one-step scheme for gradient systems, unpublished. http://www.math.utah.edu/~eyre/research/methods/stable.ps.
29. Evans J.A., Bazilevs Y., Babuška I., Hughes T.J.R. n-widths, sup infs, and optimality ratios for the k-version of the isogeometric finite element method. Computer Methods in Applied Mechanics and Engineering 2009, 198:1726-1741.
30. Feng X., Wu H. A posteriori error estimates for finite element approximations of the Cahn-Hilliard equation and the Hele-Shaw flow. Journal of Computational Mathematics 2008, 26:767-796.
31. Fonseca I., Morini M. Surfactants in foam stability: a phase field model. Archive for Rational Mechanics 2007, 183:411-456.
32. Frieboes H.B., Lowengrub J.S., Wise S., Zheng X., Macklin P., Bearer E.L., Cristini V. Computer simulation of glioma growth and morphology. NeuroImage 2007, 37:59-70.
33. Fried E., Gurtin M.E. Dynamic solid-solid transitions with phase characterized by an order parameter. Phys. D 1994, 72:287-308.
34. Furihata D. A stable and conservative finite difference scheme for the Cahn-Hilliard equation. Numerische Mathematik 2001, 87:675-699.
35. Gomez H., Calo V.M., Bazilevs Y., Hughes T.J.R. Isogeometric analysis of the Cahn-Hilliard phase-field model. Computer Methods in Applied Mechanics and Engineering 2008, 197:4333-4352.
36. Gomez H., Hughes T.J.R., Nogueira X., Calo V.M. Isogeometric analysis of the isothermal Navier-Stokes-Korteweg equations. Computer Methods in Applied Mechanics and Engineering 2010, 199:1828-1840.
37. Gurtin M.E. Generalized Ginzburg-Landau and Cahn-Hilliard equations based on a microforce balance. Physica D 1996, 92:178-192.
38. Gustafsson K. Control-theoretic techniques for stepsize selection in implicit Runge-Kutta methods. ACM Transactions on Mathematical Software 1994, 20:496-517.
39. Harten A. On the symmetric form of systems of conservation laws with entropy. Journal of Computational Physics 1983, 49:151-164.
40. He L., Liu Y. A class of stable spectral methods for the Cahn-Hilliard equation. Journal of Computational Physics 2009, 228:5101-5110.
41. Hu Z., Wise S.M., Wang C., Lowengrub J.S. Stable and efficient finite-difference nonlinear multigrid schemes for the phase field crystal equation. Journal of Computational Physics 2009, 228:5323-5339.
42. Huber J.E., Fleck N.A., Landis C.M., McMeeking R.M. A constitutive model for ferroelectric polycrystals. Journal of the Mechanics and Physics of Solids 2009, 47:1663-1697.
43. Hughes T.J.R. The Finite Element Method: Linear Static and Dynamic Finite Element Analysis 2000, Dover Publications, Mineola, NY.
44. Hughes T.J.R., Cottrell J.A., Bazilevs Y. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Computer Methods in Applied Mechanics and Engineering 2005, 194:4135-4195.
45. Hughes T.J.R., Franca L.P., Mallet M. A new finite element formulation for computational fluid dynamics: I. Symmetric forms of the compressible Euler and Navier-Stokes equations and the second law of thermodynamics. Computer Methods in Applied Mechanics and Engineering 1986, 54:223-234.
46. R Hughes T.J., Reali A., Sangalli G. Duality and unified analysis of discrete approximations in structural dynamics and wave propagation: comparison of p-method finite elements with k-method NURBS. Computer Methods in Applied Mechanics and Engineering 2008, 197:4104-4124.
47. Jansen K.E., Whiting C.H., Hulbert G.M. A generalized-α method for integrating the filtered Navier-Stokes equations with a stabilized finite element method. Computer Methods in Applied Mechanics and Engineering 1999, 190:305-319.
48. Jeong J.-H., Goldenfeld N., Dantzig J.A. Phase field model for three-dimensional dendritic growth with fluid flow. Physical Review E 2001, 64:1-14.
49. Kay D., Welford R. A multigrid finite element solver for the Cahn-Hilliard equation. Journal of Computational Physics 2006, 212:288-304.
50. Kim J. A numerical method for the Cahn-Hilliard equation with a variable mobility. Communications in Nonlinear Science and Numerical Simulation 2007, 12:1560-1571.
51. Kim Y-T., Provatas N., Goldenfeld N., Dantzig J.A. universal dynamics of phase field models for dendritic growth. Physical Review E 1999, 59:2546-2549.
52. Kobayashi R. A numerical approach to three-dimensional dendritic solidification. Experimental Mathematics 1994, 3:59-81.
53. L.D. Landau, V.I. Ginzburg, On the theory of superconductivity, in: Collected papers of L.D. Landau, D. ter Haar, ed., Pergamon Oxford (1965) 626-633.
54. Landis C.M. Fully coupled multi-axial, symmetric constitutive laws for polycrystalline ferroelectric ceramics. Journal of the Mechanics and Physics of Solids 2002, 50:127-152.
55. Lang J. Two-dimensional fully adaptive solutions of reaction-diffusion equations. Applied Numerical Mathematics 1995, 18:223-240.
56. Lipton S., Evans J.A., Bazilevs Y., Elguedj T., Hughes T.J.R. Robustness of isogeometric structural discretizations under severe mesh distortion. Computer Methods in Applied Mechanics and Engineering 2010, 199:357-373.
57. Mahapatra D.R., Melnik R.V.N. Finite element analysis of phase transformation dynamics in shape memory alloys with a consistent Landau-Ginzburg free energy model. Mechanics of Advanced Materials and Structures 2006, 13:443-455.
58. Mahapatra D.R., Melnik R.V.N. Finite element approach to modelling evolution of 3D shape memory materials. Mathematics and Computers in Simulation 2007, 76:141-148.
59. de Mello E.V.L., Teixeira da Silveira Filho O. Numerical study of the Cahn-Hilliard equation in one, two and three dimensions. Physica A 2005, 347:429-443.
60. Meng X.N., Laursen T.A. Energy consistent algorithms for dynamic finite deformation plasticity. Computer Methods in Applied Mechanics and Engineering 2002, 191:1639-1675.
61. Oden J.T., Hawkins A., Prudhomme S. General diffuse-interface theories and an approach to predictive tumor growth modeling. Mathematical Models and Methods in Applied Sciences 2010, 20:477-517.
62. Penrose O., Fife P.C. Thermodynamically consistent models of phase field type for the kinetics of phase transition. Physica D 1990, 43:44-62.
63. Piegl L., Tiller W. The NURBS Book, Monographs in Visual Communication 1997, Springer-Verlag.
64. Rajagopal A., Fischer P., Kuhl E., Steinmann P. Natural element analysis of the Cahn-Hilliard phase-field model. Computational Mechanics 2010, 46:471-493.
65. Rogers D.F. An Introduction to NURBS With Historical Perspective 2001, Academic Press.
66. Romero I. Algorithms for coupled problems that preserve symmetries and the laws of thermodynamics: Part I: Monolithic integrators and their application to finite strain thermoelasticity. Computer Methods in Applied Mechanics and Engineering 2010, 199:1841-1858.
67. Romero I. Algorithms for coupled problems that preserve symmetries and the laws of thermodynamics: Part II: Fractional step methods. Computer Methods in Applied Mechanics and Engineering 2010, 199:2235-2248.
68. Saad Y., Schultz M.H. GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM Journal of Scientific and Statistical Computing 1986, 7:856-869.
69. Shakib F., Hughes T.J.R., Johan Z. A new finite element formulation for computational fluid-dynamics. 10. The compressible Euler and Navier-Stokes equations. Computer Methods in Applied Mechanics and Engineering 1991, 89:141-219.
70. Tadmor E. Skew-self adjoint form for systems of conservation laws. Journal of Mathematical Analysis and Applications 1984, 103:428-442.
71. Tremaine S. On the origin of irregular structure in Saturn's rings. Astronomical Journal 2003, 125:894-901.
72. van der Houwen P.J., Sommeijer B.P., Couzy W. Embedded diagonally implicit Runge-Kutta algorithms on parallel computers. Mathematics of Computation 1992, 58:135-159.
73. Vollmayr-Lee B.P., Rutemberg A.D. Fast and accurate coarsening simulation with an unconditionally stable time step. Physical Review E 2003, 68:066703.
74. Wells G.N., Kuhl E., Garikipati K. A discontinuous Galerkin method for the Cahn-Hilliard equation. Journal of Computational Physics 2006, 218:860-877.
75. Xu C., Tang T. Stability analysis of large time-stepping methods for epitaxial growth models. SIAM Journal of Numerical Analysis 2006, 44:1759-1779.
76. van der Waals J.D. The thermodynamic theory of capillarity under the hypothesis of a continuous variation of density. Journal of Statistical Physics 1979, 20:197-244.