TABULATE_SUBSTANCE
This command is for the TABULATE_REACTION (TAB) module.
Tabulate thermodynamic properties of a substance (with a given chemical formula but it is not certain in which phase/state it may form) or a pure stoichiometric phase (getting data from a specific substance database such as SSUB, or from a specific solution database such as SSOL and TCFE), or of a solution phase with a fixed composition (getting data from a specific solution database such as SSOL and TCFE), under a constant pressure and various temperatures.
In case of a pure stoichiometric phase or a solution phase with a fixed composition, you must already define the system and get thermodynamic data from an appropriate solution database in the TDB module before using this TAB command. The phase name can be given either upper or lower or mixed cases, and can be abbreviated until it is unique, e.g. fcc
, BCC
, Liq
, cem, Al5Fe2
, etc. as the SSOL database is used for the tabulations, or GAS
, FE-S
, wustite
, Fe2O3-hem
, etc. as the SSUB database is used. The module then optionally goes through a specific sublattice or all the possible sublattices for the specified phase, as well as through the whole framework of a currently defined system (including all the defined elements/species, as well as the default defined vacancy and electron if necessary for the current database), when it prompts for inputs of site fractions for various species possibly occupying such sublattices. Therefore, you can tabulate thermodynamic properties for a pure phase, or for a solution phase as a pure end-member (with a composition definition for the corresponding non-interacting species on all necessary sublattices in the end-member), or for a solution phase as a real solution (with a composition definition for related interacting species on all necessary sublattices in the solution).
Moreover, the variation for one of the tabulated properties for a substance, a pure stoichiometric phase, or a solution phase with a fixed composition can be plotted against temperature as a graph or saved as an experimental file (with an extension of .EXP).
Syntax |
tabulate_substance |
Prompt |
Substance (phase): <Name of the species or phase> In case of a substance, give its chemical formula, e.g. Fe, H2, C1H6, FeC1, CaC1O3, MgSiO3, etc. When the TABULATE_SUBSTANCE <substance> command sequence is used for the first time in the TAB module, it always uses the SSUB Substances database as the default. If SSUB is not available and a specific database has not been switched on yet, it may prompt to specify an appropriate substance or solution database (e.g. typing Before other prompts, a list of the used database, defined elements, getting-data sequence, references, etc. is displayed on screen. From such information, you know what the TAB module is performing. In case of a pure stoichiometric phase or a solution phase, give its phase name, e.g. FCC, CEMENTITE, LIQUID, SLAG, AQUEOUS, GAS, Al5Fe2, Fe-S, Wustite, Fe2O3-Hematite, etc. REMEMBER: It is important that you have already defined the system and gotten thermodynamic data from an appropriate solution database in the TDB module before using this TAB_SUB command. It automatically prompts for other options and necessary inputs for the definition of the composition of the specified phase. For a pure stoichiometric phase defined with one single sublattice (such as Fe-S, Wustite, and Fe2O3-Hematite), no further composition definition is needed. For a solution phase defined with one sublattice (such as AQUEOUS solution, GAS mixture and SLAG solution), it requests n-1 site-fraction inputs for the phase constituents (if there are n species in the defined phase within the whole framework of currently defined system, including all the defined elements/species, as well as the default defined vacancy and electron if necessary for the current database); the nth species is automatically assigned as the rest. The sum of input site fractions must not exceed unity. For example, the following prompts and inputs can be seen for a LIQUID solution phase (from the SSOL database) within the Fe-Cr-Ni-C-N-O framework (Note that by default, the unprompted Ni species is assigned as the rest): |
Fraction of constituent (return for prompt): <Return> C /1/: .05 CR /1/: .1 FE /1/: <RETURN> |
|
Sum of fractions exceed unity, please reenter FE /1/: .8 N /1/: .005 For a pure stoichiometric phase defined with two or more sublattices (such as Al5Fe2 inter-metallic stoichiometric phase) and for a solution phase defined with two or more sublattices (such as FCC alloy solution, and ION_LIQ ionic liquid solution, Al5Fe4 inter-metallic solution phase), it first asks for which sublattice the site-fractions of constituents shall be specified: the default value 0 for all possible sublattices, and a given positive number for a specific sublattice (this number must be reasonable for the currently specified phase, i.e. it must be smaller than the total sublattice number of the phase). Then, for the given sublattice or for all sublattices, it prompts for all the necessary inputs of site-fractions for the possible constituents on each sublattice n-1 times if there are n species on such a sublattice in the defined phase; the nth species on each of the sublattices are automatically assigned as the rest. For example, these prompts and inputs can be seen for an FCC solution phase (from the SSOL database) within the Fe-Cr-Ni-C-N-O framework (Note that by default, the unprompted |
|
Specify sublattice (0 for all) /0/: FRACTIONS IN SUBLATTICE 1 CR /1/: .1 FE /1/: .8 NI /1/: .0995 FRACTIONS IN SUBLATTICE 2 C /1/: .05 N /1/: .05 |
|
Pressure /100000/: <Pressure of interest, in pa> Specify the constant pressure condition in Pa. |
|
Low temperature limit /298.15/: <T-low, in K> Specify the starting temperature K. |
|
High temperature limit /298.15/: <T-high, in K> Specify the ending temperature K. |
|
Step in temperature /100/: <T-step> Specify the temperature step for the tabulation. |
|
Output file /Screen/: <Return for screen, or type a file name> If you type If typing a file name, the table is both displayed on screen and saved as an *.EXP or a *.TAB file, and the program further prompts: |
|
Graphical output /Y/:
|
|
Plot column ? /2/: <1 or 2 or 3 or 4 or 5 or 6> Specify which property (as column number) to be plotted as the Y-axis (versus temperature as X-axis) as a graph on the screen. Simultaneously, all tabulated properties and the Y-axis setting (i.e. the plotted column) for the graph is written into an *.EXP file using the DATAPLOT format. The default column 2 is heat capacity, 3 enthalpy, 4 entropy and 5 Gibbs energy, and the additional column 6 is a user-entered function. The table displays on the screen. The plot then displays on the screen, followed by a POST: prompt. The POST module (postprocessor) automatically opens, and all types of the POST-module commands to refine the plotted diagram may be used. Possibilities include scaling the X/Y-axis, changing the X/Y-axis texts, etc. The command BACK or EXIT at the POST: prompt always takes you back to the TAB module. |

For the FCC
alloy phase as a non-ideal solution with a fixed composition
[Fe0.80,Cr0.10,Ni0.0995,O0.005][C0.05,N0.05,VA0.90]
The following table is obtained by typing Return for screen
at the prompt Output file /SCREEN/
:
O U T P U T F R O M T H E R M O - C A L C
Phase : FCC Pressure : 100000.00
Specie: CR1/--2
***************************************************************
T Cp H S G
(K) (Joule/K) (Joule) (Joule/K) (Joule)
***************************************************************
298.15 2.70517E+01 6.23824E+03 4.40241E+01 -6.88755E+03
300.00 2.70889E+01 6.28832E+03 4.41916E+01 -6.96915E+03
400.00 2.87304E+01 9.08420E+03 5.22235E+01 -1.18052E+04
500.00 2.99904E+01 1.20222E+04 5.87742E+01 -1.73649E+04
600.00 3.10889E+01 1.50770E+04 6.43408E+01 -2.35275E+04
700.00 3.21116E+01 1.82375E+04 6.92106E+01 -3.02100E+04
800.00 3.30994E+01 2.14982E+04 7.35633E+01 -3.73524E+04
900.00 3.40742E+01 2.48569E+04 7.75182E+01 -4.49094E+04
1000.00 3.50483E+01 2.83130E+04 8.11586E+01 -5.28456E+04
1100.00 3.60268E+01 3.18667E+04 8.45449E+01 -6.11327E+04
1200.00 3.70143E+01 3.55187E+04 8.77219E+01 -6.97476E+04
1300.00 3.80149E+01 3.92700E+04 9.07241E+01 -7.86713E+04
1400.00 3.90311E+01 4.31222E+04 9.35784E+01 -8.78875E+04
1500.00 4.00649E+01 4.70768E+04 9.63064E+01 -9.73827E+04
1600.00 4.11174E+01 5.11358E+04 9.89256E+01 -1.07145E+05
1700.00 4.21896E+01 5.53010E+04 1.01450E+02 -1.17165E+05
1800.00 4.33555E+01 5.95774E+04 1.03894E+02 -1.27433E+05
1900.00 4.58528E+01 6.40379E+04 1.06306E+02 -1.37943E+05
2000.00 4.75402E+01 6.87138E+04 1.08704E+02 -1.48694E+05

For a phase, as either an end-member or real solution, the listed species name is irrelevant.
For the H2 species as a pure substance (using the SSUB5 database), the following table is obtained by typing Return for screen
at the prompt Output file /SCREEN/
, whilst the figure is generated by typing 5 at the prompt Plot column /2/:
O U T P U T F R O M T H E R M O - C A L C
Phase : GAS Pressure : 100000.00
Specie: H2
*************************************************************** T Cp H S G
(K) (Joule/K) (Joule) (Joule/K) (Joule)
*************************************************************** 298.15 2.88369E+01 3.17684E-06 1.30680E+02 -3.89622E+04
300.00 2.88473E+01 5.33580E+01 1.30858E+02 -3.92042E+04
400.00 2.91591E+01 2.95686E+03 1.39209E+02 -5.27268E+04
500.00 2.92650E+01 5.87874E+03 1.45729E+02 -6.69856E+04
600.00 2.93441E+01 8.80908E+03 1.51071E+02 -8.18336E+04
700.00 2.94579E+01 1.17488E+04 1.55602E+02 -9.71730E+04
800.00 2.96320E+01 1.47027E+04 1.59547E+02 -1.12935E+05
900.00 2.98786E+01 1.76776E+04 1.63050E+02 -1.29068E+05
1000.00 3.02043E+01 2.06810E+04 1.66214E+02 -1.45533E+05
1100.00 3.05319E+01 2.37171E+04 1.69108E+02 -1.62302E+05
1200.00 3.09281E+01 2.67897E+04 1.71781E+02 -1.79348E+05
1300.00 3.13615E+01 2.99040E+04 1.74274E+02 -1.96652E+05
1400.00 3.18115E+01 3.30625E+04 1.76614E+02 -2.14197E+05
1500.00 3.22641E+01 3.62664E+04 1.78824E+02 -2.31970E+05
1600.00 3.27094E+01 3.95151E+04 1.80921E+02 -2.49958E+05
1700.00 3.31406E+01 4.28078E+04 1.82917E+02 -2.68151E+05
1800.00 3.35526E+01 4.61426E+04 1.84823E+02 -2.86539E+05
1900.00 3.39415E+01 4.95175E+04 1.86648E+02 -3.05113E+05
2000.00 3.43045E+01 5.29300E+04 1.88398E+02 -3.23866E+05