Abstract
We develop an algorithm and computational implementation for simulation of problems that combine Cahn–Hilliard type diffusion with finite strain elasticity. We have in mind applications such as the electro-chemo-mechanics of lithium ion (Li-ion) batteries. We concentrate on basic computational aspects. A staggered algorithm is proposed for the coupled multi-field model. For the diffusion problem, the fourth order differential equation is replaced by a system of second order equations to deal with the issue of the regularity required for the approximation spaces. Low order finite elements are used for discretization in space of the involved fields (displacement, concentration, nonlocal concentration). Three (both 2D and 3D) extensively worked numerical examples show the capabilities of our approach for the representation of (i) phase separation, (ii) the effect of concentration in deformation and stress, (iii) the effect of strain in concentration, and (iv) lithiation. We analyze convergence with respect to spatial and time discretization and found that very good results are achievable using both a staggered scheme and approximated strain interpolation.
Similar content being viewed by others
Notes
Also called spinodal decomposition.
References
Areias P Simplas. http://www.simplas-software.com
Cahn JW, Hilliard JE (1958) Free energy of a nonuniform system. I. Interfacial free energy. J Chem Phys 28:258–267
Di Leo CV, Rejovitzky E, Anand L (2014) A Cahn–Hilliard-type phase-field theory for species diffusion coupled with large elastic deformations: application to phase-separating Li-ionelectrode materials. J Mech Phys Solids 70:1–29
Elliot CM, Songmu Z (1986) On the Cahn=-Hilliard equation. Arch Ration Mech Anal 96(4):339–357
Korelc J (2002) Multi-language and multi-environment generation of nonlinear finite element codes. Eng Comput 18(4):312–327
Lim GT, Reddy JN (2003) On canonical bending relationships for plates. Int J Solids Struct 40:3039–3067
McDowell MT, Lee SW, Nix WD, Cui Y (2013) 25th Anniversary article: understanding the lithiation of silicon and other alloying anodes for lithium-ion batteries. Adv Mater 25(36):4966–4985
Miehe C, Mauthe S, Ulmer H (2014) Formulation and numerical exploitation of mixed variational principles for coupled problems of Cahn–Hilliard-type and standard diffusion in elastic solids. Int J Numer Methods Eng 99:737–762
Peerlings RHJ, de Borst R, Brekelmans WAM, de Vree JHP (1996) Gradient enhanced damage for quasi-brittle materials. Int J Numer Methods Eng 39:3391–3403
Sethuraman VA, Srinivasan V, Bower AF, Guduru PR (2010) In situ measurements of stress–potential coupling in lithiated silicon. J Electrochem Soc 157:A1253–A1261
Truesdell C (1985) The elements of continuum mechanics. Springer, New York Corrected Second Printing
Ubachs RLJM, Schreurs PJG, Geers MGD (2004) A nonlocal diffuse interface model for microstructure evolution of tin-lead solder. J Mech Phys Solids 52:1763–1792
Van der Ven A, Bhattacharya J, Belak AA (2013) Understanding Li diffusion in Li-intercalation compounds. Acc Chem Res 46(5):1216–1225
Van der Vorst HA (1992) A fast and smoothly converging variant of Bi-CG for the solution of nonsymmetric linear systems. SIAM J Sci Stat Comput 13(2):631–644
Wells GN, Kuhl E, Garikipati K (2006) A discontinuous Galerkin method for the Cahn–Hilliard equation. J Comput Phys 218:860–877
Wolfram Research Inc. (2007) Mathematica. Wolfram Research, Champaign
Acknowledgments
We gratefully acknowledge the Financial support by means of ERC-CoG COMBAT (Computational Modeling and Design of Lithium-ion batteries). We also gratefully acknowledge the useful discussions with Dr. Jahed Naghipoor who pointed out the practical importance of swelling during lithiation.
Author information
Authors and Affiliations
Corresponding author
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
About this article
Cite this article
Areias, P., Samaniego, E. & Rabczuk, T. A staggered approach for the coupling of Cahn–Hilliard type diffusion and finite strain elasticity. Comput Mech 57, 339–351 (2016). https://doi.org/10.1007/s00466-015-1235-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00466-015-1235-1