Abstract
Since CaO-Al2O3-based mold fluxes are one of the most important mold flux systems in metallurgic processes, it is important to explore their structure characteristics and viscosity. Molecular dynamics simulation is performed to study the effect of w(CaO)/w(Al2O3) ratio on both the structural and viscosity properties of CaO-Al2O3-based mold fluxes. A systematic analysis of the structure and thermodynamics on CaO-Al2O3-based mold fluxes is carried out, and it is well known that the viscosity of mold fluxes is related to the structure. The results show that the formation of stable structures of Si-O in the mold fluxes was beneficial to reduce the probability of structural interconnection, degree of polymerization, and viscosity of the molten slag. In the cationic structure, the contents of Ca-O-Al and Ca-O-Si are more stable, the interconnection of the Ca-O-Al and Ca-O-Si network weakens, and the viscosity decreases. The tetrahedra [AlO4] and [SiO4] have similar structures, but they exhibit different thermodynamic and physical properties. Viscosity test shows that CaO/Al2O3 = 0.88–2 continuously increased, when the cosolvent content Li2O = 1%–4%, CaO-Al2O3-based mold flux viscosity decreased, the degree of network structure polymerization decreased, and the complex structure depolymerized. Increasing the water content in the cosolvent is beneficial to reduce the viscosity of the crystallizer.
1. Introduction
High-manganese and high-aluminum steel is widely used in military applications, automotive manufacturing, and electrical components due to its excellent mechanical properties [1–5]; the high content of Mn and Al reduces electricity consumption, improves electromagnetic properties, and endows high strength and low density. During conventional casting production, traditional CaO-SiO2-based mold fluxes can cause slab surface defects and interrupt the casting process by reacting with Mn and Al, decreasing SiO2 content and increasing Al2O3 content. As a result, the physical and chemical properties of traditional CaO-SiO2-based mold fluxes are severely deteriorated. Alternatively, low-reactivity CaO-Al2O3-based mold fluxes have been proposed to replace the traditional CaO-SiO2-based mold fluxes [6–10]. In this paper, the structure of CaO-Al2O3-based mold fluxes is studied from the perspective of molecular dynamics.
The thermodynamic properties of CaO-Al2O3-based mold fluxes have been widely studied. For instance, Shi et al. showed that the viscosity increased slightly with the CaO/Al2O3 ratio from 2 to 3, and this increase became significant with further increasing CaO/Al2O3 ratio to 4 [11]. However, the physicochemical properties of mold fluxes have been studied from the viewpoint of thermodynamics regardless of the structure of the mold fluxes. For instance, Wang et al. have studied the influence of B2O3 addition on the physicochemical properties of CaO/Al2O3 mold fluxes, but the structure of mold fluxes has not been taken into account [12]. Zhang et al. have studied the effect of different molten mold fluxes on the viscosity of CaO-Al2O3-based mold fluxes at 1400°C [13]. However, the relationship between physicochemical properties and structure of the mold fluxes has not been explored. Similarly, Wang et al. have studied the influence of CaO/Al2O3 ratio on melting temperature and viscosity of low-reactivity CaO-Al2O3-based mold fluxes [14]. Zhou et al. have studied the viscosity variation of CaO-Al2O3-based molten mold flux in the temperature of 900°C to 1300°C [15]. Shiro et al. have studied the influence of CaO and Al2O3 content on the viscosity of CaO-Al2O3-based mold fluxes in the temperature range of 1550°C to 1650°C [16], but the structure of molding fluxes has not been considered. Some researchers pay attention to the effects of viscosity and structure [17–20]; it should be noted that the macroscopic properties of the continuous casting mold fluxes depend on the structure. Therefore, the effect of the CaO/Al2O3 ratio in a broad range on the viscosity and structure of CaO-Al2O3-based mold fluxes must be widely studied systematically.
Therefore, the current work investigates the structure of CaO-Al2O3-B2O3-Na2O-BaO-SiO2-Li2O mold flux from the viewpoint of the molecular dynamics. Moreover, density functional theory (DFT) calculations are employed to analyze the structure and thermodynamic properties from the viewpoint of quantum mechanics and calculate the bond lengths of Li-O, Al-O, Ca-O, and Si-O. Furthermore, the thermodynamic behavior of [AlO4] and [SiO4] structures is also studied to obtain Gibbs free energy. Moreover, in order to measure the viscosity of mold fluxes, these results theoretically provide novel insights into the structure and viscosity of mold fluxes. The aim of this study is to investigate the effect of CaO/Al2O3 ratio on the CaO-Al2O3-B2O3-Na2O-BaO-SiO2-Li2O mold flux structure and viscosity and to provide further fundamental understanding of the effect of this parameter on flux properties which will be used for developing CaO-Al2O3-B2O3-Na2O-BaO-SiO2-Li2O mold flux.
2. Molecular Dynamics Computation
2.1. Model Establishment of Basic Theory
The vibrational analysis evaluation can be used to compute temperature-dependent enthalpy (H), entropy (S), vibrational free energy (GV), and constant-pressure heat capacity (CP). The heat capacity (CP) can be given as [21]where Ctrans refers to the translational heat capacity, Crot represents the rotational heat capacity, Cvib denotes the vibrational heat capacity, R represents gas constant (8.314 J/mol·k), K corresponds to the Boltzmann constant, h denotes Planck’s constant, T represents the absolute temperature, and H refers to enthalpy [22]. The enthalpy (H) can be given aswhere Htrans refers to the translational enthalpy, Hrot represents the rotational enthalpy, and Hvib denotes the vibrational enthalpy. The entropy (S) can be given aswhere Strans refers to translational entropy, Srot represents the rotational entropy, Svib denotes the vibrational entropy, p denotes pressure, σ indicates symmetry number, c refers to molar concentration, and IA, IB, and IC represent the moments of inertia. can be given aswhere Gv refers to the vibrational free energy, E (0k) represents the total energy, Gproduct refers to the vibrational free energy product, and Greaction corresponds to the vibrational free energy of the reaction.
2.2. Calculation Method
Based on the application of density general letter theory in quantum modernity, this paper establishes the structural calculation model and calculates the physical and chemical properties of the basic skeleton structure and complex structure. At the same time, it calculates the binding energy, total energy, and ∆G of the structure; the structural calculation model calculates the most stable basic structure and complex structure through structural optimization, so the steps establish a calculation model as follows:(1)For stable basic structures and stable complex structures, it is necessary to perform vibration frequency calculations and perform geometric optimization based on statistical thermodynamics to obtain microstructures. DFT (density functional theory) calculations were utilized to compute the structure of the mold fluxes and analyze the thermodynamic properties. Finally, based on optimizing the working vibration frequency of the structure, the structure was obtained by statistical thermodynamic method.(2)The structural calculation model is based on quantum mechanics as the basic theory. Based on the BPW91 system in density functional theory, the basic structure and complex structure are optimized. The model uses DMol3, combines LST (linear synchronous transit)/QST (quadratic synchronous transit) with CGM (conjugate gradient methods), and then adopts the transition method. The most stable structure was searched for.(3)The optimization calculation of the structural calculation model is based on DFT's DMol3 quantum mechanics module, which can accurately obtain the parameters of molecular structure and physical and chemical properties. The stabilized structures were used, as reactants, to explore the transition states via complete LST (linear synchronous transit) and QST (quadratic synchronous transit) methods. The SCF (self-consistent field) method was used with a total energy.
3. Structure Properties
3.1. Analysis of Structures and Thermodynamic Properties
The structure should be studied because the physicochemical properties of CaO-Al2O3-SiO2 slag are inseparable from the structure [13, 23–26]. Here, Si-O and Al-O exist as an aluminosilicate structure. The bond length and bond energy of Al-O, Ca-O, Si-O, and Li-O are summarized in Table 1. In addition, the bond length in the Al-O unit (1.65 Å) is larger than that in the Si-O unit (1.53 Å), the bond length in the Li-O unit (2.05 Å) is larger than that in the Si-O unit (1.53 Å), and the bond energy of Si-O bond (−853.48 kJ/mol) is lower than that of Al-O and Li-O. The lower the binding energy and total energy, the more stable the structure. Therefore, Si-O is more stable than Al-O and Li-O, the basic framework structure of mold flux cosolvent SiO2 is Si-O, and adding the Si-O structure to the slag structure is beneficial to the stability of the physical and chemical properties of the mold slag, thus reducing the possibility of the simple structure in the slag being transformed into a complex structure.
Figure 1 shows the change in ∆G of anionic structure and cationic structure with respect to temperature; the structures become less stable with higher temperatures. The reaction that produces a O-Si-O structure has the greatest ΔG trend, while the reaction that produces Si-O-Si structures has lesser ΔG tendency; therefore, it is easier to generate O-Si-O structure in the mold fluxes.

The thermodynamic behavior of the structure mainly determines the structure of mold fluxes. Figure 2 shows the total energy of the anionic structure and cationic structure. Herein, the results indicate structure diversity and different physical characteristics. The results indicate that O-Al-O is more stable than Al-O-Al. Ca-O-Si structure is stable, it is difficult to form complex structure, and the polymerization of slag is reduced, resulting in weak connection of the molten fluxes. The viscosity of the mold flux decreases with CaO content increasing, because the free O2- ions released by O2- ions can enter cosolvent Al2O3 bond and form the O-Al-O structure, which causes the relatively reduced connections of the structure of the mold flux and the reduction in mold flux viscosity. Figure 3 shows the change in ∆G of different structures with respect to temperature; it can be readily observed that the ∆G increases with increasing temperature, the structures become less stable with higher temperatures, and the increasing temperature reformed the flux network through the breakages of structure consisting complex structure linkages, as well as the reduction in molten slag viscosity. If the structure consisting of complex increases and hinders the depolymerization of the structure, the viscosity of mold fluxes is increased due to complex structure.


3.2. Analysis of Tetrahedral Structure and Thermodynamic Properties
It should be noted that [AlO4] and [SiO4] are the basic tetrahedron structure units of the CaO-Al2O3-based mold fluxes. In the aluminate melt structure, Al2O3 exhibits a tetrahedral structure and octahedral structure, as Al2O3 has amphoteric characteristics in mold flux. It exists as a [AlO4] network structure in alkaline melt. Moreover, the addition of SiO2 promotes the incorporation of Si into Al2O3 structure. The reason is that the pickup of Al2O3 will increase Al-O complex structure and form the bond of Al-O-Si to raise the polymerization degree of the molten slag, which results in greater viscosity flow resistance and increased viscosity. In terms of molten flux, CaO decomposes to form Ca ions and free oxygen, where the former compensates for charge difference between [AlO4] and [SiO4], and the later forms a depolymerization structure with nonbridging oxygen; the viscosity is decreased due to the low polymerization.
As is shown in Table 2, in addition, the bond length in the [AlO4] tetrahedral unit is larger than that in the [SiO4] tetrahedral unit; however, the total energy and bond energy of [SiO4] tetrahedral are smaller than [AlO4], and the bond angle of [AlO4] and [SiO4] structures is 109.47°. The results indicate that Si-O is more stable than Al-O; despite the fact that both [AlO4] and [SiO4] exhibit the same tetrahedral structure and basic network, the stability of [AlO4] tetrahedral network is not as good as [SiO4], and the possibility of [SiO4] structure transition to complex structure is reduced. This indicate that [SiO4] is more stable than [AlO4]; hence, the [SiO4] difficulty forming complex structures, which causes the weaker link of the molten fluxes. As shown in Figure 4, although the structures of the [AlO4] and [SiO4] are similar, [AlO4] requires more thermal energy to depolymerize than [SiO4], both [AlO4] and [SiO4] exhibited different thermodynamic and physical properties, and [SiO4] is more stable than [AlO4]. The effect of aluminosilicate structures becomes important because of high Al2O3 content; in the aluminosilicate structure, the [AlO4] tetrahedron structures increased, the stability of the structures is reduced, and the viscosity is increased due to the high polymerization. Hence, the viscosity is increased due to the high polymerization and hindered ionic mobility.

4. Experimental
4.1. Sample Preparation and Viscosity Measurement
Table 3 shows the chemical composition of the CaO-Al2O3-based mold fluxes before heating; reagent grade powders of CaO, Al2O3, B2O3, Na2O, BaO, SiO2, and Li2O were used as raw materials, and chemical reagents are pure reagents (>95.5%). CaO was calcined at 800°C in a muffle furnace for 8 h, and other powders were placed at 500°C. Then 1wt%, 2wt%, 3wt%, and 4wt% Li2O were added to the CaO-Al2O3-based mold fluxes, and four experimental samples were obtained.
The viscous behavior of the CaO-Al2O3-based mold fluxes with various concentrations of Li2O has been studied using the viscosity measurement method to understand the effects on the viscosity with these additives, and the viscosity of mold fluxes was measured by using a rotating cylindrical viscometer. The schematic diagram is shown in Figure 5. The Mo spindle was rotating at a fixed speed of 15r/min, a Pt/Rh thermocouple was placed right beneath the bottom of a high-purity graphite crucible, 150g flux sample was placed in the crucible and heated to 1350°C, it decreased with a rate of 5°C/min, and the viscosity of experimental flux was continuously recorded.

4.2. Effect of CaO/Al2O3 Ratio on Flux Viscosity
Figure 6 shows a diagram of network changes. As the CaO content increases, the viscosity of the mold flux decreases, which is because the free Ca+ released by CaO can break O-Al-O and Al-O-Al bond and form stable structure the Ca-O-Al and Ca-Al-O structure, which results in a relatively reduced connection of the structure and a decrease in the viscosity of the mold fluxes. The CaO and Al2O3 content is very important to control the viscosity of mold fluxe, while this effect would be attenuated with the excessive addition of CaO and Al2O3, although the viscosity can be increased by adding excessive addition of Al2O3. Figure 7 shows the viscosity of mold fluxes versus w(CaO)/w(Al2O3) in the temperature range from 1200°C to 1325°C and the effect of CaO/Al2O3 ratio on the viscosity of sample A1 to A4. The results show that the viscosity of the slags decreased with the increase in CaO/Al2O3 ratio. It can be seen from Figure 7 that the viscosity of the melt decreased with increasing CaO/Al2O3 ratio from 0.88 to 2. The continuous addition of CaO into the aluminates instigated the weakening of the structures of [AlO4], and it should be noted that the content of Ca+ ions in the mold fluxes is enough for charge compensation. When CaO content increases, the content of Ca+ ions increased and destroyed the structure.


Figure 8 shows the effect of sample on the viscosity of the mold fluxes; the temperature rise reduced the viscosity of the mold fluxes, which reduces the viscosity and facilitates structure depolymerization. Figure 9 shows the viscosity of samples A1 to A4 in the mold fluxes with various Li2O contents. Meanwhile, Li2O could also stabilize the structural unit of Li-O structures by link aluminate. Li2O could release O2- ions to break the complex structure into simple structure. The schematic diagram for the changes of the [AlO6] is shown in Figure 10; the stability of [AlO4] tetrahedral network is not as good as [AlO6], and the possibility of [AlO6] structure transition to complex structure is reduced. Moreover, increasing CaO/Al2O3 = 0.88–2 leads to an increase in the fraction of [AlO6] octahedron but a decrease in the fraction of [AlO4] tetrahedron, and the higher content of CaO destroys the structure of [AlO4] and, consequently, reduces the viscosity. The changes of the [SiO4] are shown in the schematic diagram in Figure 11; the O2- from CaO reacts with the electropositive Si ion from SiO2 and generates stable structure [SiO4], the complex structure transforms into a simple structure, the stability of [AlO4] tetrahedral network is not as good as [SiO4], the possibility of [SiO4] structure transition to complex structure is reduced, the alkaline oxide determines the structural properties of the mold fluxes, which are also influenced by the structure, and the structure of [AlO4] high melting substances leads to the formation of the complex structure, increase the degree of polymerization, and increases the viscosity of the mold fluxes. It should be noted that the viscosity of the mold fluxes is also affected by the polymerization, and the microstructure of the high polymerization can lead to a high viscosity, if it can be structurally depolymerized, resulting in a decrease in viscosity.




As can be seen, as the CaO/Al2O3 content decreases, the viscosity of the slag gradually increases, and as the Li2O content decreases, the viscosity of the slag increases. Within the experimental temperature range of 1225°C to 1300°C, the viscosity of the slag decreased with addition of Li2O. The Li2O = 1–2 significantly lowered the viscosity, and further addition of Li2O = 2–4 significantly lowered the viscosity. The change in Li2O concentration releases Li ions from the mold fluxes, the complex structure is transformed into a simple structure, and the viscosity of mold fluxes is influenced by the complex anions. Therefore, the viscosity of the mold fluxes can be controlled by adjusting the polymerization of the structure. It should be noted that the dissociation of simple structures increases the fluidity of the melt, as excessive O2- does not have a significant effect on the less polymerized network structure but rather leads to an increase in viscosity.
4.3. Viscous Flow Activation Energy
The mold flux temperature is higher than the break temperature, which can be considered as a Newtonian fluid, so the viscosity changes with temperature follow the Arrhenius equation:where T is temperature, K; A is constant; Ea is activation energy, kJ/mol; and R is perfect gas constant, 8.314 J/mol·K.
The optical basicity (Λ) can characterize the ‘usability’ of providing O2-, which explains the degree of polymerization of the mold fluxes. The corrected optical basicity (Λcorr) was proposed by Mills to charge balance the Al3+ ions or the Si4+ chain or ring [28]. In the calculation shown in formula (7), corrected optical basicity proposed by Mills (Λcorr), Λ is the optical basicity;X is the mole fraction; n is the number of slag oxygen atoms in component molecule; and Λcorr is the corrected optical basicity.
Formula (6) shows the linear relationship between lnη and 1/T in the Newtonian fluid region, to obtain activation energy. Table 4 shows change trend of melting temperature and break temperature of slags A1, A2, A3, and A4 corresponding to the change trend of activation energy. In addition, the molten slag is in a completely liquid state and is a Newtonian fluid, and the activation energy can be estimated using the Arrhenius equation. Break temperature is closely related to activation energy and bond breakage of composite negative ions in mold fluxes at elevated temperature; the polymerization degree decreases and the migration energy increases. The viscosity of the mold fluxes is usually reduced as the corrected optical alkalinity increases. This is higher than those by Li al [27]. In the experimental research of mold fluxes, the different network structures of the mold fluxes of other researchers and authors lead to different activation energy.
5. Conclusions
This work investigated the variations of structure, thermodynamic properties, and viscosity of CaO-Al2O3-based mold fluxes. The following conclusions can be drawn:(1)The Si-O is more stable than Li-O, and Al-O is not easily linked to other structures, relatively weakening the mold flux structure and lowering the degree of polymerization. The B-O is more stable than Al-O, which is not easily linked to other structures, relatively weakening the mold flux structure and lowering the degree of polymerization, resulting in relatively reduced connections of the structure of the mold flux and the reduction in mold flux viscosity. While this effect would be attenuated with the excessive addition of Al2O3, the complex structure is more likely to exist in the CaO-Al2O3-based mold fluxes.(2)The viscous behavior of the CaO-Al2O3-based mold flux system with various concentrations of Li2O has been studied using the viscosity measurement method. The increased Li2O in the mold fluxes destroys the network structure and reduces the polymerization degree and the transition from complex structure to simple structure. The viscosity of the mold fluxes usually decreases with the increase of corrected optical basicity. This report shows the importance of exploring and establishing the relationship between the performance and structure of the mold fluxes from a structural perspective.
Data Availability
The data (figures) that support the findings of this study are available upon request from the corresponding author.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors are grateful for the financial support by the National Key R&D Project of China (no. 2017YFC0805100) and Science and Technology Projects of Liaoning Province (no. 2018307003). The project was supported by the State Key Laboratory of Metallic Materials for Offshore Equipment and Applications (no. 2017-01-04).