MAP Logo

Materials Algorithms Project
Kinetic Theory: Program Library

[A logo showing the University of Cambridge Crest]


Program MAP_KINETIC_MOVING_BOUNDARY

Provenance of Source Code

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.

Top | Next

Purpose

To track the position of the interface in a two phase diffusion-controlled moving boundary problem.

Top | Next | Prev

Specification

Language: C/C++
Product form: Source code.

No specific input parameters (change the first few lines of the source code instead).

Top | Next | Prev

Description

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.

Top | Next | Prev

References

  1. T.C. Illingworth, I.O. Golosnoy, V. Gergely, T.W. Clyne Numerical Modelling of Transient Liquid Phase Bonding and Other Diffusion Controlled Phase Changes, J. Mat. Sci., in press.
  2. T.C. Illingworth, I.O. Golosnoy Numerical Solutions of Diffusion-Controlled Moving Boundary Problems which Conserve Solute J. Comp. Phys., in press.

Top | Next | Prev

Parameters

Input parameters

None. Simply change the first few lines of the main routine in ConservativeIFF.

Output parameters

A file called 'results' is generated upon exectution of program; it describes how the position of the interface varies with time.

Top | Next | Prev

Error Indicators

None.

Top | Next | Prev

Accuracy

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.

Top | Next | Prev

Further Comments

Main routine is in ConservativeIFF.cpp. Other source files contain the subroutines that are called and header files.

Top | Next | Prev

Example

1. Program text

Run as supplied.

2. Program data

No input parameters.

3. Program results

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

Top | Next | Prev

Auxiliary Routines

None.

Top | Next | Prev

Keywords

diffusion, phase-change, kinetics, moving boundary

Top | Next | Prev

Download

Download source code

Top | Prev