doi:10.1016/j.jcp.2007.08.016
Copyright © 2007 Elsevier Inc. All rights reserved.
Multiscale modeling of alloy solidification using a database approach
aMaterials Process Design and Control Laboratory, Sibley School of Mechanical and Aerospace Engineering, 188 Frank H.T. Rhodes Hall, Cornell University, Ithaca, NY 14853-3801, USA
Received 20 March 2007;
revised 14 August 2007;
accepted 16 August 2007.
Available online 31 August 2007.
References and further reading may be available for this article. To view references and further reading you must
purchase this article.
Abstract
A two-scale model based on a database approach is presented to investigate alloy solidification. Appropriate assumptions are introduced to describe the behavior of macroscopic temperature, macroscopic concentration, liquid volume fraction and microstructure features. These assumptions lead to a macroscale model with two unknown functions: liquid volume fraction and microstructure features. These functions are computed using information from microscale solutions of selected problems. This work addresses the selection of sample problems relevant to the interested problem and the utilization of data from the microscale solution of the selected sample problems. A computationally efficient model, which is different from the microscale and macroscale models, is utilized to find relevant sample problems. In this work, the computationally efficient model is a sharp interface solidification model of a pure material. Similarities between the sample problems and the problem of interest are explored by assuming that the liquid volume fraction and microstructure features are functions of solution features extracted from the solution of the computationally efficient model. The solution features of the computationally efficient model are selected as the interface velocity and thermal gradient in the liquid at the time the sharp solid–liquid interface passes through. An analytical solution of the computationally efficient model is utilized to select sample problems relevant to solution features obtained at any location of the domain of the problem of interest. The microscale solution of selected sample problems is then utilized to evaluate the two unknown functions (liquid volume fraction and microstructure features) in the macroscale model. The temperature solution of the macroscale model is further used to improve the estimation of the liquid volume fraction and microstructure features. Interpolation is utilized in the feature space to greatly reduce the number of required sample problems. The efficiency of the proposed multiscale framework is demonstrated with numerical examples that consider a large number of crystals. A computationally intensive fully-resolved microscale analysis is also performed to evaluate the accuracy of the multiscale framework.
Keywords: Multiscale modeling; Solidification; Database approach; Level set method; Crystal growth
Fig. 1. Schematic of average volume to obtain microstructure features at two points. The features Λ are defined by statistical averaging of the results of appropriately defined microscale directional solidification problems.
Fig. 2. Schematic of solution features of model M (VM(x) and
) as applied to the domain of the problem of interest.
Fig. 3. (Left) Domain for computation (outer rectangle) and domain for performing averaging (inner rectangle). (Right) Schematic of the process for obtaining microstructure features.
Fig. 4. Schematic of the multiscale framework (steps indicated with the dark arrows).
Fig. 5. Schematic of the overall algorithm.
Fig. 6. Schematic of using interpolation to reduce the number of the required sample problem solutions. The top shows how we can use the time consuming fully-resolved model to evaluate f and Λ corresponding to an arbitrary feature. The bottom shows all obtained solution features (gray dots), a mesh generated in the feature space and how interpolation is used to obtain f and Λ based on results of sample problems corresponding to nodes of the mesh (black dots).
Fig. 8. Contour of VM (left two plots) and
(right two plots) for temperature boundary condition Tb = 50 exp(−t/10) − 40 (the first and the third plots) and Tb = 100 exp(−t/20) − 90 (the second and fourth plots).
Fig. 9. (Left) Obtained features of model M, (VM and
) with temperature boundary condition Tb = 50 exp(−t/10) − 40 (red square symbols) and Tb = 100 exp(−t/20) − 90 (green triangle symbols). (Right) (VM,
) of 64 sample problems selected for applying the fully-resolved model. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)
Fig. 10. Obtained microstructure of 64 sample runs. Each rectangle corresponds to a feature (VM,
) shown in the right plot of Fig. 9. The temperature gradient increases from left to right on a given row (constant growth velocity) and the growth velocity increases from bottom to top on a given column (constant temperature gradient).
Fig. 11. Obtained relation between the liquid volume fraction and temperature from sample problems with various features FM (64 sample problems are computed but only 9 are plotted here).
Fig. 12. Predicted field of
at the first, second, and third iteration for the case with Tb = 100 exp(−t/20) − 90. Contour line with value 0.7 demonstrates the predicted location of columnar to equiaxed transition (CET).
Fig. 13. Field of
(the 14 microstructures shown correspond to the closest microstructure in the database). (Left) Tb = 50 exp(−t/10) − 40. (Right) Tb = 100 exp(−t/20) − 90.
Fig. 14. Field of
(the 14 microstructures correspond to the closest microstructure in the database). (Left) Tb = 50 exp(−t/10) − 40. (Right) Tb = 100 exp(−t/20) − 90.
Fig. 15. Comparison of the predicted microstructures using the database approach with the microstructures obtained from solving the problem in the whole domain using the microscale model. (Left) Tb = 50 exp(−t/10) − 40. (Right) Tb = 100 exp(−t/20) − 90. For each plot (left or right): the picture in the middle is the fully-resolved result; the dark line in the middle picture is the predicted location of CET transition using the database approach
; the 14 pictures (around the middle picture) are the closest microstructure in the database based on features FM at selected locations.
Fig. 16. (Left) Predicted temperature field and liquid volume fraction contours with values 0.95 and 0.05 at time 130 for the case with Tb = 50 exp(−t/10) − 40. (Right) Obtained microstructure (using microscale model) and liquid volume fraction contours with value 0.95 and 0.05 (using the database approach) at time 130 for the case with Tb = 50 exp(−t/10) − 40.
Fig. 17. Microstructure at time 81.6 for the case with Tb = 100 exp(−t/20) − 90 using different sampling of potential nucleation sites.
Fig. 18. (Left 3 plots) Microscopic temperature at time 81.6 for the case with Tb = 100 exp(−t/20) − 90 using different sampling of potential nucleation sites. (Right plot) Averaged microscopic temperature.
Fig. 19. Temperature field at time 81.6 for the case with Tb = 100 exp(−t/20) − 90. Contour line shows the position where the temperature is −2. (Left) Temperature field obtained from the microscale model by averaging among three computation results. (Middle) Predicted temperature field from the database approach. (Right) Predicted temperature using the Lever rule.
Fig. 20. Schematic of the solidification of an Al–Cu alloy. The very small rectangle inside the elliptic shape is used to demonstrate the relevant size of the domain of the sample problem versus the size of the domain of the problem of interest. It is magnified by 50 times in the right plot.
Fig. 21. Obtained GM (left) and VM (right) fields using model M. Units of axes, GM and VM are mm, K/μm and μm/s, respectively.
Fig. 22. (Left) Obtained solute features in the VM and GM coordinates. (Right) VM and GM for the sample runs. Units of GM and VM are K/μm and μm/s, respectively.
Fig. 23. Demonstration of sample problem domain with periodic boundary conditions applied at the top and bottom sides. The bottom half is the computational domain, the top half is just a copy of the solution from the bottom.
Fig. 24. Sample problem results using the fully-resolved model.
Fig. 25. Relation of volume fraction and temperature.
Fig. 26. Temperature field at time 40 s. (Left) Predicted by Lever rule. (Right) Predicted by the database approach.
Fig. 27. GM and VM fields (left half is result of the first iteration, right half is result after three iterations) and nearest microstructure in database at location A (95 mm, 75 mm), B (90 mm, 75 mm), C (75 mm, 75 mm), D (60 mm, 80 mm), E (90 mm, 10 mm), F (80 mm, 20 mm), G (65 mm, 35 mm), H (50 mm, 50 mm).