[MAP Logo]

Materials Algorithms Project
Program Library


  1. Provenance of code.
  2. Purpose of code.
  3. Specification.
  4. Description of subroutine's operation.
  5. References.
  6. Parameter descriptions.
  7. Error indicators.
  8. Accuracy estimate.
  9. Any additional information.
  10. Example of code
  11. Auxiliary subroutines required.
  12. Keywords.
  13. Download source code.
  14. Links.

Provenance of Source Code

P. E. J. Rivera Díaz del Castillo
Phase Transformations and Complex Properties Group,
Department of Materials Science and Metallurgy,
University of Cambridge,
Cambridge, U.K.

Top | Next


To obtain the variation of dimensionless supersaturation with Péclet number and the ratio of the plate tip radius and the critical radius for nucleation.

Top | Next | Prev


This program reads the value of Q, and a table with the variation of M1 and M2 with Péclet number from the input data file 'indat11.d' (see MAP_KINETIC_PARAMETER_S1), and writes the output in the file 'res1.d'.
Product form:Source code

Top | Next | Prev


The program obtains the variation of dimesnonless supersaturation with Péclet number and the ratio of the plate tip radius and the critical radius for nucleation by simultaneously solving the diffusion equation expressed in parabolic coordinates and an equation to account for mass balance at the parabolic interface. The equations can be found in the references.

Top | Next | Prev


  1. P. E. J. Rivera Díaz del Castillo and H. K. D. H. Bhadeshia, 'Growth of needle and plate shaped particles: theory for small supersaturations, maximum velocity hypothesis', Materials Science and Technology, 17 25-29 2001.
  2. P. E. J. Rivera Díaz del Castillo, Kinetics of Precipitation Reactions, Ph.D. Thesis, University of Cambridge.
  3. R. Trivedi, 'The role of interfacial energy and interface kinetics during the growth of precipitate plates and needles', Metallurgical Transactions, 1 921-927 1970.
Top | Next | Prev


Input parameters

Q - double precision
parameter to characterise the relative importance of interface kinetics and diffusion effects.

PEC - double precision array of size 6053
Péclet number.

N1, N2 - double precision arrays of size 6053
M1, M2 functions for plate growth.

Output parameters

O - double precision
Dimensionless supersaturation.

RRR - double precision
Ratio between plate tip radius and critical radius for nucleation.

P - double precision
Péclet number.

Top | Next | Prev

Error Indicators


Top | Next | Prev


No information supplied.

Top | Next | Prev

Further Comments

The program is set to write the results to files 'res1.d', a sample version of that file is provided.

Top | Next | Prev


The given example is for Q=0.01.

1. Program data

See input data file 'indat11.d'.

2. Program results

See sample results file 'res1.d'.
Top | Next | Prev

Auxiliary Routines

The subroutines called by this program are :-


Top | Next | Prev


capillarity, interface kinetics, plate, growth, Trivedi

Top | Next | Prev


Download source code

Top | Prev

MAP originated from a joint project of the National Physical Laboratory and the University of Cambridge.
MAP Website administration / map@msm.cam.ac.uk

Top | Index | MAP Homepage Valid HTML 3.2!