Simulation of cellular-dendritic solidification structures of binary alloys in three-dimensional growth using a multiparticle diffusion-limited aggregation model
Abstract
Solidification in binary alloy melts was simulated with two different three-dimensional models, that differ in the way the role of the solid/liquid interface energy on the solidification is expressed: either surface rearrangement or a surface curvature-dependent attachment probability was applied. In the melt, the simultaneous diffusion of all diffusing particles (atoms) was taken into account. Cellular –dendritic growth modes, as observed in practice, could be well simulated. The surface rearrangement model is essential for exhibiting details (as facetting) and the effect of the next-nearest neighbour interaction.
References
[1] W.W. Mullins, R.F. Sekerka: J. Appl. Phys. 35 (1967) 444.10.1063/1.1713333Search in Google Scholar
[2] J. Langer: Rev. Mod. Phys. 52 (1980) 1.10.1103/RevModPhys.52.1Search in Google Scholar
[3] W. Kurz, D.J. Fisher: Acta Metall. 29 (1981) 11.10.1016/0001-6160(81)90082-1Search in Google Scholar
[4] R. Trivedi, W. Kurz: Acta Metall. Mater. 42 (1994) 15.10.1016/0956-7151(94)90044-2Search in Google Scholar
[5] H. Müller –Krumbhaar, J.S. Langer: Acta Metall. 29 (1981) 145.10.1016/0001-6160(81)90095-XSearch in Google Scholar
[6] D.E. Porter, K.E. Easterling: Phase Transformations in Metals and Alloys, Chapmann and Hall, London (1981).Search in Google Scholar
[7] J. Crank: Mathematics of Diffusion, Oxford University Press (1957).Search in Google Scholar
[8] L.E. Reichl: A Modern Course in Statistical Physics, Edward Arnold (1980).Search in Google Scholar
[9] T.A. Witten, L.M. Sander: Phys. Rev. Lett. 47 (1981) 1400.10.1103/PhysRevLett.47.1400Search in Google Scholar
[10] T.A. Witten, L.M. Sander: Phys. Rev. B 27 (1983) 5656.10.1103/PhysRevB.27.5686Search in Google Scholar
[11] T.A.M. Haemers, D.G. Rickerby, E.J. Mittemeyer: Modelling Simul. Mater. Sci. Eng. 7 (1999) 23310.1088/0965-0393/7/2/007Search in Google Scholar
[12] R.F. Voss: Phys. Rev. B 30 (1984) 334.10.1103/PhysRevB.30.334Search in Google Scholar
[13] R.F. Voss: J. Stat. Phys. 36 (1984) 861.10.1007/BF01012945Search in Google Scholar
[14] A. Das, E.J. Mittemeijer: Phil. Mag A 81 (2001) 2725.10.1080/01418610108216666Search in Google Scholar
[15] A. Das, E.J. Mittemeijer: Z. f. Metallkunde 93 (2002) 459.10.3139/146.020459Search in Google Scholar
[16] T. Vicsek: Phys. Rev. Lett. 53 (1984) 2281.10.1103/PhysRevLett.53.2281Search in Google Scholar
[17] T. Vicsek: Phys. Rev. A 32 (1985) 3048.10.1103/PhysRevA.32.3084Search in Google Scholar
[18] R. Tao, M.A. Novotny, K. Kaski: Phys. Rev. A 38 (1986) 1018.Search in Google Scholar
[19] Y. Saito, T. Ueta: Phys. Rev. A 40 (1989) 3408.10.1103/PhysRevA.40.3408Search in Google Scholar
[20] A. Das, E.J. Mittemeijer: Met. & Mat. Trans. A 31 (2000) 2049.10.1007/s11661-000-0232-xSearch in Google Scholar
[21] A. Karma, W.J. Rapell: Phys. Rev. Lett. 77 (1996) 4050.10.1103/PhysRevLett.77.4050Search in Google Scholar PubMed
Appendix: Discretization of the Laplacian in an f. c. c. lattice
In one dimension, the discretization of the second derivative of ρi,t is (in the following, the subscripts i and t have been omitted.
where Δx is the distance between two lattice points.
In two dimensions this becomes
and for three dimensions, it follows
In order to arrive at an expression for Δ2ρx,y,z for the f. c. c. lattice where the individual ρ values in the discretized expression (as Eq. (A-1)) pertain to ρ values at nearest neighbour lattice sites, we proceed as follows.
First, recognizing that p nearest neighbours in the f. c. c. lattice are on a distance
where positions with coordinates (x, y, z) do not correspond with sites of the f. c. c. lattice. To get (partial) coincidence with sites of the f. c. c. lattice sites, Eq. (A-2) can be rotated over each of the three axes. Rotation of 45° over the x axis results in
with only the terms (ρ values) on the x axis not on a site of the f. c. c. lattice.
Rotation of 45° over the z axis results in
with only the terms (ρ values) on the z axis not on a site of the f. c. c. lattice.
Rotation of 45° over the y axis results in
with only the terms (ρ values) on the y axis not on a site of the f. c. c. lattice.
Now, to finally get an Equation with all terms (ρ values) pertaining to nearest neighbour sites of the f. c. c. lattice, first, Eqs. (A-3), (A-4) and (A-5) are added, secondly, Eq. (A-2) is subtracted from this sum. Then, after dividing the left- and right-hand sides of the Equation by 2, the final result is
which for an f. c. c. lattice represents the discretized Laplacian of ρi,tin terms of nearest neighbour values of ρi,t
© 2004 Carl Hanser Verlag, München
Articles in the same Issue
- Frontmatter
- Editorial
- Editorial
- Werner-Köster-Preis
- den Werner-Köster-Preis 2003
- Festabend zum 85. Geburtstag der DGM
- Festabend zum 85. Geburtstag der DGM
- Articles Applied
- NiAl – material for the application in gas turbines
- High-temperature creep behavior and hot-pressing consolidation of NiAl
- Direct Strip Casting of Magnesium
- Thermo-mechanical analysis of cast/mould interaction in casting processes
- Aluminum–Lithium alloy development for thixoforming
- On the influence of heavy warm reduction on the microstructure and mechanical properties of a medium-carbon ferritic –pearlitic steel
- Thermodynamic assessment of the Dy–Mg system
- Influence of pre-straining on mechanical properties of HSLA steel by using ball indentation technique
- Fabrication and application of co-continuous Al2O3/Al composite by reactive infiltration of molten Al into the preform of SiO2 with cordierite addition
- Articles Basic
- Simulation of cellular-dendritic solidification structures of binary alloys in three-dimensional growth using a multiparticle diffusion-limited aggregation model
- Microstructure and age-hardening effects of aluminium alloys with additions of scandium and zirconium
- Notifications / Mitteilungen
- Personal/Personelles
- Materials Week
- News/Aktuelles
- Books/Bücher
- Conferences / Konferenzen
Articles in the same Issue
- Frontmatter
- Editorial
- Editorial
- Werner-Köster-Preis
- den Werner-Köster-Preis 2003
- Festabend zum 85. Geburtstag der DGM
- Festabend zum 85. Geburtstag der DGM
- Articles Applied
- NiAl – material for the application in gas turbines
- High-temperature creep behavior and hot-pressing consolidation of NiAl
- Direct Strip Casting of Magnesium
- Thermo-mechanical analysis of cast/mould interaction in casting processes
- Aluminum–Lithium alloy development for thixoforming
- On the influence of heavy warm reduction on the microstructure and mechanical properties of a medium-carbon ferritic –pearlitic steel
- Thermodynamic assessment of the Dy–Mg system
- Influence of pre-straining on mechanical properties of HSLA steel by using ball indentation technique
- Fabrication and application of co-continuous Al2O3/Al composite by reactive infiltration of molten Al into the preform of SiO2 with cordierite addition
- Articles Basic
- Simulation of cellular-dendritic solidification structures of binary alloys in three-dimensional growth using a multiparticle diffusion-limited aggregation model
- Microstructure and age-hardening effects of aluminium alloys with additions of scandium and zirconium
- Notifications / Mitteilungen
- Personal/Personelles
- Materials Week
- News/Aktuelles
- Books/Bücher
- Conferences / Konferenzen