Materials Algorithms Project
Program Library
- Provenance of code.
- Purpose of code.
- Specification.
- Description of subroutine's operation.
- References.
- Parameter descriptions.
- Error indicators.
- Accuracy estimate.
- Any additional information.
- Example of code
- Auxiliary subroutines required.
- Keywords.
- Download source code.
- Links.
H.K.D.H. Bhadeshia,
Phase Transformations Group,
Department of Materials Science and Metallurgy,
University of Cambridge,
Cambridge, U.K.
Top |
Next
To calculate the carbon concentration (mole fraction) at the T0 phase
boundary at a specified temperature.
Top |
Next |
Prev
Language: | FORTRAN
|
Product form: | Source code |
SUBROUTINE MAP_STEEL_AXTO(H,S,H1,S1,XTO,T,W,W1,F,
& AJ,AJ1,STORE,JN,J,R)
DOUBLE PRECISION H,S,H1,S1,XTO,T,W,W1,F,AJ,AJ1,
& STORE,R
DOUBLE PRECISION MAP_STEEL_FTO1, MAP_STEEL_G91
INTEGER J, JN
Top |
Next |
Prev
MAP_STEEL_AXTO calculates the carbon concentration (mole fraction) at the T0 phase boundary at a specified temperature. Newton's iterative method is used to refine an initial guess.
Top |
Next |
Prev
- H.K.D.H. Bhadeshia and D.V. Edmonds, Acta Metallurgica,
28, (1980), 1265-1273.
Top |
Next |
Prev
Input parameters
- H - real
- H is the excess partial molar enthalpy of solution of carbon in austenite
(in joules per mole, Jmol-1).
- S - real
- S is the excess partial molar entropy of solution of carbon in austenite
(in joules per mole, Jmol-1).
- H1 - real
- H1 is the excess partial molar enthalpy of solution of carbon in ferrite
(in joules per mole, Jmol-1).
- S1 - real
- S1 is the excess partial molar entropy of solution of carbon in ferrite (in
joules per mole, Jmol-1).
- T - real
- T is the temperature (in kelvin).
- W - real
- W is the carbon-carbon interaction energy in austenite (in joules per mole,
Jmol-1).
- W1 - real
- W1 is the carbon-carbon interaction energy in ferrite (in joules per
mole, Jmol-1).
- F - real
- F is the energy (in joules).
- AJ - real
- AJ = 1-e{W/RT}, where R is the universal gas constant.
- AJ1 - real
- AJ1 = 1-e{W1/RT}, where R is the universal gas
constant.
- STORE - real
- STORE is the stored energy (in joules per mole,
Jmol-1).
- JN - integer
- JN is the maximum number of iterations.
- R - real
- R is the universal gas constant (in joules per mole per kelvin,
Jmol-1K-1).
Output parameters
- XTO - real
- XTO is the carbon concentration at the T0 boundary.
- J - integer
- J is the number of iterations performed.
Top |
Next |
Prev
None.
Top |
Next |
Prev
No information supplied.
Top |
Next |
Prev
See also subroutine MAP_STEEL_OMEGA and function MAP_STEEL_ENERGY.
Top |
Next |
Prev
1. Program text
DOUBLE PRECISION H,S,H1,S1,XTO,T,W,W1,F,AJ,AJ1,STORE,R
INTEGER J,JN
INCLUDE 'map_constants_gas.f'
READ (5,*) H,S,H1,S1,T,W,W1
READ (5,*) F,AJ,AJ1,STORE
READ (5,*) JN
CALL MAP_STEEL_AXTO(H,S,H1,S1,XTO,T,W,W1,F,AJ,AJ1,STORE,
& JN,J,R)
WRITE (6,10) XTO
WRITE (6,20) J
10 FORMAT ('Carbon concentration at T0 boundary is ',F8.5)
20 FORMAT (I4,' iterations were required for convergence')
STOP
END
2. Program data
None supplied.
3. Program results
None supplied.
Top |
Next |
Prev
MAP_STEEL_FTO1
MAP_STEEL_G91
Top |
Next |
Prev
carbon, T0 boundary
Top |
Next |
Prev
Download source code
Top |
Prev