Hostname: page-component-7479d7b7d-k7p5g Total loading time: 0 Render date: 2024-07-11T17:38:20.689Z Has data issue: false hasContentIssue false
  • English
  • Français

An Efficient Rescaling Algorithm for Simulating the Evolution of Multiple Elastically Stressed Precipitates

Published online by Cambridge University Press:  03 June 2015

Amlan K. Barua ,
Shuwang Li ,
Hualong Feng ,
Xiaofan Li  and
John Lowengrub
Amlan K. Barua*
Affiliation:
Department of Applied Mathematics, Illinois Institute of Technology, Chicago 60616, USA
Shuwang Li*
Affiliation:
Department of Applied Mathematics, Illinois Institute of Technology, Chicago 60616, USA
Hualong Feng*
Affiliation:
Department of Applied Mathematics, Illinois Institute of Technology, Chicago 60616, USA
Xiaofan Li*
Affiliation:
Department of Applied Mathematics, Illinois Institute of Technology, Chicago 60616, USA
John Lowengrub*
Affiliation:
Department of Mathematics, University of California at Irvine, Irvine, CA 92697, USA
*
Email:abarua@iit.edu
Corresponding author.Email:sli@math.iit.edu
Email:hfeng@iit.edu
Email:lix@iit.edu
Email:lowengrb@math.uci.edu
Get access

Abstract

In this paper, we propose a space-time rescaling scheme for computing the long time evolution of multiple precipitates in an elastically stressed medium. The algorithm is second order accurate in time, spectrally accurate in space and enables one to simulate the evolution of precipitates in a fraction of the time normally used by fixed-frame algorithms. In particular, we extend the algorithm recently developed for single particle by Li et al. (Li, Lowengrub and Leo, J. Comput. Phys., 335 (2007), 554) to the multiple particle case, which involves key differences in the method. Our results show that without elasticity there are successive tip splitting phenomena accompanied by the formation of narrow channels between the precipitates. In presence of applied elastic field, the precipitates form dendrite-like structures with the primary arms aligned in the principal directions of the elastic field. We demonstrate that when the far-field flux decreases with the effective radius of the system, tip-splitting and dendrite formation can be suppressed, as in the one particle case. Depending on the initial position of the precipitates, we further observe that some precipitates grow while others may shrink, even when a positive far field flux is applied.

Keywords

74N2574B15Diffusionelasticitydendriteside branchingboundary integral method
Type
Research Article
Information
Communications in Computational Physics , Volume 14 , Issue 4 , October 2013 , pp. 940 - 959
Copyright
Copyright © Global Science Press Limited 2013

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Akaiwa, N., Thornton, K. and Voorhees, P. W., Large-scale simulations of microstructural evolution of elastically stressed solids, J. Comput. Phys., 173 (2004), 61.CrossRefGoogle Scholar
[2]Ardell, A. J. and Nicholson, R. B., On modulated structure of aged ni-al alloys, Acta. Metall., 14 (1966), 1225.Google Scholar
[3]Barua, A., Li, S., Li, X. and Lowengrub, J., Self similar evolution of a precipitate in inhomogeneous elastic media, J. Cryst. Growth, 351 (2012), 62.Google Scholar
[4]Bates, P. W., Chen, X. and Deng, X., A numerical scheme for the two phase mullins-sekerka problem, Electron. J. Differential Equations, 11 (1995), 27.Google Scholar
[5]Belytschko, T., Gracie, R. and Ventura, G., A review of extended/generalized finite element methods for material modeling, Modell. Simul. Mater. Sci. Eng., 17(4) (2009).Google Scholar
[6]Boisse, J., Lecoq, N., Patte, R. and Zapolsky, H., Phase-field simulation of coarsening of gamma precipitates in an ordered gamma matrix, Acta Mater., 55(18) (2007), 6151–6158.CrossRefGoogle Scholar
[7]Brebbia, C., Telles, J. and Wrobel, L., Boundary Element Techniques, Springer-Verlag, Berlin, 1984.Google Scholar
[8]Cheng, H., Greengard, L. and Rokhlin, V., A fast adaptive multipole algorithm in three dimensions, J. Comput. Phys., 155(2) (1999), 468–498.Google Scholar
[9]Cristini, V. and Lowengrub, J. S., Three-dimensional crystal growth: nonlinear simulation and control of the mullins sekerka instability, J. Cryst. Growth, 266 (2004), 552.Google Scholar
[10]Feng, H., Barua, A., Li, S. and Li, X., A fast tree code algorithm for elastic kernel functions, in review, 2012.Google Scholar
[11]Greenbaum, A., Greengard, L. and McFadden, G. B., Laplace equation and dirichlet-neuman map in multiply connected domains, J. Comput. Phys., 105 (1993), 267.Google Scholar
[12]Greengard, L. and Helsing, J., On the numerical evaluation of elastostatic fields in locally isotropic two-dimensional composites, J. Mech. Phys. Solids, 46(8) (1998), 1441–1462.Google Scholar
[13]Greengard, L. and Rokhlin, V., A Fast algorithm for particle simulations, J. Comput. Phys., 73(2) (1987), 325–348.Google Scholar
[14]Hou, T., Lowengrub, J. S. and Shelley, M., Removing the stiffness from interfacial flows with surface tension, J. Comput. Phys., 114 (1994), 312.CrossRefGoogle Scholar
[15]Jou, H.-J., Leo, P. H. and Lowengrub, J. S., Microsructural evolution in inhomogenous elastic media, J. Comput. Phys., 131 (1997), 109.Google Scholar
[16]Koyama, T., Simulation of microstructural evolution based on the phase-field method and its applications to material development, J. Jpn. Inst. Met., 73(12) (2009), 891–905.CrossRefGoogle Scholar
[17]Kwon, Y., Thornton, K. and Voorhees, P. W., The topology and morphology of bicontinuous interfaces during coarsening, EPL, 86(4) (2009).Google Scholar
[18]Kwon, Y., Thornton, K. and Voorhees, P. W., Morphology and topology in coarsening of domains via non-conserved and conserved dynamics, Philos. Mag., 90(1-4) (2010), 317–335.CrossRefGoogle Scholar
[19]Lashuk, I., Chandramowlishwaran, A., Langston, H., Nguyen, T.-A., Sampath, R., Shringarpure, A., Vuduc, R., Ying, L., Zorin, D. and Biros, G., A massively parallel adaptive fast multipole method on heterogeneous architectures, Commun. ACM, 55(5) (2012), 101–109.Google Scholar
[20]Leo, P. H. and Jou, Herng Jeng, Shape evolution of an initially circular precipitate growing by diffusion in an applied stress field, Acta Metall., 41 (1993), 2271.Google Scholar
[21]Leo, P. H., Lowengrub, J. S. and Nie, Q., Microstructural evolution in orthotropic elastic media, J. Comput. Phys., 157 (2000), 44.Google Scholar
[22]Leo, P. H. and Sekerka, R. F., The effect of elastic fields on the morphological stability of a precipitate grown from solid solution, Acta Metall., 37 (1989), 3139.Google Scholar
[23]Li, B., Lowengrub, J. S., Raetz, A. and Voigt, A., Geometric evolution laws for thin crystalline films: modeling and numerics, Commun. Comput. Phys., 6(3) (2009), 433–482.Google Scholar
[24]Li, S. and Li, X., A boundary integral method for computing the dynamics of an epitaxial island, SIAM J. Sci. Comput., 33 (2011), 3282.Google Scholar
[25]Li, S., Lowengrub, J., Fontana, J. and Palffy-Muhoray, P., Control of viscous fingering patterns in a radial hele-shaw cell, Phys. Rev. Lett., 102 (2009), 174501.Google Scholar
[26]Li, S., Lowengrub, J. and Leo, P., A rescaling scheme with application to the long time simulation of viscous fingering in a hele-shaw cell. J. Comput. Phys., 225 (2007), 554.Google Scholar
[27]Li, X., Thornton, K., Nie, Q., Voorhees, P. W. and Lowengrub, J. S., Two- and three-dimensional equilibrium morphology of a misfitting particle and the gibbs-thomson effect, Acta Mater., 52(20) (2004), 5829–5843.Google Scholar
[28]Lindsay, K. and Krasny, R., A particle method and adaptive treecode for vortex sheet motion in three-dimensional flow, J. Comput. Phys., 172(2) (2001), 879–907.CrossRefGoogle Scholar
[29]Liu, Y. J., A fast multipole boundary element method for 2D multi-domain elastostatic problems based on a dual BIE formulation, Comput. Mech., 42(5) (2008), 761–773.Google Scholar
[30]Maheshwari, A. and Ardell, A. J., Morphological evolution of coherent misfitting precipitates in anisotropic elastic media, Phys. Rev. Lett., 70 (1993), 2305.Google Scholar
[31]Mikhlin, S. G., Integral Equation and Their Applications to Certian Problems of Mechanics, Pergammon, New York, 1957.Google Scholar
[32]Mullins, W. and Sekerka, R. F., Morphological stability of a particle growing by diffusion of heat flow, J. Appl. Phys., 34 (1963), 323.CrossRefGoogle Scholar
[33]Perez, D. and Lewis, L. J., Multiscale model for microstructure evolution in multiphase materials: application to the growth of isolated inclusions in presence of elasticity, Phys. Rev. E, 74(3) (2006).Google Scholar
[34]Qin, R. S. and Bhadeshia, H. K., Phase field method, Mat. Sci. Technol., 26(7) (2010), 803–811.CrossRefGoogle Scholar
[35]Sampath, R. S. and Biros, G., A parallel geometric multigrid methos for finite elements on octree meshes, SIAM J. Sci. Comput., 32(3) (2010), 1361–1392.Google Scholar
[36]Sohn, J.-S., Tseng, Y.-H., Li, S., Lowengrub, J. and Voigt, A., Dynamics of multicomponent vesicles in a viscous fluid, J. Comput. Phys., 229 (2010), 119.Google Scholar
[37]Steinbach, I., Phase-field models in materials science, Modell. Simul. Mater. Sci. Eng., 17(7) (2009).Google Scholar
[38]Su, C. H. and Voorhees, P. W., Dynamics of precipitate evolution in elastically stressed solids. i. inverse coarsening, Acta Mater., 44 (1996), 1987.Google Scholar
[39]Su, C. H. and Voorhees, P. W., Dynamics of precipitate evolution in elastically stressed solids. ii. particle alignment, Acta Mater., 44 (1996), 2001.Google Scholar
[40]Thompson, M. E., Su, C. S. and Voorhees, P. W., Equilibrium shape of a misfitting precipitate, Acta Metall., 42 (1994), 2107.Google Scholar
[41]Thornton, K., Akaiwa, N. and Voorhees, P. W., Large-scale simulations of ostwald ripening in elastically stressed solids. coarsening kinetics and particle size distribution, Acta Mater., 52 (2004), 1365.Google Scholar
[42]Thornton, K., Akaiwa, N. and Voorhees, P. W., Large-scale simulations of ostwald ripening in elastically stressed solids. developement of microstructure, Acta Mater., 52 (2004), 1353.Google Scholar
[43]Thornton, K., Argen, J. and Voorhees, P. W., Modelling the evolution of phase boundaries in solids at the meso- and nano-scales, Acta Mater., 51 (2003), 5675.Google Scholar
[44]Voorhees, P. W., Boisvert, R. F. and Johnson, W. C., On the morphological development of second phase particles in elastically stressed solids, Acta Metall., 40 (1992), 2979.Google Scholar
[45]Voorhees, P. W., Boisvert, R. F., McFadden, G. B. and Meiron, D. I., Numerical simulation of morphological development during ostwald ripening, Acta Metall., 36 (1988), 207.Google Scholar
[46]Wada, R., Kobayashi, N., Fujii, T., Onaka, S. and Kato, M., Effects of elastic interaction energy and interface energy on morphology and distribution of precipitates: two-dimensional analysis based on anisotropic elasticity, J. Jpn. Inst. Met., 71(8) (2007), 587–591.CrossRefGoogle Scholar
[47]Wang, Q. X, Variable order revised binary treecode, J. Comput. Phys., 200(1) (2004), 192–210.Google Scholar
[48]Wang, Y. and Li, J., Phase field modeling of defects and deformation, Acta Mater., 58(4) (2010), 1212–1235.Google Scholar
[49]Ying, L., Biros, G. and Zorin, D., A kernel-independent adaptive fast multipole method in two and three dimensions, J. Comput. Phys., 196 (2004), 591.Google Scholar