Adaptive phase-field modeling of three-dimensional thermally-induced brittle fracture in COMSOL
Source abstract
Abstract
Abstract This work proposes an adaptive thermomechanical phase-field framework for simulating three-dimensional (3D) thermally-induced fracture in brittle materials. Adaptive mesh refinement is performed via the COMSOL application programming interface (API), using a user-defined phase-field threshold as the refinement indicator. COMSOL’s built-in adaptive meshing produces a locally refined mesh without hanging nodes. A staggered iteration scheme is adopted to ensure computational robustness of the fracture phase-field framework. To further improve computational efficiency in 3D fracture problems, the restarted generalized minimum residual (GMRES) method, combined with smoothed aggregation algebraic multigrid (SAAMG), is adopted as an iterative linear system solver. Through the analysis of three numerical examples, the accuracy and efficiency of the present framework can be validated. The open-source model enables users without programming experience to simulate fractures in a multiphysics environment.
Source: openaire · Quality: candidate unverified · Source record ↗
Research context
OpenAlex topics: Numerical methods in engineering · Advanced Numerical Methods in Computational Mathematics · Solidification and crystal growth phenomena
SciVX discovery tags are terminology-based. OpenAlex topics are source-provided classifications.
References (81)
Show source-provided references
1. Agathos K, Bordas SP, Chatzi E (2019) Improving the conditioning of XFEM/GFEM for fracture mechanics problems through enrichment quasi-orthogonalization. Comput Methods Appl Mech Eng 346:1051–107310.1016/j.cma.2018.08.007 ↗
2. Agathos K, Chatzi E, Bordas SP (2016) Stable 3D extended finite elements with higher order enrichment for accurate non planar fracture. Comput Methods Appl Mech Eng 306:19–4610.1016/j.cma.2016.03.023 ↗
3. Agathos K, Chatzi E, Bordas SP, Talaslidis D (2016) A well-conditioned and optimally convergent XFEM for 3D linear elastic fracture. Int J Numer Meth Eng 105:643–67710.1002/nme.4982 ↗
4. Ambati M, Gerasimov T, De Lorenzis L (2015) A review on phase-field models of brittle fracture and a new fast hybrid formulation. Comput Mech 55:383–40510.1007/s00466-014-1109-y ↗
5. Amor H, Marigo J-J, Maurini C (2009) Regularized formulation of the variational brittle fracture with unilateral contact: Numerical experiments. J Mech Phys Solids 57:1209–122910.1016/j.jmps.2009.04.011 ↗
6. Badnava H, Msekh MA, Etemadi E, Rabczuk T (2018) An h-adaptive thermo-mechanical phase field model for fracture. Finite Elem Anal Des 138:31–4710.1016/j.finel.2017.09.003 ↗
7. Becker R, Bordas SP, Chouly F, Omnes P (2024) A short perspective on a posteriori error control and adaptive discretizations. Adv Appl Mech 58:1–36
8. Bordas SP, Menk A, Natarajan S (2023) Partition of unity methods, John Wiley & Sons10.1002/9781118535875 ↗
9. Bordas S, Duflot M (2007) Derivative recovery and a posteriori error estimate for extended finite elements. Comput Methods Appl Mech Eng 196:3381–339910.1016/j.cma.2007.03.011 ↗
10. Bordas S, Duflot M, Le P (2008) A simple error estimator for extended finite elements. Commun Numer Methods Eng 24:961–97110.1002/cnm.1001 ↗
11. Borden MJ, Hughes TJ, Landis CM, Verhoosel CV (2014) A higher-order phase-field model for brittle fracture: Formulation and analysis within the isogeometric analysis framework. Comput Methods Appl Mech Eng 273:100–11810.1016/j.cma.2014.01.016 ↗
12. Borden MJ, Verhoosel CV, Scott MA, Hughes TJ, Landis CM (2012) A phase-field description of dynamic brittle fracture. Comput Methods Appl Mech Eng 217:77–9510.1016/j.cma.2012.01.008 ↗
13. Bourdin B, Francfort GA, Marigo J-J (2000) Numerical experiments in revisited brittle fracture. J Mech Phys Solids 48:797–82610.1016/S0022-5096(99)00028-9 ↗
14. Bourdin B, Marigo J-J, Maurini C, Sicsic P (2014) Morphogenesis and propagation of complex cracks induced by thermal shocks. Phys Rev Lett 112:01430110.1103/PhysRevLett.112.014301 ↗
15. Boz A, Demir O (2024) Fracture analysis for nozzle cracks in nuclear reactor pressure vessel using FCPAS. Nucl Eng Technol 56:2292–230610.1016/j.net.2024.01.040 ↗
16. Bui HP, Tomar S, Courtecuisse H, Cotin S, Bordas SP (2017) Real-time error control for surgical simulation. IEEE Trans Biomed Eng 65:596–60710.1109/TBME.2017.2695587 ↗
17. Chen W-X, Wu J-Y (2022) Phase-field cohesive zone modeling of multi-physical fracture in solids and the open-source implementation in Comsol Multiphysics. Theoret Appl Fract Mech 117:10315310.1016/j.tafmec.2021.103153 ↗
18. Chu DY, Li X, Liu ZL (2017) Study the dynamic crack path in brittle material under thermal shock loading by phase field modeling. Int J Fract 208:115–13010.1007/s10704-017-0220-4 ↗
19. De Vaucorbeil A, Nguyen VP, Sinaie S, Wu JY (2020) Material point method after 25 years: Theory, implementation, and applications. Adv Appl Mech 53:185–39810.1016/bs.aams.2019.11.001 ↗
20. Ding JL, Yu TT, Fang WH, Natarajan S (2024) An adaptive phase field modeling of fatigue crack growth using variable-node elements and explicit cycle jump scheme. Comput Methods Appl Mech Eng 429:11720010.1016/j.cma.2024.117200 ↗
21. Duflot M, Bordas S (2008) A posteriori error estimation for extended finite elements by an extended global recovery. Int J Numer Meth Eng 76:1123–113810.1002/nme.2332 ↗
22. Francfort GA, Marigo J-J (1998) Revisiting brittle fracture as an energy minimization problem. J Mech Phys Solids 46:1319–134210.1016/S0022-5096(98)00034-9 ↗
23. González-Estrada OA, Ródenas JJ, Bordas SP, Nadal E, Kerfriden P, Fuenmayor FJ (2015) Locally equilibrated stress recovery for goal oriented error estimation in the extended finite element method. Computers & Structures 152:1–1010.1016/j.compstruc.2015.01.015 ↗
24. Gupta A, Krishnan UM, Mandal TK, Chowdhury R, Nguyen VP (2022) An adaptive mesh refinement algorithm for phase-field fracture models: Application to brittle, cohesive, and dynamic fracture. Comput Methods Appl Mech Eng 399:11534710.1016/j.cma.2022.115347 ↗
25. He J-N, Yu TT, Hirshikesh (2025) Adaptive phase-field simulation of crack propagation in brittle materials: COMSOL implementation, Theoretical and Applied Fracture Mechanics , 10507710.1016/j.tafmec.2025.105077 ↗
26. He J-N, Kumaresan T, Yu TT, Fang WH, Natarajan S (2024) Numerical implementation and comparison study on simulating thermo-elastic fracture using adaptive phase-field method combined with BFGS algorithm and AM algorithm. Theoret Appl Fract Mech 133:10465010.1016/j.tafmec.2024.104650 ↗
27. He JN, Yu TT, Fang WH, Natarajan S (2024) Dynamic crack growth in orthotropic brittle materials using an adaptive phase-field modeling with variable-node elements. Compos Struct 337:11806810.1016/j.compstruct.2024.118068 ↗
28. He J-N, Yu TT, Fang WH, Natarajan S (2024) An adaptive dynamic phase-field modeling with variable-node elements for thermoelastic fracture in orthotropic media. Theoret Appl Fract Mech 133:10455510.1016/j.tafmec.2024.104555 ↗
29. Hirshikesh H, Pramod A, Waisman H, Natarajan S (2021) Adaptive phase field method using novel physics based refinement criteria. Comput Methods Appl Mech Eng 383:11387410.1016/j.cma.2021.113874 ↗
30. Honda S, Suzuki T, Nishikawa T, Awaji H, Akimune Y, Hirosaki N (2002) Estimation of thermal shock properties for silicon nitride having high thermal conductivity. J Ceram Soc Jpn 110:38–4310.2109/jcersj.110.38 ↗
31. Jin Y, González-Estrada O, Pierard O, Bordas S (2017) Error-controlled adaptive extended finite element method for 3D linear elastic crack propagation. Comput Methods Appl Mech Eng 318:319–34810.1016/j.cma.2016.12.016 ↗
32. Khelejani AS, Jahanshahi A, Shamshirsaz M (2025) Experimental and numerical analysis of cyclic ratchetting and low cycle fatigue behaviour in meander-shaped interconnects for stretchable electronics applications. Microelectron Reliab 167:11563210.1016/j.microrel.2025.115632 ↗
33. Krishnan UM, Gupta A, Chowdhury R (2022) Adaptive phase-field modeling of brittle fracture using a robust combination of error-estimator and markers. Eng Fract Mech 274:10875810.1016/j.engfracmech.2022.108758 ↗
34. Kumar A, Bourdin B, Francfort GA, Lopez-Pamies O (2020) Revisiting nucleation in the phase-field approach to brittle fracture. J Mech Phys Solids 142:10402710.1016/j.jmps.2020.104027 ↗
35. Lancioni G, Royer-Carfagni G (2009) The variational approach to fracture mechanics. A practical application to the French Panthéon in Paris Journal of Elasticity 95:1–30
36. Li WD, Ambati M, Nguyen-Thanh N, Du H, Zhou K (2023) Adaptive fourth-order phase-field modeling of ductile fracture using an isogeometric-meshfree approach. Comput Methods Appl Mech Eng 406:11586110.1016/j.cma.2022.115861 ↗
37. Li PD, Li WD, Fan HD, Wang QY, Zhou K (2023) A phase-field framework for brittle fracture in quasi-crystals. Int J Solids Struct 279:11238510.1016/j.ijsolstr.2023.112385 ↗
38. Li DY, Li PD, Li WD, Li WG, Zhou K (2022) Three-dimensional phase-field modeling of temperature-dependent thermal shock-induced fracture in ceramic materials. Eng Fract Mech 268:10844410.1016/j.engfracmech.2022.108444 ↗
39. Li WD, Nguyen-Thanh N, Zhou K (2020) Phase-field modeling of brittle fracture in a 3D polycrystalline material via an adaptive isogeometric-meshfree approach. Int J Numer Meth Eng 121:5042–506510.1002/nme.6509 ↗
40. Li DY, Pang Y, Lu T, Liu ZS, Chen SW (2022) Numerical analysis of thermal shock cracking behaviors of ceramics based on the force-heat equivalence energy density principle. Frontiers in Materials 8:82532710.3389/fmats.2021.825327 ↗
41. Liu Y, Zhao HT, Liu K, Zhao ZJ, Feng M, Chen J (2025) Multi-scale modeling: Analysis and design of thermal-mechanical coupling behavior of integrated thermal protection systems, Chinese Journal of Aeronautics , 10385510.1016/j.cja.2025.103855 ↗
42. Li YC, Yu TT, Xing C, Natarajan S (2023) Modeling quasi-static and dynamic thermo-elastic coupled brittle fracture using an adaptive isogeometric hybrid phase-field method. Finite Elem Anal Des 224:10399310.1016/j.finel.2023.103993 ↗
43. Mandal TK, Nguyen VP, Wu J-Y (2019) Length scale and mesh bias sensitivity of phase-field models for brittle and cohesive fracture. Eng Fract Mech 217:10653210.1016/j.engfracmech.2019.106532 ↗
44. Mandal TK, Nguyen VP, Wu J-Y, Nguyen-Thanh C, de Vaucorbeil A (2021) Fracture of thermo-elastic solids: Phase-field modeling and new results with an efficient monolithic solver. Comput Methods Appl Mech Eng 376:11364810.1016/j.cma.2020.113648 ↗
45. Manson S, Smith R (1955) Theory of thermal shock resistance of brittle materials based on Weibull’s statistical theory of strength. J Am Ceram Soc 38:18–2710.1111/j.1151-2916.1955.tb14547.x ↗
46. Melenk JM, Babuška I (1996) The partition of unity finite element method: basic theory and applications. Comput Methods Appl Mech Eng 139:289–31410.1016/S0045-7825(96)01087-0 ↗
47. Meng L, Lee HW, Ashkpour A, Sales CM, Farnam YA, Hubler MH, Najafi AR (2024) Bound-constrained optimization using Lagrange multiplier for a length scale insensitive phase field fracture model. Eng Fract Mech 310:11049610.1016/j.engfracmech.2024.110496 ↗
48. Menk A, Bordas SP (2010) Numerically determined enrichment functions for the extended finite element method and applications to bi-material anisotropic fracture and polycrystals. Int J Numer Meth Eng 83:805–82810.1002/nme.2858 ↗
49. Menk A, Bordas SP (2011) Crack growth calculations in solder joints based on microstructural phenomena with X-FEM. Comput Mater Sci 50:1145–115610.1016/j.commatsci.2010.11.014 ↗
50. Miehe C, Hofacker M, Welschinger F (2010) A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits. Comput Methods Appl Mech Eng 199:2765–277810.1016/j.cma.2010.04.011 ↗
51. Miehe C, Schaenzel LM, Ulmer H (2015) Phase field modeling of fracture in multi-physics problems. Part I. Balance of crack surface and failure criteria for brittle crack propagation in thermo-elastic solids Computer Methods in Applied Mechanics and Engineering 294:449–485
52. Nguyen V-H, Trinh M-C, Jun H (2024) Fracture behavior of thermal mismatch in functionally graded materials using phase-field modeling. Eng Fract Mech 310:11043310.1016/j.engfracmech.2024.110433 ↗
53. Nie WJ, Yang XG, Zhang KY, Li J, Zhang QC, Yuan W (2025) Test rig and experimental Investigation on the blade containment of aero-engine turbine casing, Engineering Failure Analysis , 10972610.1016/j.engfailanal.2025.109726 ↗
54. Notay Y (2012) Aggregation-based algebraic multigrid for convection-diffusion equations. SIAM J Sci Comput 34:A2288–A231610.1137/110835347 ↗
55. Olson LN, Schroder JB (2011) Smoothed aggregation multigrid solvers for high-order discontinuous Galerkin methods for elliptic problems. J Comput Phys 230:6959–697610.1016/j.jcp.2011.05.009 ↗
56. Panetier J, Ladeveze P, Chamoin L (2010) Strict and effective bounds in goal-oriented error estimation applied to fracture mechanics problems solved with XFEM. Int J Numer Meth Eng 81:671–70010.1002/nme.2705 ↗
57. Pang Y, Li PD, Li DY, Zhou XD, Fan HD, Wang QY (2024) Phase-field modeling of thermal shock fracture in functionally graded materials. Eng Fract Mech 307:11028610.1016/j.engfracmech.2024.110286 ↗
58. Papšík R, Ševeček O, Schlacher J, Bermejo R (2024) Prediction of thermal shock induced cracking in multi-material ceramics using a stress-energy criterion. Eng Fract Mech 303:11012110.1016/j.engfracmech.2024.110121 ↗
59. Piska R, El-Borgi S, Nafees M, Rajagopal A, Reddy J (2024) A thermodynamically consistent phase field model for brittle fracture in graded coatings under thermo-mechanical loading. Theoret Appl Fract Mech 131:10441410.1016/j.tafmec.2024.104414 ↗
60. Remmers JJ, Borst Rd, Needleman A (2003) A cohesive segments method for the simulation of crack growth. Comput Mech 31:69–7710.1007/s00466-002-0394-z ↗
61. Ricardo LFF, Leguillon D, Parry G, Doitrand A (2020) Modeling the thermal shock induced cracking in ceramics. J Eur Ceram Soc 40:1513–152110.1016/j.jeurceramsoc.2019.11.071 ↗
62. Sicsic P, Marigo J-J, Maurini C (2014) Initiation of a periodic array of cracks in the thermal shock problem: a gradient damage modeling. J Mech Phys Solids 63:256–28410.1016/j.jmps.2013.09.003 ↗
63. Sołowski WT, Berzins M, Coombs WM, Guilkey JE, Möller M, Tran QA, Adibaskoro T, Seyedan S, Tielen R, Soga K (2021) Material point method: overview and challenges ahead. Adv Appl Mech 54:113–20410.1016/bs.aams.2020.12.002 ↗
64. Svolos L, Bronkhorst CA, Waisman H (2020) Thermal-conductivity degradation across cracks in coupled thermo-mechanical systems modeled by the phase-field fracture method. J Mech Phys Solids 137:10386110.1016/j.jmps.2019.103861 ↗
65. Svolos L, Mourad HM, Bronkhorst CA, Waisman H (2021) Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality. Eng Fract Mech 245:10755410.1016/j.engfracmech.2021.107554 ↗
66. Tanné E, Li TY, Bourdin B, Marigo J-J, Maurini C (2018) Crack nucleation in variational phase-field models of brittle fracture. J Mech Phys Solids 110:80–9910.1016/j.jmps.2017.09.006 ↗
67. Tian R, Wen LF (2015) Improved XFEM-An extra-dof free, well-conditioning, and interpolating XFEM. Comput Methods Appl Mech Eng 285:639–65810.1016/j.cma.2014.11.026 ↗
68. Wang T, Han HY, Wang YF, Ye X, Huang GY, Liu ZL, Zhuang Z (2022) Simulation of crack patterns in quasi-brittle materials under thermal shock using phase field and cohesive zone models. Eng Fract Mech 276:10888910.1016/j.engfracmech.2022.108889 ↗
69. Wang T, Zhang YC, Han HY, Wang L, Ye X, Zhuang Z (2024) Phase field modeling of crack propagation in three-dimensional quasi-brittle materials under thermal shock. Eng Fract Mech 302:11007010.1016/j.engfracmech.2024.110070 ↗
70. Wang YT, Zhou XP, Kou MM (2019) An improved coupled thermo-mechanic bond-based peridynamic model for cracking behaviors in brittle solids subjected to thermal shocks. European Journal of Mechanics-A/Solids 73:282–30510.1016/j.euromechsol.2018.09.007 ↗
71. Wu J-Y, Cervera M (2018) A novel positive/negative projection in energy norm for the damage modeling of quasi-brittle solids. Int J Solids Struct 139:250–26910.1016/j.ijsolstr.2018.02.004 ↗
72. Wu J-Y, Nguyen VP, Zhou H, Huang YL (2020) A variationally consistent phase-field anisotropic damage model for fracture. Comput Methods Appl Mech Eng 358:11262910.1016/j.cma.2019.112629 ↗
73. Xia CC, Xu YJ, Zhou SW, Qin SK, He XH (2023) Fracture initiation and propagation in the lined underground caverns for compressed air energy storage: Coupled thermo-mechanical phase-field modeling. Comput Geotech 157:10532910.1016/j.compgeo.2023.105329 ↗
74. Yue Q, Wang Q, Rabczuk T, Zhou W, Zhuang XY, Chang XL (2024) A thermo-mechanical phase-field model for mixed-mode fracture and its application in rock-like materials. Int J Rock Mech Min Sci 183:10590710.1016/j.ijrmms.2024.105907 ↗
75. Yu YF, Hou C (2025) A space-time adaptive solution scheme for phase field fracture model. Computers & Structures 316:10790410.1016/j.compstruc.2025.107904 ↗
76. Zhang TC, Hirshikesh, Yu TT, Ding JL, Natarajan S (2024) Low-cycle fatigue crack growth in brittle materials: Adaptive phase-field modeling with variable-node elements, Computer Methods in Applied Mechanics and Engineering 425 , 11691710.1016/j.cma.2024.116917 ↗
77. Zhang X, Chen Z, Liu Y (2016) The material point method: a continuum-based particle method for extreme loading cases, Academic Press,10.1016/B978-0-12-407716-4.00003-X ↗
78. Zhao SM, Li PD, Wang T, Tan Y, Fan HD, Wang QY (2023) A phase-field model for thermo-elastic fracture in quasicrystals. Eng Fract Mech 289:10943210.1016/j.engfracmech.2023.109432 ↗
79. Zhou SW, Rabczuk T, Zhuang XY (2018) Phase field modeling of quasi-static and dynamic crack propagation: COMSOL implementation and case studies. Adv Eng Softw 122:31–4910.1016/j.advengsoft.2018.03.012 ↗
80. Zhou SW, Zhuang XY, Rabczuk T (2019) Phase-field modeling of fluid-driven dynamic cracking in porous media. Comput Methods Appl Mech Eng 350:169–19810.1016/j.cma.2019.03.001 ↗
81. Zou QM (2023) GMRES algorithms over 35 years. Appl Math Comput 445:12786910.1016/j.amc.2023.127869 ↗