no way to compare when less than two revisions

Differences

This shows you the differences between two versions of the page.


Previous revision
Next revision
tutorials:t11 [2018/08/01 17:51] – [Step 4: Add the equilibrium solidification path] pwarczok
Line 1: Line 1:
 +===== T11: Simulation of solidification of 0.7C 3Mn steel =====
 +
 +//This tutorial was tested on\\
 +MatCalc version 6.01 rel 1.003\\
 +license: free\\
 +database: mc_fe.tdb//
 +
 +==== Complimentary files ====
 +
 +Click {{:tutorials:t11:script:t11_2013.mcs|here}} to view the script for this tutorial. 
 +
 + 
 +
 +==== Contents: ====
 +
 +  * Scheil calculation
 +  * Back-Diffusion
 +  * Composition set
 +  * Solid-solid transformation
 +
 +The **Scheil-Gulliver method** allows calculating the **fraction and composition of all phases during solidification** step by step from the liquidus temperature to the temperature where solidification of the residual liquid phase occurs ((**References**\\
 +[1] E. Kozeschnik, W. Rindler and B. Buchmayr, //„Scheil-Gulliver simulation with partial redistribution of fast diffusers and simultaneous solid-solid phase transformations“//, Int. J. Mater. Res., 98 (9), 2007, 826-831.\\
 +[2] W. Rindler, E. Kozeschnik and B. Buchmayr, //“Computer simulation of the brittle temperature range (BTR) for hot cracking in steels”//, Steel Res., 2000, 71 (11), 460-465.\\
 +[3] E. Kozeschnik, //“A Scheil-Gulliver Model with Back-Diffusion Applied to the Micro Segregation of Chromium in Fe-Cr-C Alloys”//, Met. Mater. Trans., 2000, 31A, 1682-1684.)). The particular residual fraction at final solidification is dependent on the cooling rate. Generally, the fraction is higher, the higher the cooling rate is. At slow cooling rates, the liquid film can become very thin and the residual liquid enriches more than in the case of higher fraction residual liquid at solidification.\\
 +According to the Scheil-Gulliver hypothesis, illustrated in the figure below (Fig.1), a solidifying liquid with the initial composition C<sub>0</sub> is slightly undercooled. Consequently, a certain fraction of solid is formed with a composition of C<sub>S,1</sub> with the remaining liquid phase with the composition C<sub>L,1</sub>. Both phases are in **local equilibrium following the lever rule**. From that instant on, the composition of the solid phase with regard to the substitutional atoms is frozen due to the Scheil-Gulliver assumption of very slow (and therefore negligible) diffusion in the solid phase. The interstitial elements boron, carbon and nitrogen are assumed to be nevertheless highly mobile in the liquid phase as well as the solid phase. Therefore, these elements are assumed to be able to establish equilibrium between the liquid and solid phases due to back-diffusion of enriched solute atoms from the liquid phase into the solid phase. This process occurs repeatedly, enriching the solid phase as well as the liquid phase with solute during cooling (Fig.1:C<sub>L,i</sub> and C<sub>S,i</sub>).
 +
 +{{:tutorials:t11:img:t11_figure1v2.png|MatCalc figure}}
 +
 +''Fig.1: Scheme of Scheil solidification of a hypothetic Fe-C alloy. During solidification the actual liquid phase, beginning with C0, is undercooled and solidifies according to the lever rule enriching the liquid phase with solute atoms''
 + 
 +When solidification is finished, the microstructure of the solid phase shows a composition gradient from the regions where the solidification started (centre of the dendrite) to the areas where the last liquid solidified (outer shell of the dendrite). The composition of the residual liquid at final solidification corresponds to the composition of the interdendritic regions. It should be noted that the Scheil-Gulliver analysis yields an upper limit of segregation since the substitutional elements also have a finite mobility and have some potential of back-diffusion.\\
 +Scheil-type calculations in multi-component systems have proven to be a useful method to simulate solidification processes. The present example shows how to carry out this type of simulation with MatCalc and demonstrates the effect of carbon back-diffusion on the solidification process. Moreover, the influence of solid-solid phase transformations, i.e. the peritectic transformation, is analyzed.
 + 
 +===== Step 1: Setup the thermodynamic system (see also Tutorial T2) =====
 +
 +Create a new workspace file. From a suitable database (**mc_fe.tdb**) define the elements **Fe**, **Mn** and **C** and the phases **liquid**, **BCC_A2** (ferrite), **FCC_A1** (austenite) and **Cementite**.
 +
 +{{:tutorials:t11:img:t11_database_6011003.png?650| MatCalc Databases}}
 +
 +Enter the system composition in weight percent as listed in the subsequent figure selecting **'Global Composition ...'** or pressing the **F7** key.
 +
 +{{:tutorials:t11:img:t11_composition_2016.png| MatCalc Composition}}
 +
 +Set initial values with **'Calc'->'Set start values'** or **Ctrl+Shift+F**. Calculate equilibrium at 1600°C. The results in the **'Phase summary'** window are
 + 
 +<code>
 +LIQUID  *        act 1,00000e+000 dfm: +0,00000e+000
 +### inactive ###
 +BCC_A2        - OK - 0,00000e+000  dfm: -1,02368e+003
 +FCC_A1        - OK - 0,00000e+000  dfm: -1,08270e+003
 +CEMENTITE    - OK - 0,00000e+000  dfm: -1,28107e+004
 +</code>
 + 
 +
 +===== Step 2: Carry out a Scheil calculation =====
 +
 +Since we expect austenite to be the first phase to form on solidification, let us first look for the solubility temperature of this phase. Select **'Calc - Search phase boundary ...'** or press **Ctrl+Shift+T** and select **FCC_A1** as target phase.
 +
 +{{:tutorials:t11:img:t11_search_boundary_2016.png| MatCalc search boundary}}
 +
 +Press **'Go'**:
 + 
 +<code>Tsol 'FCC_A1': 1471,09 C (1744,24 K) iter: 4, time used: 0,02 s</code>
 + 
 +The first solid phase becomes stable at **1471.09°C**. So let us start with the Scheil simulation at **1500°C** and go down to **1000°C** in steps of **10**. Open the Scheil-calculation dialog with **'Calc - Scheil calculation ...'** or press **Ctrl+Alt+T**.
 +
 +{{:tutorials:t11:img:t11_scheil_calculation_6011003.png| MatCalc scheil calculation}}
 +
 +
 +Press **'Go'** to start the simulation. In order to be able to remove the equilibrium content of solid phases from the system after each temperature step, MatCalc needs to create **copies of all phases except the dependent** (liquid) **phase**. MatCalc therefore asks
 +
 +{{:tutorials:t11:img:t11_warning_2016.png| MatCalc warning}}
 +
 +Click **'Yes'**. MatCalc then creates 3 new phases with the name of the original phase plus the suffix **'_S'**. This suffix denotes that the corresponding phase is a **'solid'** phase. The result in the **'Output'** window looks as follows
 + 
 +<code>
 +Checking solid phases ... - OK -
 +Searching initial equilibrium ...
 +1, 0,05 s, 1500,00 C (1773,15 K), its 2, f=1,00000000,  LIQUID
 +...
 +3, 0,05 s, 1480,00 C (1753,15 K), its 2, f=1,00000000,  LIQUID
 +4, 0,05 s, 1470,00 C (1743,15 K), its 5, f=0,97018194,  LIQUID FCC_A1
 +...
 +12, 0,08 s, 1390,00 C (1663,15 K), its 6, f=0,25232525,  LIQUID FCC_A1
 +13, 0,08 s, 1380,00 C (1653,15 K), its 6, f=0,22793997,  LIQUID FCC_A1
 +14, 0,08 s, 1370,00 C (1643,15 K), its 6, f=0,20743835,  LIQUID FCC_A1
 +15, 0,08 s, 1360,00 C (1633,15 K), its 6, f=0,18999540,  LIQUID FCC_A1
 +...
 +39, 0,17 s, 1120,00 C (1393,15 K), its 5, f=0,05685971,  LIQUID FCC_A1
 +40, 0,17 s, 1110,00 C (1383,15 K), its 11, f=0,00000000,  FCC_A1 CEMENTITE
 +fraction of phase 'LIQUID' smaller than 0,01. Calc finished
 +Steps: 41, CalcTime: 0,17 s
 +AktStepVal: 1383,150000
 +- OK -
 +</code>
 +
 +The current value of the fraction liquid at each temperature step is displayed with **'f=xxx'** in each line. The final eutectic is reached at T=1110°C. When the fraction liquid comes **below 0,01 (minimum liquid fraction)**, the dependent phase is dissolved and the calculation is finished.\\
 +Let's now look at the result.\\
 +Create a **XY-data** plot that shows the fraction of residual liquid as a function of temperature. The correct variable is **F$LIQUID**, you have to change the default x-data to **T$C** in the options window to show °C instead of K. Rescale the **x-axis** from **1100°C to 1500°C**. A few other settings were made until the result looks as follows
 +
 +{{:tutorials:t11:img:t11_plot1_result_scheil_2_2016.png?650| MatCalc plot}}
 +
 +===== Step 3: Add a Scheil calculation with back-diffusion of carbon =====
 +
 +The solidification of the current steel as calculated by the Scheil model predicts final solidification at too low temperature. In reality, **carbon atoms are fast enough** not only in the liquid but also in the solid to be able **to equilibrate between the solid and liquid phases**. Therefore, to get more realistic simulation results, we must allow for **back-diffusion of carbon**. First, in order not to loose the results of the previous simulation, rename the current buffer to **'Scheil'** and create a new one with the name **'Scheil with BD of C'**. Open the Scheil calculation dialog with **'Calc - Scheil calculation ...'** or press **Ctrl+Alt+T**. Highlight carbon in the list box and press **'Toggle'**.
 +
 +{{:tutorials:t11:img:t11_scheil_calc_toggle_2016.png| MatCalc scheil calculation}}
 +
 +Press **'Go'** to start the simulation. After each temperature step, MatCalc sets up a paraequilibrium calculation, where all elements without back-diffusion have fixed composition variables and **only the elements with back-diffusion are unconstrained**. By that means, carbon is always brought back into equilibrium with regard to the **solid and liquid** phases after a regular Scheil simulation step is carried out.\\
 +We are going to display the Scheil curve with and without back-diffusion in one window. Therefore, first, lock the previous Scheil calculation, which is still displayed in the diagram window and name it **'Scheil'**.
 +
 +{{:tutorials:t11:img:t11_options_window_2016.png| MatCalc options}}
 +
 +Then, we make sure that the just created buffer (**'Scheil with BD of C'**) is selected as current buffer for the plot in the options window and drag and drop the **'F$LIQUID'** variable into the plot (note that alternatively you can use the 'Lock and duplicate series' option accessible via the right mouse button or the 'View' menu). Rename the new series to **'Scheil with BD of C'**. The new diagram looks like the following figure.
 +
 +{{:tutorials:t11:img:t11_plot2_result_scheil_bd_2_2016.png?650| MatCalc plot}}
 +
 +===== Step 4: Add the equilibrium solidification path =====
 +
 +Let us finally add the solidification temperatures for thermodynamic equilibrium conditions. To speed up the calculations, we **suspend all 'solid' phases which were needed for the Scheil calculations**. Open the **'Phase status'** dialog (**F8**). Set the **'suspended'** flag for each **'solid'** phase.
 +
 +{{:tutorials:t11:img:t11_phase_status_2016.png?650| MatCalc phase status}}
 +
 +
 +Create a new buffer with the name **'Equilibrium'**. Open the stepped calculation dialog (**Ctrl+T**). Calculate a stepped calculation with temperature as the variable between 1500°C and 1350°C in steps of 5. Press **'Go'** to start the calculation.
 +
 +{{:tutorials:t11:img:t11_step_equilibrium_6011003.png| MatCalc step equilibrium}}
 +
 +
 +Lock the last series ('Scheil with BD of C'). Change the buffer to be used in the diagram window to **'Equilibrium'**and drag and drop the **'F$LIQUID'** variable into the plot. Rename the current series and the diagram looks like this:
 +
 +{{:tutorials:t11:img:t11_plot3_scheil_calc_final_2_2016.png?650| MatCalc plot}}
 +
 +
 +On changing the **'y-axis'** type to **'log'** and scaling from 0.01, the results finally display as
 +
 +{{:tutorials:t11:img:t11_plot4_scheil_calc_final_log_2_2016.png?650| MatCalc plot}}
 +
 +
 +From this diagram, and with the assumption that all residual liquid freezes (solidifies) when 1% residual liquid is reached, the predicted solidification temperatures using the classical Scheil model, Scheil with back-diffusion and full thermodynamic equilibrium can be read as 1110, 1335 and 1380°C, respectively. Probably, 1335°C comes closest to the real, experimentally observed solidification temperature for this alloy composition.
 +
 +
 +
 +===== Consecutive articles =====
 +
 +The tutorial is continued in article [[tutorials:T12 | T12 - Using the MatCalc console]]
 +
 +Go to [[:tutorials|MatCalc tutorial index]].
  
tutorials/t11.txt · Last modified: 2023/08/02 22:09 by pwarczok
 
Recent changes RSS feed Donate Powered by PHP Valid XHTML 1.0 Valid CSS Driven by DokuWiki