![]() |
Materials Algorithms Project
|
T.C. Illingworth and I.O. Golosnoy
Department of Materials Science and Metallurgy
University of Cambridge
Cambridge
CB2 3QZ
U.K.
E-mail: iog20@cam.ac.uk
Added to MAP: February 2005.
To track the position of the interface in a two phase diffusion-controlled moving boundary problem.
Language: | C/C++ |
Product form: | Source code. |
Even under isothermal conditions, the position of an interface between two different phases of a binary alloy can move if there are concentration gradients present in one (or both) phases. This code models diffusion in either phase and tracks thee way in which the interface moves for conditions where the interfacial concentrations are constant (i.e. time-invariant).
The algorithm can be used for either planar or spherical geometries and relies on an implicit front-fixing (IFF) algorithm. In contrast to some previous models, it has been derived in such a way as to ensure that it conserves solute.
None. Simply change the first few lines of the main routine in ConservativeIFF.
A file called 'results' is generated upon exectution of program; it describes how the position of the interface varies with time.
None.
Conserves solute in every case. Results indicate that the algorithm is first order accuracy in both space and time.
This code constitutes implements a fully implicit algorithm - there is concequently no limitation on the size of timestep that can be used to generate stable solutions.
Main routine is in ConservativeIFF.cpp. Other source files contain the subroutines that are called and header files.
Run as supplied.
No input parameters.
Generates a file (results.txt) which contains the following text:
Time | Interface Position |
0 | 1. |
0.1 | 1.00005 |
0.2 | 1.0001 |
0.3 | 1.00015 |
0.4 | 1.0002 |
0.5 | 1.00025 |
0.6 | 1.0003 |
0.7 | 1.00035 |
0.8 | 1.0004 |
0.9 | 1.00045 |
1 | 1.00049 |
None.
diffusion, phase-change, kinetics, moving boundary