ADD_INITIAL_EQUILIBRIUM
This command is for the POLY module.
Add initial equilibrium points from which a phase diagram is calculated (through the MAP command).
For each initial equilibrium the first stage of the mapping varies the value of the given axis variable condition, changing it in the given direction. When a phase change is found it is saved as a start equilibrium, and if the symbol ">" was used, it continues the search until the end of the axis variable.
If no axis variable conditions differ between initial equilibria then the map results in several overlaid diagrams, in accordance with the selected conditions and the values.
The ADD_INITIAL_EQUILIBRIUM default command implicitly makes the MAP command search for phase changes according to settings set by ADVANCED_OPTIONS STEP_AND_MAP command. Default for STEP_AND_MAP are three initial equilibrium equidistantly along each axis variable, along the edge of the axis variable rectangle.
By default these default start points are used in MAP. If any ADD_INITIAL_EQUILIBRIUM is used these default start points are not used. They are again used if the command ADD default is used together with another ADD_INITIAL_EQUILIBRIUM command.
The ADD command with the default direction scans along the axis variables and generates start points each time the scan procedure crosses a phase boundary. In addition, it generates four start points, scanning cross the middle of each axis, if there is any solubility line that does not reach the axes. At the MAP command, a search for lines in the diagram is made along each direction of the axis variables in the diagram.
In this way, it should guarantee that all possible phase boundary lines in a phase diagram are found. It may take a little longer to execute than using the minimum number of start points, as some lines may be calculated more than once. But the POLY module remembers all node points and subsequently stops calculations along a line when it finds a known node point.
It is also possible to create a sequence of start points from one initial equilibria by appending a > after the direction at the ADD command. For example:
Direction /Default/: 2>
Direction /Default/: -2>
This generates one start point for each set of phase change in the positive direction of the axis 2 (or negative direction of the axis 2); this ensures finding all possible phase boundary lines (not just the first one) along such an axis direction.
This is particularly useful when you have a phase diagram with several lines with no intersection. It is thus possible to calculate e.g. an isopleth for a much more limited composition range. It is also useful for calculating CVD diagrams and Pourbaix diagrams.
Syntax |
ADD_INITIAL_EQUILIBRIUM |
Prompt |
Direction /Default/: <Direction code> Direction code(s):
The direction is important when the initial equilibrium point is in a single-phase region or when the diagram is an isopleth (tie-lines not in the plane of calculation). In such cases, the program will search for a line in the diagram (i.e., a line where the amount of a phase is zero) in the given direction. |
Normally, this command is not needed to calculate a property diagram (with the STEP_WITH_OPTIONS → NORMAL
command-sequence), as the STEP procedure starts from the current equilibrium state already calculated. In many cases, the ADD_INITIAL_EQUILIBRIUM command is not required for MAP commands either.
In order to calculate a simple property diagram, set the equilibrium conditions and the stepping axis variable, and then give a STEP command. But if a phase diagram has disconnected lines, ADD_INITIAL_EQUILIBRIUM may still be needed to add two or more initial equilibria so that the MAP calculation starts from such initial equilibria at the specified directions to find all phase boundary lines.
This command becomes unnecessary in most cases, as the MAP routines that by default use the Global Minimization Technique handles all the initial equilibrium points in a robust and automatic way. Therefore, you do not need to have a good guess of the starting point and to add any initial equilibrium point prior to the calculations of various types of phase diagrams. However, if preferred, use this command to add any initial equilibrium points in certain directions, the POLY module uses the specified starting point(s) and corresponding user-added initial equilibrium point(s) for a mapping calculation; in such a case, the automatic procedure in the rewritten MAP routines are not enforced while the Global Minimization Technique can still be applied. Similar situations occur when a specific database, in which there are some definitions of initial-equilibrium adding direction(s) for accessed binary and/or ternary subsystems in its ACCESSED_SYSTEM section, is used in the BINARY or TERNARY module for calculating a binary or ternary phase diagram.
All the initial equilibrium points generated by ADD_INITIAL_EQUILIBRIUM (previously and presently; saved in the current POLY workspace) can be easily listed on screen with LIST_INITIAL_EQUILIBRIA. A certain initial equilibrium point (including its conditions and equilibrium results) can be loaded into the current equilibrium, if needed, with LOAD_INITIAL_EQUILIBRIUM. Any specific or all of the initial equilibrium points can be deleted from current POLY workspace, if desired, with DELETE_INITIAL_EQUILIBRIUM.