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

Title: Unconditionally stable, second-order accurate schemes for solid state phase transformations driven by mechano-chemical spinodal decomposition

Abstract

Here, we consider solid state phase transformations that are caused by free energy densities with domains of non-convexity in strain-composition space; we refer to the non-convex domains as mechano-chemical spinodals. The non-convexity with respect to composition and strain causes segregation into phases with different crystal structures. We work on an existing model that couples the classical Cahn-Hilliard model with Toupin’s theory of gradient elasticity at finite strains. Both systems are represented by fourth-order, nonlinear, partial differential equations. The goal of this work is to develop unconditionally stable, second-order accurate time-integration schemes, motivated by the need to carry out large scale computations of dynamically evolving microstructures in three dimensions. We also introduce reduced formulations naturally derived from these proposed schemes for faster computations that are still second-order accurate. Although our method is developed and analyzed here for a specific class of mechano-chemical problems, one can readily apply the same method to develop unconditionally stable, second-order accurate schemes for any problems for which free energy density functions are multivariate polynomials of solution components and component gradients. Apart from an analysis and construction of methods, we present a suite of numerical results that demonstrate the schemes in action.

Authors:
 [1];  [1];  [1]
  1. Univ. of Michigan, Ann Arbor, MI (United States)
Publication Date:
Research Org.:
Univ. of Michigan, Ann Arbor, MI (United States). PRedictive Integrated Structural Materials Science (PRISMS) Center
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES)
OSTI Identifier:
1324034
Alternate Identifier(s):
OSTI ID: 1399034
Grant/Contract Number:  
SC0008637
Resource Type:
Accepted Manuscript
Journal Name:
Computer Methods in Applied Mechanics and Engineering
Additional Journal Information:
Journal Name: Computer Methods in Applied Mechanics and Engineering; Journal ID: ISSN 0045-7825
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; 36 MATERIALS SCIENCE; unconditional stability; mechano-chemistry; Cahn-Hilliard; strain gradient elasticity; spinodal decomposition; phase-field

Citation Formats

Sagiyama, Koki, Rudraraju, Shiva, and Garikipati, Krishna. Unconditionally stable, second-order accurate schemes for solid state phase transformations driven by mechano-chemical spinodal decomposition. United States: N. p., 2016. Web. doi:10.1016/j.cma.2016.09.003.
Sagiyama, Koki, Rudraraju, Shiva, & Garikipati, Krishna. Unconditionally stable, second-order accurate schemes for solid state phase transformations driven by mechano-chemical spinodal decomposition. United States. https://doi.org/10.1016/j.cma.2016.09.003
Sagiyama, Koki, Rudraraju, Shiva, and Garikipati, Krishna. Tue . "Unconditionally stable, second-order accurate schemes for solid state phase transformations driven by mechano-chemical spinodal decomposition". United States. https://doi.org/10.1016/j.cma.2016.09.003. https://www.osti.gov/servlets/purl/1324034.
@article{osti_1324034,
title = {Unconditionally stable, second-order accurate schemes for solid state phase transformations driven by mechano-chemical spinodal decomposition},
author = {Sagiyama, Koki and Rudraraju, Shiva and Garikipati, Krishna},
abstractNote = {Here, we consider solid state phase transformations that are caused by free energy densities with domains of non-convexity in strain-composition space; we refer to the non-convex domains as mechano-chemical spinodals. The non-convexity with respect to composition and strain causes segregation into phases with different crystal structures. We work on an existing model that couples the classical Cahn-Hilliard model with Toupin’s theory of gradient elasticity at finite strains. Both systems are represented by fourth-order, nonlinear, partial differential equations. The goal of this work is to develop unconditionally stable, second-order accurate time-integration schemes, motivated by the need to carry out large scale computations of dynamically evolving microstructures in three dimensions. We also introduce reduced formulations naturally derived from these proposed schemes for faster computations that are still second-order accurate. Although our method is developed and analyzed here for a specific class of mechano-chemical problems, one can readily apply the same method to develop unconditionally stable, second-order accurate schemes for any problems for which free energy density functions are multivariate polynomials of solution components and component gradients. Apart from an analysis and construction of methods, we present a suite of numerical results that demonstrate the schemes in action.},
doi = {10.1016/j.cma.2016.09.003},
journal = {Computer Methods in Applied Mechanics and Engineering},
number = ,
volume = ,
place = {United States},
year = {Tue Sep 13 00:00:00 EDT 2016},
month = {Tue Sep 13 00:00:00 EDT 2016}
}

Journal Article:

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

Save / Share:

Works referenced in this record:

Free Energy of a Nonuniform System. I. Interfacial Free Energy
journal, February 1958

  • Cahn, John W.; Hilliard, John E.
  • The Journal of Chemical Physics, Vol. 28, Issue 2
  • DOI: 10.1063/1.1744102

Theories of elasticity with couple-stress
journal, January 1964

  • Toupin, R. A.
  • Archive for Rational Mechanics and Analysis, Vol. 17, Issue 2
  • DOI: 10.1007/BF00253050

The Global Dynamics of Discrete Semilinear Parabolic Equations
journal, December 1993

  • Elliott, C. M.; Stuart, A. M.
  • SIAM Journal on Numerical Analysis, Vol. 30, Issue 6
  • DOI: 10.1137/0730084

Unconditionally Gradient Stable Time Marching the Cahn-Hilliard Equation
journal, January 1998


An Energy-Stable and Convergent Finite-Difference Scheme for the Phase Field Crystal Equation
journal, January 2009

  • Wise, S. M.; Wang, C.; Lowengrub, J. S.
  • SIAM Journal on Numerical Analysis, Vol. 47, Issue 3
  • DOI: 10.1137/080738143

Numerical schemes for a three component Cahn-Hilliard model
journal, December 2010

  • Boyer, Franck; Minjeaud, Sebastian
  • ESAIM: Mathematical Modelling and Numerical Analysis, Vol. 45, Issue 4
  • DOI: 10.1051/m2an/2010072

A gradient stable scheme for a phase field model for the moving contact line problem
journal, February 2012


Numerical Analysis of a Continuum Model of Phase Transition
journal, October 1991

  • Du, Qiang; Nicolaides, R. A.
  • SIAM Journal on Numerical Analysis, Vol. 28, Issue 5
  • DOI: 10.1137/0728069

Second order schemes and time-step adaptivity for Allen–Cahn and Cahn–Hilliard models
journal, October 2014

  • Guillén-González, Francisco; Tierra, Giordano
  • Computers & Mathematics with Applications, Vol. 68, Issue 8
  • DOI: 10.1016/j.camwa.2014.07.014

Fully Discrete Finite Element Approximations of the Navier--Stokes--Cahn-Hilliard Diffuse Interface Model for Two-Phase Fluid Flows
journal, January 2006

  • Feng, Xiaobing
  • SIAM Journal on Numerical Analysis, Vol. 44, Issue 3
  • DOI: 10.1137/050638333

Energy law preserving C0 finite element schemes for phase field models in two-phase flow computations
journal, August 2011


An energy law preserving C0 finite element scheme for simulating the kinematic effects in liquid crystal dynamics
journal, December 2007


Energetic variational approach in complex fluids: Maximum dissipation principle
journal, January 2010

  • Hyon, Yunkyong; Young Kwak, Do; Liu, Chun
  • Discrete & Continuous Dynamical Systems - A, Vol. 26, Issue 4
  • DOI: 10.3934/dcds.2010.26.1291

Conservative multigrid methods for Cahn–Hilliard fluids
journal, January 2004

  • Kim, Junseok; Kang, Kyungkeun; Lowengrub, John
  • Journal of Computational Physics, Vol. 193, Issue 2
  • DOI: 10.1016/j.jcp.2003.07.035

Conservative multigrid methods for ternary Cahn-Hilliard systems
journal, January 2004

  • Kang, Kyungkeun; Kim, Junseok; Lowengrub, John
  • Communications in Mathematical Sciences, Vol. 2, Issue 1
  • DOI: 10.4310/CMS.2004.v2.n1.a4

Stabilized second-order convex splitting schemes for Cahn-Hilliard models with application to diffuse-interface tumor-growth models: STABILIZED SECOND-ORDER CONVEX SPLITTING SCHEMES FOR CAHN-HILLIARD MODELS
journal, September 2013

  • Wu, X.; van Zwieten, G. J.; van der Zee, K. G.
  • International Journal for Numerical Methods in Biomedical Engineering, Vol. 30, Issue 2
  • DOI: 10.1002/cnm.2597

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


A discontinuous Galerkin method for the Cahn–Hilliard equation
journal, November 2006

  • Wells, Garth N.; Kuhl, Ellen; Garikipati, Krishna
  • Journal of Computational Physics, Vol. 218, Issue 2
  • DOI: 10.1016/j.jcp.2006.03.010

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

Three-dimensional isogeometric solutions to general boundary value problems of Toupin’s gradient elasticity theory at finite strains
journal, August 2014

  • Rudraraju, S.; Van der Ven, A.; Garikipati, K.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 278
  • DOI: 10.1016/j.cma.2014.06.015

A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening
journal, June 1979


Three-dimensional isogeometric solutions to general boundary value problems of Toupin’s gradient elasticity theory at finite strains
journal, August 2014

  • Rudraraju, S.; Van der Ven, A.; Garikipati, K.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 278
  • DOI: 10.1016/j.cma.2014.06.015

Mechanochemical spinodal decomposition: a phenomenological theory of phase transformations in multi-component, crystalline solids
journal, June 2016


An energy-stable convex splitting for the phase-field crystal equation
text, January 2014


Works referencing / citing this record:

PRISMS: An Integrated, Open-Source Framework for Accelerating Predictive Structural Materials Science
journal, August 2018


Verification of Convergence Rates of Numerical Solutions for Parabolic Equations
journal, June 2019

  • Jeong, Darae; Li, Yibao; Lee, Chaeyoung
  • Mathematical Problems in Engineering, Vol. 2019
  • DOI: 10.1155/2019/8152136

On the Thermodynamics of the Swift-Hohenberg Theory
text, January 2016