DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials

Abstract

This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.

Authors:
 [1]; ORCiD logo [1]
  1. Columbia Univ., New York, NY (United States)
Publication Date:
Research Org.:
Columbia Univ., New York, NY (United States)
Sponsoring Org.:
USDOE Office of Nuclear Energy (NE), Nuclear Energy University Program (NEUP); USDOE
OSTI Identifier:
1801240
Alternate Identifier(s):
OSTI ID: 1592511
Grant/Contract Number:  
NE0008534
Resource Type:
Accepted Manuscript
Journal Name:
Computer Methods in Applied Mechanics and Engineering
Additional Journal Information:
Journal Volume: 362; Journal ID: ISSN 0045-7825
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
42 ENGINEERING; fast Fourier transform; multi-phase-field fracture; higher-order phase field fracture; anisotropic cracks; polycrystal; salt

Citation Formats

Ma, Ran, and Sun, WaiChing. FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials. United States: N. p., 2020. Web. doi:10.1016/j.cma.2019.112781.
Ma, Ran, & Sun, WaiChing. FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials. United States. https://doi.org/10.1016/j.cma.2019.112781
Ma, Ran, and Sun, WaiChing. Mon . "FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials". United States. https://doi.org/10.1016/j.cma.2019.112781. https://www.osti.gov/servlets/purl/1801240.
@article{osti_1801240,
title = {FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials},
author = {Ma, Ran and Sun, WaiChing},
abstractNote = {This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.},
doi = {10.1016/j.cma.2019.112781},
journal = {Computer Methods in Applied Mechanics and Engineering},
number = ,
volume = 362,
place = {United States},
year = {Mon Jan 20 00:00:00 EST 2020},
month = {Mon Jan 20 00:00:00 EST 2020}
}

Journal Article:

Citation Metrics:
Cited by: 28 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

A fast numerical scheme for computing the response of composites using grid refinement
journal, April 1999

  • Eyre, D. J.; Milton, G. W.
  • The European Physical Journal Applied Physics, Vol. 6, Issue 1
  • DOI: 10.1051/epjap:1999150

Computational thermomechanics of crystalline rock, Part I: A combined multi-phase-field/crystal plasticity approach for single crystal simulations
journal, August 2018


A unified variational eigen-erosion framework for interacting brittle fractures and compaction bands in fluid-infiltrating porous media
journal, May 2017


Transition between crack patterns in quenched glass plates
journal, March 1993


Meta-modeling game for deriving theory-consistent, microstructure-based traction–separation laws via deep reinforcement learning
journal, April 2019


Spectral database constitutive representation within a spectral micromechanical solver for computationally efficient polycrystal plasticity modelling
journal, May 2017

  • Eghtesad, Adnan; Zecevic, Miroslav; Lebensohn, Ricardo A.
  • Computational Mechanics, Vol. 61, Issue 1-2
  • DOI: 10.1007/s00466-017-1413-4

A phase-field formulation for fracture in ductile materials: Finite deformation balance law derivation, plastic degradation, and stress triaxiality effects
journal, December 2016

  • Borden, Michael J.; Hughes, Thomas J. R.; Landis, Chad M.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 312
  • DOI: 10.1016/j.cma.2016.09.005

Isogeometric analysis of the Cahn–Hilliard phase-field model
journal, September 2008

  • Gómez, Héctor; Calo, Victor M.; Bazilevs, Yuri
  • Computer Methods in Applied Mechanics and Engineering, Vol. 197, Issue 49-50
  • DOI: 10.1016/j.cma.2008.05.003

Effects of spatial heterogeneity and material anisotropy on the fracture pattern and macroscopic effective toughness of Mancos Shale in Brazilian tests
journal, August 2017

  • Na, SeonHong; Sun, WaiChing; Ingraham, Mathew D.
  • Journal of Geophysical Research: Solid Earth, Vol. 122, Issue 8
  • DOI: 10.1002/2016JB013374

A FFT solver for variational phase-field modeling of brittle fracture
journal, June 2019

  • Chen, Yang; Vasiukov, Dmytro; Gélébart, Lionel
  • Computer Methods in Applied Mechanics and Engineering, Vol. 349
  • DOI: 10.1016/j.cma.2019.02.017

A numerical method for computing the overall response of nonlinear composites with complex microstructure
journal, April 1998


Phase field modeling of ductile fracture at finite strains: A variational gradient-extended plasticity-damage theory
journal, September 2016


Phase-field modeling and simulation of fracture in brittle materials with strongly anisotropic surface energy: PHASE-FIELD MODELING AND SIMULATION OF STRONGLY ANISOTROPIC FRACTURE
journal, July 2014

  • Li, Bin; Peco, Christian; Millán, Daniel
  • International Journal for Numerical Methods in Engineering, Vol. 102, Issue 3-4
  • DOI: 10.1002/nme.4726

A higher-order phase-field model for brittle fracture: Formulation and analysis within the isogeometric analysis framework
journal, May 2014

  • Borden, Michael J.; Hughes, Thomas J. R.; Landis, Chad M.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 273
  • DOI: 10.1016/j.cma.2014.01.016

A mixed-mode phase field fracture model in anisotropic rocks with consistent kinematics
journal, December 2018

  • Bryant, Eric C.; Sun, WaiChing
  • Computer Methods in Applied Mechanics and Engineering, Vol. 342
  • DOI: 10.1016/j.cma.2018.08.008

A multiscale approach for modeling progressive damage of composite materials using fast Fourier transforms
journal, January 2014

  • Spahn, J.; Andrä, H.; Kabel, M.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 268
  • DOI: 10.1016/j.cma.2013.10.017

Forbidden Directions for the Fracture of Thin Anisotropic Sheets: An Analogy with the Wulff Plot
journal, April 2013


A strong discontinuity approach on multiple levels to model solids at failure
journal, January 2013

  • Linder, Christian; Raina, Arun
  • Computer Methods in Applied Mechanics and Engineering, Vol. 253
  • DOI: 10.1016/j.cma.2012.07.005

Simulation of micromechanical behavior of polycrystals: finite elements versus fast Fourier transforms
journal, August 2009

  • Prakash, A.; Lebensohn, R. A.
  • Modelling and Simulation in Materials Science and Engineering, Vol. 17, Issue 6
  • DOI: 10.1088/0965-0393/17/6/064010

Experimental validation of a phase-field model for fracture
journal, January 2017

  • Pham, K. H.; Ravi-Chandar, K.; Landis, C. M.
  • International Journal of Fracture, Vol. 205, Issue 1
  • DOI: 10.1007/s10704-017-0185-3

Revisiting brittle fracture as an energy minimization problem
journal, August 1998


An updated Lagrangian LBM–DEM–FEM coupling model for dual-permeability fissured porous media with embedded discontinuities
journal, February 2019


Sparse and scalable eigenstrain-based reduced order homogenization models for polycrystal plasticity
journal, November 2017


A nonlocal multiscale discrete-continuum model for predicting mechanical behavior of granular materials: A NONLOCAL MULTISCALE DISCRETE-CONTINUUM MODEL FOR GRANULAR MATERIALS
journal, November 2015

  • Liu, Yang; Sun, WaiChing; Yuan, Zifeng
  • International Journal for Numerical Methods in Engineering, Vol. 106, Issue 2
  • DOI: 10.1002/nme.5139

Phase field modeling of directional fracture in anisotropic polycrystals
journal, February 2015


Phase-field modeling of anisotropic brittle fracture including several damage mechanisms
journal, July 2018

  • Bleyer, Jeremy; Alessi, Roberto
  • Computer Methods in Applied Mechanics and Engineering, Vol. 336
  • DOI: 10.1016/j.cma.2018.03.012

Thermodynamically consistent phase-field models of fracture: Variational principles and multi-field FE implementations
journal, August 2010

  • Miehe, C.; Welschinger, F.; Hofacker, M.
  • International Journal for Numerical Methods in Engineering, Vol. 83, Issue 10
  • DOI: 10.1002/nme.2861

Fourier-based schemes with modified Green operator for computing the electrical response of heterogeneous media with accurate local fields: FOURIER-BASED SCHEMES WITH MODIFIED GREEN OPERATOR
journal, February 2014

  • Willot, François; Abdallah, Bassam; Pellegrini, Yves-Patrick
  • International Journal for Numerical Methods in Engineering, Vol. 98, Issue 7
  • DOI: 10.1002/nme.4641

Crack kinking in a variational phase-field model of brittle fracture with strongly anisotropic surface energy
journal, April 2019


Multi-phase-field modeling of anisotropic crack propagation for polycrystalline materials
journal, April 2017

  • Nguyen, Thanh-Tung; Réthoré, Julien; Yvonnet, Julien
  • Computational Mechanics, Vol. 60, Issue 2
  • DOI: 10.1007/s00466-017-1409-0

Progressive damage analysis of 3D braided composites using FFT-based method
journal, May 2018


Applications of semi-implicit Fourier-spectral method to phase field equations
journal, February 1998


Necessary and sufficient elastic stability conditions in various crystal systems
journal, December 2014


Large-scale 3D random polycrystals for the finite element method: Generation, meshing and remeshing
journal, April 2011

  • Quey, R.; Dawson, P. R.; Barbe, F.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 200, Issue 17-20
  • DOI: 10.1016/j.cma.2011.01.002

Finite-deformation irreversible cohesive elements for three-dimensional crack-propagation analysis
journal, March 1999


Fracture of solar-grade anisotropic polycrystalline Silicon: A combined phase field–cohesive zone model approach
journal, March 2018

  • Paggi, M.; Corrado, M.; Reinoso, J.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 330
  • DOI: 10.1016/j.cma.2017.10.021

A multiscale mean-field homogenization method for fiber-reinforced composites with gradient-enhanced damage models
journal, August 2012

  • Wu, L.; Noels, L.; Adam, L.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 233-236
  • DOI: 10.1016/j.cma.2012.04.011

A multiscale multi-permeability poroplasticity model linked by recursive homogenizations and deep learning
journal, June 2018


Multiscale method for characterization of porous microstructures and their impact on macroscopic effective permeability
journal, May 2011

  • Sun, W. C.; Andrade, J. E.; Rudnicki, J. W.
  • International Journal for Numerical Methods in Engineering, Vol. 88, Issue 12
  • DOI: 10.1002/nme.3220

The Variational Approach to Fracture
journal, March 2008

  • Bourdin, Blaise; Francfort, Gilles A.; Marigo, Jean-Jacques
  • Journal of Elasticity, Vol. 91, Issue 1-3
  • DOI: 10.1007/s10659-007-9107-3

A finite element perspective on nonlinear FFT-based micromechanical simulations: A FINITE ELEMENT PERSPECTIVE ON NONLINEAR FFT-BASED SIMULATIONS
journal, March 2017

  • Zeman, J.; de Geus, T. W. J.; Vondřejc, J.
  • International Journal for Numerical Methods in Engineering, Vol. 111, Issue 10
  • DOI: 10.1002/nme.5481

A finite element method for crack growth without remeshing
journal, September 1999


Two-scale FE–FFT- and phase-field-based computational modeling of bulk microstructural evolution and macroscopic material behavior
journal, June 2016

  • Kochmann, Julian; Wulfinghoff, Stephan; Reese, Stefanie
  • Computer Methods in Applied Mechanics and Engineering, Vol. 305
  • DOI: 10.1016/j.cma.2016.03.001

Computational homogenization of elasticity on a staggered grid: COMPUTATIONAL HOMOGENIZATION OF ELASTICITY ON A STAGGERED GRID
journal, October 2015

  • Schneider, Matti; Ospald, Felix; Kabel, Matthias
  • International Journal for Numerical Methods in Engineering, Vol. 105, Issue 9
  • DOI: 10.1002/nme.5008

A damage model for crack prediction in brittle and quasi-brittle materials solved by the FFT method
journal, January 2012

  • Li, Jia; Tian, Xiao-Xiao; Abdelmoula, Radhi
  • International Journal of Fracture, Vol. 173, Issue 2
  • DOI: 10.1007/s10704-011-9671-1

A study of the influence of REV variability in double-scale FEM ×DEM analysis: A study of the influence of REV variability in double-scale FEM ×DEM analysis
journal, February 2016

  • Shahin, Ghassan; Desrues, Jacques; Pont, Stefano Dal
  • International Journal for Numerical Methods in Engineering, Vol. 107, Issue 10
  • DOI: 10.1002/nme.5202

A phase-field model for cohesive fracture: A PHASE-FIELD MODEL FOR COHESIVE FRACTURE
journal, July 2013

  • Verhoosel, Clemens V.; de Borst, René
  • International Journal for Numerical Methods in Engineering, Vol. 96, Issue 1
  • DOI: 10.1002/nme.4553

Ab initio study of the elastic properties of sodium chloride at high pressure
journal, May 2010


Plane strain finite element analysis of pressure sensitive plasticity with strong discontinuity
journal, May 2001


FFT-based homogenization of hypoelastic plasticity at finite strains
journal, June 2019

  • Ma, Ran; Truster, Timothy J.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 349
  • DOI: 10.1016/j.cma.2019.02.037

Pressure and fluid-driven fracture propagation in porous media using an adaptive finite element phase field model
journal, June 2016

  • Lee, Sanghyun; Wheeler, Mary F.; Wick, Thomas
  • Computer Methods in Applied Mechanics and Engineering, Vol. 305
  • DOI: 10.1016/j.cma.2016.02.037

An integrated fast Fourier transform-based phase-field and crystal plasticity approach to model recrystallization of three dimensional polycrystals
journal, March 2015

  • Chen, L.; Chen, J.; Lebensohn, R. A.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 285
  • DOI: 10.1016/j.cma.2014.12.007

An augmented-Lagrangian method for the phase-field approach for pressurized fractures
journal, April 2014

  • Wheeler, M. F.; Wick, T.; Wollner, W.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 271
  • DOI: 10.1016/j.cma.2013.12.005

A phase field model of deformation twinning: Nonlinear theory and numerical simulations
journal, April 2011


Generalized finite element methods for three-dimensional structural mechanics problems
journal, June 2000


FFT-based interface decohesion modelling by a nonlocal interphase
journal, April 2018

  • Sharma, Luv; Peerlings, Ron H. J.; Shanthraj, Pratheek
  • Advanced Modeling and Simulation in Engineering Sciences, Vol. 5, Issue 1
  • DOI: 10.1186/s40323-018-0100-0

On micro-cracking, inelastic dilatancy, and the brittle-ductile transition in compact rocks: A micro-mechanical study
journal, May 2008


Coupled phase-field and plasticity modeling of geological materials: From brittle fracture to ductile flow
journal, March 2018

  • Choo, Jinhyun; Sun, WaiChing
  • Computer Methods in Applied Mechanics and Engineering, Vol. 330
  • DOI: 10.1016/j.cma.2017.10.009

A phase field approach to simulate intralaminar and translaminar fracture in long fiber composite materials
journal, July 2019


A review on phase-field models of brittle fracture and a new fast hybrid formulation
journal, December 2014

  • Ambati, Marreddy; Gerasimov, Tymofiy; De Lorenzis, Laura
  • Computational Mechanics, Vol. 55, Issue 2
  • DOI: 10.1007/s00466-014-1109-y

Cracking and damage from crystallization in pores: Coupled chemo-hydro-mechanics and phase-field modeling
journal, June 2018

  • Choo, Jinhyun; Sun, WaiChing
  • Computer Methods in Applied Mechanics and Engineering, Vol. 335
  • DOI: 10.1016/j.cma.2018.01.044

A multiscale DEM-LBM analysis on permeability evolutions inside a dilatant shear band
journal, May 2013


An optimization-based phase-field method for continuous-discontinuous crack propagation: A phase-field method for continuous-discontinuous crack propagation
journal, July 2018

  • Geelen, Rudy J. M.; Liu, Yingjie; Dolbow, John E.
  • International Journal for Numerical Methods in Engineering, Vol. 116, Issue 1
  • DOI: 10.1002/nme.5911

The role of interfaces in enhancing the yield strength of composites and polycrystals
journal, May 2005


Formulation of nonlocal damage models based on spectral methods for application to complex microstructures
journal, October 2015


Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients
journal, October 2010

  • Zeman, Jan; Vondřejc, Jaroslav; Novák, Jan
  • Journal of Computational Physics, Vol. 229, Issue 21
  • DOI: 10.1016/j.jcp.2010.07.010

Mixed Arlequin method for multiscale poromechanics problems: Mixed Arlequin method for multiscale poromechanics problems
journal, February 2017

  • Sun, WaiChing; Cai, Zhijun; Choo, Jinhyun
  • International Journal for Numerical Methods in Engineering, Vol. 111, Issue 7
  • DOI: 10.1002/nme.5476

N-site modeling of a 3D viscoplastic polycrystal using Fast Fourier Transform
journal, August 2001


Phase field modeling of fracture in anisotropic brittle solids
journal, December 2017


A numerical spectral approach for solving elasto-static field dislocation and g-disclination mechanics
journal, November 2014

  • Berbenni, Stéphane; Taupin, Vincent; Djaka, Komlan Sénam
  • International Journal of Solids and Structures, Vol. 51, Issue 23-24
  • DOI: 10.1016/j.ijsolstr.2014.08.009

Aspects of the computational testing of the mechanical properties of microheterogeneous material samples
journal, January 2001

  • Zohdi, T. I.; Wriggers, P.
  • International Journal for Numerical Methods in Engineering, Vol. 50, Issue 11
  • DOI: 10.1002/nme.146

Weakening of rock salt by water during long-term creep
journal, December 1986

  • Urai, Janos L.; Spiers, Christopher J.; Zwart, Hendrik J.
  • Nature, Vol. 324, Issue 6097
  • DOI: 10.1038/324554a0

A multilevel finite element method (FE2) to describe the response of highly non-linear structures using generalized continua
journal, July 2003


Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement
journal, October 2005

  • Hughes, T. J. R.; Cottrell, J. A.; Bazilevs, Y.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 194, Issue 39-41
  • DOI: 10.1016/j.cma.2004.10.008

Effective properties of composite materials with periodic microstructure: a computational approach
journal, April 1999

  • Michel, J. C.; Moulinec, H.; Suquet, P.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 172, Issue 1-4
  • DOI: 10.1016/S0045-7825(98)00227-8

A scalable 3D fracture and fragmentation algorithm based on a hybrid, discontinuous Galerkin, cohesive element method
journal, January 2011

  • Radovitzky, R.; Seagraves, A.; Tupek, M.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 200, Issue 1-4
  • DOI: 10.1016/j.cma.2010.08.014