Segregation and Microstructure in Continuous Casting Shell Brian G. Thomas and Young-Mok Won Department of Mechanical & Industrial Engineering University of Illinois at Urbana-Champaign September 25, 2000 1
Objectives Develop a fast, simple microsegregation model for the solidification of multicomponent steel alloys. Implement model into other macroscopic models such as heat flow and thermal-stress analysis (CON1D and CON2D) and apply models to continuous casting. 2
Simple Microsegregation Model -based on the Clyne-Kurz model; -extended to account for multiple components, columnar dendrite microstructure, coarsening and /γ transformation. = 1 1 βk 1 C L C o (1 βk)/(k 1) where β = 2α + 1 exp 1 α + exp 1, 2α + α + = 2(α + α C ) and α C = 0.1 α = D S t f X 2 t f = T liq T sol C R 3
Assumption for Simple Model 1. Complete diffusion in the liquid phase. 2. Local equilibrium at the solid-liquid interface. 3. The equilibrium partition coefficient of solute elements applies at the solid-liquid interface and is constant throughout solidification. 4. Nucleation undercooling effects are negligible. 5. Fluid flow effects are negligible. 4
λ SDAS, µm 1000 800 600 400 200 Secondary Dendrite Arm Spacing Model λ SDAS (µm) = (169.1-720.9 C C ) C R -0.4835 for C C 0.15 [1] [2] [1] Measured datas [2] = 143.9 C R -0.3616 C C (0.5501-1.996C) for C C > 0.15 0.03~0.1 o C/sec ~0.2 ~0.5 ~2.0 ~5.0 ~1200 [2] 0 0.0 [3] 0.2 0.4 0.6 0.8 1.0 [2] Carbon Content, wt% 0.03 o C/sec [4] 0.1 0.2 0.5 5 [5] 2.0 5.0 C R = cooling rate ( o C/sec) C C = carbon content (wt%) References 1. H. Jacobi et al. : Steel Res., 1999, v70, p. 357. 2. B. Weisgerber et al. : Steel Res., 1999, v70, p. 362. 3. M. Imagumbai et al. : ISIJ Int., 1994, v34, p.574. 4. D. Senk et al. : Steel Res., 1999, v70, p. 368. 5. A. Suzuki et al.: Nippon Kinzoku Gakkaishi, 1968, v32, p.1301.
Multicomponent Alloy Effect * Liquid temperature for a given liquid composition at the solidliquid interface T( o C) = 1536 78 %C 7.6 %Si 4.9 %Mn 34.4 %P 38 %S where %X (X = C, Si, Mn, P, S) is the liquid concentration at the solid-liquid interface. * Initial guess (equilibrium solidus temperature) => Solidus temperature, T sol, is given when = 1.0 6
Peritectic Phase Transformation Effect * Starting temperature of d/g transformation / T γ start ( o C) = T Ar4 = 1392 +1122 %C 60 %Si +12 %Mn 140 %P 160 %S where / %X (X = C,Si,Mn, P,S) = k L X C L,X * Ending temperature of d/g transformation C L,C 0.53wt%C * Solid fraction of d- and g-phase in the solid = / f γ end / f γ / γ end f start 2 * Average liquid concentration during the d/g transformation and γ = C ave L, i = C L, i + γ γ C L,i 7
Equilibrium Partition Coefficient and Diffusion Coefficient oolute Elements Element k d/l k g/l D d (cm 2 /sec) D g (cm 2 /sec) C 0.19 0.34 0.0127exp(-19450/RT) 0.0761exp(-32160/RT) Si 0.77 0.52 8.0exp(-59500/RT) 0.3exp(-60100/RT) Mn 0.76 0.78 0.76exp(-53640/RT) 0.055exp(-59600/RT) P 0.23 0.13 2.9exp(-55000/RT) 0.01exp(-43700/RT) S 0.05 0.035 4.56exp(-51300/RT) 2.4exp(-53400/RT) * R is gas constant in cal/mol and T is temperature in K. 8
1-D D Finite-Difference Model for Microsegregation local equilibrium domain A j L j 2 γ liquid 1 j = N = 100 d g L Solute diffusion calculation with FDM x j complete mixing * Diffusion equation and initial and boundary conditions for calculation. C S,i t = x D (T) C S,i S,i x C I. C. C S, i = k S / L S,i C o,i at t = 0 B.C. x = 0 at x = 0, λ / 2 SDAS 9
Validation (Case 1) Eutectic Fraction, % 10 8 6 4 2 Solidification time, sec 1000 100 10 1 0.1 measured data[1] numerical prediction[2] simple model 0 0.01 0.1 1 10 100 1000 Cooling Rate, o C/sec Al-4.9%Cu alloy system C L = 33.2%Cu k = 0.145 D S (cm 2 /sec) = 5 10-9 λ SDAS (µm) =46.6 C R -0.29 T liq ( o C) = 660-3.374 C o References 1. Sareal et al. : Metall. Trans. B, 1986, v17a, pp. 2063-73. 2. Voller et al. : Metall. Tarns. A, 1999, v30a, pp. 2183-89. 10
Validation (Case 2) 0.13%C-0.35%Si-1.52%Mn-0.016%P-0.002%S Mn in Liquid Phase, wt% 3.0 2.7 2.4 2.1 C R = 0.045 o C/sec Scheil Clyne-Kurz & Ohnaka(2α) Ohnaka(4α) simple model present FDM Brody-Flemings Lever rule C in Liquid Phase, wt% 1.0 0.8 0.6 0.4 C R = 0.25 o C/sec Scheil Clyne-Kurz & Ohnaka(2α) Ohnaka(4α) & simple model present FDM Lever rule Brody-Flemings 1.8 α = 0.3738 k = 0.77 C o = 1.52 1.5 0.5 0.6 0.7 0.8 0.9 1.0 Solid Fraction 0.0 0.5 0.6 0.7 0.8 0.9 1.0 11 0.2 Solid Fraction α = 3.773 k = 0.19 C o = 0.13
Validation (Case 3) 0.13%C-0.35%Si-1.52%Mn-0.016%P-0.002%S 2.4 0.10 Mn in Liquid Phase, wt% 2.2 2.0 1.8 1.6 0.045 0.25 cooling rate ( o C/sec) measured data[1] (360) (100) present FDM (assumed DAS) (360) (100) simple model (assumed DAS) (348) (149) simple model (SDAS in Eq.12) P in Liquid Phase, wt% 0.08 0.06 0.04 0.02 0.045 0.25 cooling rate ( o C/sec) measured data[1] (360) (100) present FDM (assumed DAS) (360) (100) simple model (assumed DAS) (348) (149) simple model (SDAS in Eq.12) 1.4 0.0 0.2 0.4 0.6 0.8 1.0 Solid Fraction 0.00 0.0 0.2 0.4 0.6 0.8 1.0 Solid Fraction Reference 1. T. Matsumiya et al. : Trans. ISIJ, 1984, v24, pp. 873-82. 12
Validation (Case 4) Temperature, o C 1550 1500 1450 1400 0.015%Si-1.05%Mn-0.0009%P-0.0008%S C R = 0.17 o C/sec +γ γ +L +γ+l γ+l Liquid =0.0 =0.75 ZST(exp.)[1] ZDT(exp.)[1] 1350 Fe-C equilibrium diagram present FDM =1.0 simple model 1300 0.0 0.2 0.4 0.6 0.8 1.0 Carbon Content, wt% Temperature, o C 1550 1500 1450 1400 0.34%Si-1.52%Mn-0.012%P-0.015%S C R = 10.0 o C/sec +γ +L γ +γ+l γ+l Liquid =0.0 =0.75 ZST(exp.)[2] 1350 ZDT(exp.)[2] Fe-C equilibrium diagram present FDM simple model =1.0 1300 0.0 0.2 0.4 0.6 0.8 1.0 Carbon Content, wt% Reference 1. G. Shin et al. : Tetsu-to-Hagane, 1992, v78, pp. 587-93. 2. E. Schmidtmann et al. : Arch. Eisenhuttenwes., 1983, v54, pp. 357-62 13
Validation (Case 5) Calculated Temperature, o C 1550 1500 1450 1400 1350 1300 ZST (55 points) ZDT (53 points) Liquidus (48 points) /γ transformation (30 points) Solidus (47 points) 1250 1250 1300 1350 1400 1450 1500 1550 Experimental Temperature, o C 1. Hot Tensile Tests - Zero Strength Temp.[1-4] - Zero Ductility Temp.[1-5] 2. Differential Thermal Analysis - Liquidus Temp.[6-9] - Solidus Temp.[6-9] - peritectic Temp. [6-8] References 1. G. Shin et al. : Tetsu-to-Hagane, 1992, v78, p. 587.. 2. D. J. Seol et al : ISIJ Int., 2000, v40, p.356. 3. E. Schmidtmann et al : Acrh. Eisenhuttenwe., 1983, v54, p. 357. 4. T. Nakagawa et al. : ISIJ Int., v35, p. 723. 5. H. G. Suzuki et al. : Trans. ISIJ, 1984, v24, p. 54. 6. A Guide to the Solidification oteels, 1977. 7. L. Ericson : Scand. J. Metall., 1977, v6, p. 116. 8. S. Kobayashi : Trans. ISIJ, 1988, v28, p. 535. 9. POSCO data. 14
Effects of C R and l SDAS on Segregation (by Simple Model) Effect of C R ( l SDAS = const.) Effect of l SDAS ( C R = const.) Phase Fraction 1.0 0.8 0.6 0.4 0.2 0.0 constant λ SDAS, (µm) C R, ( o C/sec) 1 10 100 γ L γ 0.18C (45.1) 0.8C (λ SDAS =79.0) L 0.044C (44.1) L Phase Fraction 1.0 0.8 0.6 0.4 0.2 0.0 constant C R =10 o C/sec λ SDAS, (µm) 0.044C 0.18C 0.8C 137.4 103.7 182.0 44.1 45.1 79.0 14.2 19.6 34.4 0.044C L γ L γ L 0.18C 0.8C 1200 1250 1300 1350 1400 1450 1500 1550 Temperature, o C 1200 1250 1300 1350 1400 1450 1500 1550 Temperature, o C * Composition :0.044%C, 0.18%C, 0.8%C (-0.34%Si-1.52%Mn-0.012%P-0.015%S) 15
Combined Effects of C R and l SDAS on Segregation By Simple Model By Finite Difference Model 1.0 0.044C 1.0 0.044C Phase Fraction 0.8 0.6 0.4 0.2 0.0 using λ SDAS in Eq. (12) C R, ( o C/sec) 1 10 100 γ L L γ 0.18C 0.8C L Phase Fraction 0.8 0.6 0.4 0.2 0.0 using λ SDAS in Eq. (12) C R, ( o C/sec) 1 10 100 γ L L γ 0.18C 0.8C L 1200 1250 1300 1350 1400 1450 1500 1550 Temperature, o C 1200 1250 1300 1350 1400 1450 1500 1550 Temperature, o C * Composition :0.044%C, 0.18%C, 0.8%C (-0.34%Si-1.52%Mn-0.012%P-0.015%S) 16
Calculated Solidus Temperatures using Simple Model for Plain Carbon Steels 0.044wt%C 0.18wt%C 0.8wt%C C R λ SDAS T sol λ SDAS T sol λ SDAS T sol 1. Constant secondary dendrite arm spacing 1 44.1 1491.00 45.1 1455.64 79.0 1308.44 10 44.1 1487.86 45.1 1447.13 79.0 1287.40 100 44.1 1478.39 45.1 1428.13 79.0 1234.14 2. Constant cooling rate 10 137.4 1478.55 130.7 1434.06 182.0 1254.72 10 44.1 1487.86 45.1 1447.13 79.0 1287.40 10 14.2 1490.99 19.6 1454.00 34.4 1304.75 3. Combined effects of cooling rate and secondary dendrite arm spacing 1 137.4 1487.93 130.7 1450.38 182.0 1295.43 10 44.1 1487.86 45.1 1447.13 79.0 1287.40 100 14.2 1487.78 19.6 1442.83 34.4 1277.52 17
Non-equilibrium Phase Diagram Calculated with Simple Model 1550 0.34%Si-1.52%Mn-0.012%P-0.015%S Temperature, o C 1500 1450 1400 1350 1300 +L γ cooling rate ( o C/sec) 1 10 100 +γ+l 0.0 0.2 0.4 0.6 0.8 Carbon Content, wt% 18 γ+l Liquid =1.0 =0.0 =0.75 =0.9
Conclusions 1. A simple microsegregation model based on the Clyne-Kurz model developed. 2. A new equation for l SDAS proposed. l SDAS (mm) = (169.1-720.9 C C ) C -0.4835 R for C C 0.15 = 143.9 C -0.3616 R C (0.5501-1.996C) C for C C > 0.15 3. T sol is lowered significantly with independent increases in either C R or l SDAS. 4. The effect of C R less than 100 o C/sec on phase fraction evolution is insignificant in low alloy steels with less than 0.1wt%C, or for phase fractions below 0.9 in other steels 5. Phosphorus and sulfur have a significant effect on solidus temperature due to their enhanced segregation near the final stage of solidification. 6. The simple analytical model presented easily and efficiently incorporate microsegregation phenomena into solidification calculations for use in advanced macroscopic models. 19
Applications with CON1D * Conditions for Calculation - Steel compositions : %C-0.34%Si-1.52%Mn-0.012%P-0.015%S - Casting speed : 1.524 m/min - Slab dimension : 960 mm * 132.1 mm - Working mold length : 1096 mm - Superheat : 1 o C * Predicted T liq and T sol at strand surface using simple model (after iteration) C (wt%) 0.003 0.044 0.1 0.18 0.3 0.44 0.6 0.8 t f (sec) 0.303 0.350 0.437 0.575 0.846 1.268 1.862 2.862 C R ( o C/sec) 73.5 96.6 114.1 120.5 113.0 97.8 82.5 65.1 l SDAS (mm) 20.0 14.4 9.36 18.3 27.6 35.9 40.6 40.2 T liq ( o C) 1524.8 1521.6 1517.2 1511.0 1501.6 1490.7 1478.2 1462.6 T sol ( o C) 1502.5 1487.7 1467.3 1441.7 1406.0 1366.7 1324.5 1276.4 20
Prediction of Phase Fraction (at 10 mm from surface) Solid Fraction 1.0 0.8 0.6 0.4 0.2 0.0 Simple FDM 0.003C 0.044C 0.1C 0.18C 0.3C 0.44C 0.6C 0.8C 1250 1300 1350 1400 1450 1500 1550 Temperature, o C Phase Fraction 1.0 0.8 0.6 0.4 0.2 0.0 Simple FDM 0.1C 0.18C 0.3C 0.44C γ L γ L L γ γ L 1375 1400 1425 1450 1475 1500 1525 Temperature, o C 21
Prediction of t f, C R and l SDAS Profiles 40 100 Solidification Time, sec 30 20 10 0.003C 0.044C 0.1C 0.18C 0.3C 0.44C 0.6C 0.8C 0 0 4 8 12 16 20 Average Cooling Rate, o C/sec 10 0.003C 0.044C 0.1C 0.18C 0.3C 0.44C 0.6C 0.8C 1 0 5 10 15 20 Distance from Surface, mm Distance from Surface, mm 160 Seconary Dendrite Arm Spacing, µm 140 120 100 80 60 40 20 0 0 5 10 15 20 22 0.003C 0.044C 0.1C 0.18C Distance from Surface, mm 0.3C 0.44C 0.6C 0.8C
Predicted Effect of Carbon Content on T sol Variation through Shell Thickness at Mold Exit 1502.7 0.003C 1487.9 0.044C 1467.5 0.1C T sol, o C 1502.6 1502.5 1502.4 0.18 o C T sol, o C 1487.8 1487.7 0.139 o C T sol, o C 1467.4 1467.3 0.132 o C 1502.3 0 5 10 15 20 T sol, o C 1450 1448 1446 1444 1442 0.18C Distance from Surface, mm 7.354 o C 1440 0 5 10 15 20 Distance from Surface, mm 1340 T sol, o C 1335 1330 1325 0.6C T sol, o C 11.622 o C 1320 0 5 10 15 20 Distance from Surface, mm 1487.6 0 5 10 15 20 1416 0.3C 1412 1408 1404 Distance from Surface, mm 9.783 o C 0 5 10 15 20 Distance from Surface, mm 1288 0.8C 0 5 10 15 20 23 T sol, o C 1284 1280 1276 T sol, o C 1380 1376 1372 1368 11.361 o C Distance from Surface, mm 1467.2 0 5 10 15 20 0.44C Distance from Surface, mm 11.097 o C 1364 0 5 10 15 20 Distance from Surface, mm
Microsegregation Model Implemented into CON1D new outputs: t f, C R, l SDAS and T sol profiles Future Work (1) Maximum C L,i (composition between dendrites) and minimum C S,i (composition at dendrite trunks) (2) Non-equilibrium phase diagrams (3) Stainless steel (4) Macrosegregation 24