[MAP Logo]

Materials Algorithms Project
Program Library


Program MAP_KINETIC_PARAMETER_R1

  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

Purpose

To obtain the variation of functions N1, N2, R1, R2 with Péclet number. These values can be used to calculate the needle lengthening rate using Trivedi's kinetics.

Top | Next | Prev

Specification

This is a self-contained program.
Language:FORTRAN
Product form:Source code

Top | Next | Prev

Description

The program obtains the variation of N1, N2, R1, R2 functions when the Péclet number is varied. The expressions for those functions can be found in the references.

Top | Next | Prev

References

  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, 'Growth of dendritic needles from a supercooled melt', Acta Metallurgica, 18 287-296 1970.
Top | Next | Prev

Parameters

Input parameters

none.

Output parameters

N1, N2 - double precision
N1, N2 functions for needle growth.

R1, R2 - double precision
R1, R2 functions for needle growth.

P - double precision
Péclet number.

Top | Next | Prev

Error Indicators

None.

Top | Next | Prev

Accuracy

The accuracy of the results can be set by modifying the code, but very large precision may cause a non-convergent soluton.

Top | Next | Prev

Further Comments


The program is set to write the results to files 'res1.d', 'res2.d', a sample version those files is provided.
The variation of N1 and N2 with Péclet number is tabulated in `res1.d', and the variation of R1 and R2 with Péclet number is tabulated in `res2.d'.
The program is set to calculate the values of N1, N2, R1, R2 for values of Péclet number between 1.0 and 3.0.

Top | Next | Prev

Example

See files `res1.d' and `res2.d'.

Top | Next | Prev

Auxiliary Routines

The subroutines called by this program are :-

ERROR
GAMMA
LGAMA
CHGU
PSI

Top | Next | Prev

Keywords

capillarity, interface kinetics, needle, growth, Trivedi

Top | Next | Prev

Download

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!