Next Article in Journal
Influence of Strengthened Nodes on the Mechanical Performance of Aeolian Sand–Geogrid Interface
Next Article in Special Issue
Numerical Modeling of Transient Flow Characteristics on the Top Surface of a Steel Slab Continuous Casting Strand Using a Large Eddy Simulation Combined with Volume of Fluid Model
Previous Article in Journal
Microstructural Characterization and Prior Particle Boundary (PPB) of PM Nickel-Based Superalloys by Spark Plasma Sintering (SPS)
Previous Article in Special Issue
The Effect of Mold Structure and Cooling Parameters on Heat Transfer during Billet High-Speed Continuous Casting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of Spatial Distribution of Multi-Size Bubbles in a Slab Continuous-Casting Mold Water Model

1
School of Metallurgical and Ecological Engineering, University of Science and Technology Beijing, Beijing 100083, China
2
National Engineering Research Center of Continuous Casting Technology, Central Iron and Steel Research Institute, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Materials 2023, 16(13), 4666; https://doi.org/10.3390/ma16134666
Submission received: 23 February 2023 / Revised: 22 May 2023 / Accepted: 24 June 2023 / Published: 28 June 2023
(This article belongs to the Special Issue Metallurgical Process Simulation and Optimization2nd Volume)

Abstract

:
In this paper, a fully coupled large eddy simulation model, including the volume of fluid model, the discrete phase model, the bubble-collision model, and the bubble-breakup model was used to simulate the spatial distribution of multi-size bubbles and its impact on the instantaneous two-phase flow in a slab continuous-casting mold. The influence of the bubble-interaction model on the bubbles’ three-dimensional spatial distribution and size distribution, as well as on two-phase flow was discussed. By comparison with the velocity on the meniscus and the average bubble diameter inside a continuous-casting slab water model, the appropriate numerical model was recommended to accurately simulate the two-phase flow and characteristics of discrete bubbles. The submerged entry nozzle and the area around it saw bubble coalescence and breakage more frequently than other areas. The key interaction between the bubbles was their bouncing in the deep region of the mold. In the mold, the average bubble diameter was 0.741 mm, and 44.5% of the total number of bubbles had an approximately average diameter.

1. Introduction

The discrete argon bubbles in a continuous-casting (CC) mold form surface defects after they are entrapped in the solidification front [1]. In addition, the disturbance caused by the bubbles close to the meniscus increases the probability of slag entrainment [2]. However, argon gas injection is still required during the slab CC process as one of the crucial measures to avoid nozzle clogging [3]. Therefore, the current study focuses on improving the multiphase flow (steel, slag, and air) [4] inside the mold after the argon gas injection, as shown in Figure 1. The treatment of argon gas bubbles inside the CC mold is a key and difficult point to accurately calculate the multiphase flow. One widely used method is to treat argon bubbles as the Eulerian phase, similar to molten steel. For example, Liu established a Eulerian–Eulerian (E–E) two-fluid model to investigate the vortex flow in the mold with argon gas injection [5]. Bai employed the E–E multiphase model to simulate turbulent flow with the injection of argon bubbles in a slide-gate tundish nozzle [6]. Using the E–E model, Sánchez-Pérez et al. [7] discussed the distinction between the two-phase flows with the coupled and uncoupled model in a slab mold. Singh [8] developed a three-dimensional E–E mathematical model to study bubble movement in a water model with various parameters. That model successfully predicted the influence of argon gas injection on the multiphase flow but ignored the effect of the change in bubble diameter. Recently, the homogeneous multi-site group (MUSIG) model [9,10] and the non-homogeneous MUSIG model [11] were applied in an E–E approach to include the size distribution of bubbles. The E–E approach extended with the multi-size model was useful in predicting the volume fraction of bubbles and its influence on the flow field. However, it is difficult to calculate the influence on meniscus fluctuation, especially slag entrainment.
Thus, many studies [2,12,13,14,15,16] calculate each argon bubble as a separate individual using a Lagrangian method. Two-way coupled interphase forces were used to actualize the interactions between discrete bubbles and molten steel [17]. Chen et al. systematically studied the bubble distribution and its influence on the flow field using the two-way coupled discrete-phase model (DPM) [2,15,16,18], as well as the entrapment of inclusions using the one-way coupled DPM [19]. The phase interface between the molten steel and the slag phase (or air phase) was determined using the volume of fluid (VOF) model. Wang [20] combined the VOF model, k-ε model, solidification model, and DPM to investigate the impact of the injection gas’s flow rate on the two-phase flow and solidification of steel in a slab CC mold. Yang and Zhang [21,22] used the modified bubble-coalescence-and-breakup model to consider the effect of the bubble diameter based on the above models. When applying DPM to solve the bubble motion, either a Reynolds-averaged Navier–Stokes (RANS) or a large eddy simulation (LES) turbulence model are required for the steel (slag and air) phase. The LES model was recommended to calculate the turbulent flow in a submerged entry nozzle (SEN) and mold due to the presence of asymmetric flow in the mold [23,24]. However, unfortunately, few studies have considered both the collision and the breakup of bubbles and the transient multiphase flow at the same time due to the complexity of the model.
In order to study the spatial and size distribution of bubbles and two-phase flow during the slab CC process, a three-dimensional numerical model including the LES model, VOF model, DPM, bubble-collision model, and bubble-breakup model was constructed. Three bubble-breakup models were compared with the measured bubble diameter and meniscus speed. The main aims of this paper were to propose and establish a mathematical model that considers both the transient multiphase flow inside the CC mold and the coalescence, bounce, and breakup of the bubbles.

2. Mathematical Model

2.1. Governing Equations

The water and air phases were treated as incompressible Newtonian fluids and the turbulent flow inside the slab CC mold was solved using the LES turbulent model [25]. The water–air-free surface, and the volume fraction of each phase in the mold were determined using the VOF model [26]. The DPM [27,28] was used to calculate the three-dimensional spatial distribution of air bubbles. The interactions between bubbles including the coalescence, breakup, and bounce were included to accurately simulate the interaction between the Eulerian phase and disperse bubbles. The momentum equation is defined as follows:
ρ t u i + ρ x j u i u j = p x i + x j μ + μ t u i x j + u j x i + F m o m , i
where u is the speed in m/s. ρ is the density in kg/m3. p is the pressure in Pa. μ is the viscosity in kg/(m·s). The subscript i and j represents the direction x, y, and z. The Smagorinsky–Lilly model [26] was applied in the current study and Equation (2) was used to calculate the turbulent viscosity μt. A user-defined function (UDF) was used to achieve Fmom, which was the source term caused by the interaction between the water phase and disperse air bubbles in kg/(m2·s2).
μ t = ρ L s 2 S ¯
where Ls is the mixing length and is defined as min κ d , C s V 1 / 3 . The κ is the Kármán constant and d is the closest wall’s distance in m. S is the rate-of-strain tensor for the resolved scale in s−1. V is the cell’s volume in m3. CS is the Smagorinsky constant with a value of 0.1.
The continuity equation is defined as
x i α q u q , i = 0
where the subscript q represents the q phase. αq represents the volume fraction of the q phase. ρq and uq represent the density in kg/m3 and the speed in m/s of the q phase, respectively. The density and viscosity of each computational cell are calculated using the mixing law.
The trajectory of each dispersed air bubble is solved using the DPM and the force balance on the bubble is defined as
d u b d t = F B + F D r a g + F L i f t + F P r e + F V M
where the FB, FDrag, FLift, FPre, FVM represent the interphase forces in m/s2. ub is the speed of bubbles in m/s. The expression of each interphase force is shown as follows:
F B = ρ b ρ l ρ b g
F D r a g = 3 4 μ C D Re ρ b d b 2 u l u b
F L i f t = C L ρ l ρ b u l u b × × u l
F V M = C V M ρ l ρ b u b u l d u b d t
where the CD, CL, and CVM represent the coefficients of drag, lift, and virtual mass, respectively. g is the acceleration of gravity in m/s2. Re is the Reynolds number and is defined as ρ l d b u l u b μ . The published studies have shown that these coefficients have a greater impact on interphase force and the specific value, as can be seen in reference [16]. db is the bubble diameter in m. The subscripts l and b represent the liquid phase and disperse bubble, respectively.
The collision of two different bubbles is calculated using the O’Rourke model [29]. The chance of collision is determined in this approach using the idea of a collision volume. The area of the π (r1 + r2)2 and the covered distance of the smaller bubble in one simulation time step is used to calculate the collision volume, as shown in Equation (9).
V c = π r 1 + r 2 2 u r e Δ t
where r1 and r2 represent the radius of two separate bubbles in m. ure is the relative speed of two separate air bubbles in m/s. ∆t is the time step employed in the current calculation in s.
Thus, based on the premise that the bubble has an equal chance of being anywhere in the cell, the definition of the probability of collision is the ratio of collision volume (Vc) and the volume of the computational cell (V) where the bubble is located, as shown in Equation (10).
P c = π r 1 + r 2 2 u r e Δ t V
For n small bubbles, the average number of collisions experienced by the large bubble in a time step is defined in Equation (11). In general, the average number in Equation (11) does not correspond to the actual number of collisions. In the O’Rourke model, Equation (12) demonstrates that a Poisson distribution describes the distribution of the probability of the collision number.
N ¯ = N n π r 1 + r 2 2 u r e Δ t V
P N = e N ¯ N ¯ N N !
where N is the collision number.
Two bubbles will coalescence or bounce once the collision occurred. Therefore, the collision’s outcome is determined by the collision parameter b, defined as r 1 + r 2 Y . Coalescence is the outcome of a collision when the real collision parameter b is smaller than the critical collision parameter bcrit. If it is not, there will be a bounce. The critical collision parameter bcrit is calculated as follows in the O’Rourke model:
b c r i t = r 1 + r 2 min 1.0 , 2.4 ( r 12 3 2.4 r 12 2 + 2.7 r 12 ) We
We = ρ u r e 2 d 1 d 2 σ
where Y is the random number with a value between 0 and 1. r12 is the ratio between the r1 and r2. σ is the surface tension in N/m. d1 and d2 are the diameter of two separate bubbles in m. The velocity of the new bubble after the coalescence is calculated in Equation (15). Also, the velocity of each bubble after the bounce is calculated in Equation (16) according to the basic conservation law.
u 1 = m 1 u 1 + m 2 u 2 m 1 + m 2 r 1 = 3 m 1 + m 2 4 π ρ b 1 3
u 1 = m 1 u 1 + m 2 u 2 m 1 + m 2 + m 2 u r e m 1 + m 2 b b c r i t r 1 + r 2 b c r i t u 2 = m 1 u 1 + m 2 u 2 m 1 + m 2 + m 1 u r e m 1 + m 2 b b c r i t r 1 + r 2 b c r i t
where m and u are the mass and speed of air bubbles in kg and m/s, respectively, and subscripts 1 and 2 represent bubble 1 and bubble 2, respectively.
The turbulent motion of the liquid water causes the air bubble diameter to have the highest critical value. The breakup will happen when the bubble diameter is greater than the critical value. Three breakup models including the Taylor model [30], Kelvin–Helmholtz model [31], and Stochastic model [32] were compared in the current study. In the Taylor model, the oscillation and distortion of a bubble are described by Equation (17). When the oscillation increases to a sufficient level, the bubble will fragment into a number of smaller bubbles.
d 2 X t 2 = C F ρ u r e 2 ρ b r C k σ ρ b r 3 X C d μ ρ b r 2 d X d t
where X represents the bubble’s displacement from its center position in m. CF, Ck, and Cd are the dimensionless constants and are set as 8.0, 5.0, and 1/3, respectively. The bubble is assumed to breakup when X > 0.5 r and the diameter of the newly formed bubbles is calculated according to the energy balance [33].
The Rayleigh–Taylor instabilities and the Kelvin–Helmholtz waves’ effects on the wave growth on the bubble surface are both considered in the Kelvin–Helmholtz model. Equations (18) and (19) calculate the frequency of the fastest growing wave and the matching wave number on the surface of the bubble, respectively.
Ω R T = 2 g ρ l ρ b 3 / 2 3 3 σ ρ l + ρ b
K R T = g ρ l ρ b 3 σ
Assuming that the bubble is broken up when the Rayleigh–Taylor waves have been growing for a time larger than 0.5/ΩRT, the radius of the newly formed bubble is calculated as follows:
r = π C R T K R T
where CRT is the constant with a value of 0.1.
The Stochastic model assumes that the bubble is subject to breakup when the breakup time is larger than a critical value. The breakup time is defined in Equation (21).
t b = B ρ b ρ l r u r e
where B is the constant number of the breakup model and is set as 1.73. The critical breakup time can be calculated using the critical radius as follows:
r c r i t = We crit σ ρ l u r e 2
Wecrit represents the critical Weber number, which, in the present investigation, is set as 6. The newly formed bubbles’ diameters are taken as a sample from the analytical Fokker–Planck equation for the probability distribution [32].

2.2. Calculation Domain

The current three-dimensional model was based on a slab CC water model published by Ren and Chen [16,34]. Figure 2 displays the computational domain and mesh distribution. About 0.58 million structured hexahedral meshes were used to divide up the computational domain, and the mesh close to the water–air interface was finetuned. The 510 mm × 50 mm CC mold and the bifurcated nozzle with a 40 mm immersion depth were included. The air phase with a 25 mm initial thickness was taken into consideration at the top. The air bubbles were injected at the SEN inlet. Table 1 provides a summary of the specific parameters that were employed in the present simulation.

2.3. Numerical Details

The fully coupled three-dimensional mathematical model was solved by combining commercial ANSYS FLUENT 2021 software and the UDF. The interaction between the liquid water and bubbles was achieved by considering of the interphase forces, and the effect of bubbles passing through the interface or breaking near the interface on the surface level was ignored. The interaction between the bubbles was achieved by the collision and breakup models. The changes in bubble diameter and velocity after the collision and breakup also have a direct impact on the velocity distribution of the liquid water and the surface level. At the SEN input and outflow of the domain, the velocity inlet and pressure outlet boundary conditions were used, respectively. The computation domain’s top surface was intended to be flat and insulated utilizing the free-slip condition. For other walls, the non-slip boundary condition was adopted. The air bubble diameter at the inlet of the SEN was 0.71 mm, determined through the water model [35]. The water–air interface was thought to be where bubbles were eliminated. The escape and reflection conditions were used at the outlet and other walls, respectively. For the pressure–velocity coupling solution, the PISO scheme was employed. The momentum equation was solved using the bounded central differencing scheme, and the equation of the volume fraction was calculated using the compressive method. All variables’ convergence thresholds were set to 10−4. The calculation time step of 0.0002 s was chosen to stabilize the calculation because the model was fully coupled and considered bubble coalescence, breakup, and bouncing. Each simulation was run for 130 s.

3. Comparison of Bubble Size and Spatial Distribution

3.1. Bubble Spatial Distribution

The typical bubble coalescence, breakdown, and bounce process, as predicted by the current model, is seen in Figure 3. The interaction between bubbles, particularly the probability of collision between bubbles, was enhanced by the turbulent transport of the liquid phase inside the SEN and the CC mold. With effective head-on collision, bubble diameter grows with time, and once it reaches the maximum critical diameter, it will disintegrate. The largest diameter of the bubbles that could remain stable at various points along the mold varied. The turbulence flow in the SEN was strong and the bubbles were prone to breakup. The deep section of the mold had weak turbulence flow and a bigger maximum critical bubble diameter. One of the factors contributing to the temporary asymmetric flow in the mold is the SEN’s higher frequency of bubble coalescence and breakup. The outcome of the collision tended to be a bounce if the collision was more oblique, which was the main bubble interaction in the deep part of the CC mold.
In Figure 4, various bubble-breakup models are used to depict the three-dimensional distribution of spatial of air bubbles at t = 102.0 s. The injecting flow rate was set as 90 mL/min, and the casting speed was set as 0.425 m/min. It can be seen that air bubbles were almost full of the entire mold in three bubble-breakup models. The lower recirculation flow carried the fine newly formed bubbles that broke off from the giant bubbles to the deep part of the CC mold. However, due to the different breakup mechanisms, the distribution of bubble diameter showed an obvious distinction with the Taylor model, K-H model, and Stochastic model.
The effect of the bubble-breakup model on the three-dimensional average bubble mass concentration was compared in Figure 5. The contour of the bubble concentration is shown on the left, and the bubble concentration ISO surface is shown on the right. The time-averaged results were calculated for 30 s after the two-phase flow attained a stable state. It can be seen that from the Taylor model, K-H model and Stochastic model, the bubbles’ degree of dispersion inside the CC mold increased gradually. Most of the bubbles were concentrated near the SEN with the Stochastic breakup model. The spatial distribution of bubbles with the Stochastic model in Figure 4c shows that there are many fine bubbles located in the deep part of the CC mold. However, the average bubble concentration in the mold deep region was almost zero in Figure 5c. This indicates that the number and size of fine air bubbles generated by bubble breakup with the Stochastic model were between the Taylor model and K-H model.
The effect of the bubble-breakup model on the average bubble concentration at 170 mm below the meniscus was compared in Figure 6. The bubble concentration with the Taylor model was significantly higher than that with the K-H model and Stochastic model. The bubble concentration was close to zero in the Stochastic model. The bubble concentration increased gradually from the mold center to the narrow face at the lower part of the mold. Because of the small buoyance of fine bubbles, they were more likely to move with the molten steel. Thus, the number of fine bubbles was higher near the narrow face at the deep part of the CC mold. The bubble concentration decreased from the SEN to the narrow face at the upper part of the mold due to the rise of the large bubbles. Figure 7 analyses the mass concentration of air bubbles along the casting direction quantitatively in more detail. From the Taylor model to the K-H model to the Stochastic model, the mass concentration of air bubbles in the vicinity of the SEN steadily rose. The average bubble concentration did, however, steadily decline close to the 1/4 mold width and narrow face. The biggest bubble concentration was found close to the SEN, and it was almost ten times closer to the narrow face.

3.2. Distribution of Bubble Diameter

The impact of the bubble-breakup model on the distribution of the bubble diameter is depicted in Figure 8. With the Taylor breakup model, the percentage of bubbles less than 0.25 mm reached 70%, whereas the percentage in the K-H model and the Stochastic model was between 15% and 20%. Another diameter peak under the Taylor model was 0.75 mm and reached 12.6%. For the K-H model and Stochastic model, the peak bubble diameter and number proportion was 0.45 mm and 33.3%, and 0.75 mm and 44.5%, respectively. Compared with the Taylor model and K-H model, the Stochastic model had the largest number proportion of bubbles with a 1–2 mm diameter. The number proportion of bubbles greater than 2 mm was very small.
The bubble’s average diameter with respect to time is shown in Figure 9. As can be observed, the bubble’s average diameter changed slightly with time under the current transient simulation and remained basically unchanged after the flow reached a stable state. The average bubble diameter in the entire mold was 0.253 mm, 0.440 mm, and 0.741 mm with the Taylor model, K-H model, and Stochastic model, respectively. The average bubble diameter obtained by the Stochastic model was in the best agreement with the measured values. The area between the narrow face and the SEN was divided into 10 equal zones, as indicated in Figure 10, and the bubble average diameter in each zone was computed. Due to their high buoyancy, the large size bubbles produced by the bubble coalescence climbed to the meniscus swiftly. Thus, from the SEN to the mold edge, the average bubble diameter dropped. With the Stochastic model, the bubble’s average diameter was 0.85 mm close to the SEN and 0.4 mm close to the narrow face.

3.3. Removal Position of Bubbles

In order to compare the bubble removal position at the meniscus, the position was recorded and outputted in real-time during the calculation. The removal position for differently sized bubbles and breakup models is shown in Figure 11. It is clear that from the removal position, the horizontal penetration distance of bubbles less than 1.7 mm can span the full width of the mold. Bubbles larger than 3.0 mm had a horizontal penetration distance that did not surpass half the mold’s width. With the K-H model, the removal position for bubbles larger than 3.0 mm was discovered close to the narrow face. This can be brought on by the big bubbles that coalesced in the deep part of the mold (Figure 3b).

4. Comparison of the Instantaneous Flow Field

The stochastic breakup model is employed to depict the instantaneous two-phase flow in Figure 12 at a specific time. The air injection rate was 90 mL/min, and the casting speed was 0.425 m/min. The flow field in the mold had an instantaneous asymmetrical distribution. At this time, there was a region with a large velocity near the meniscus on the right of the mold, resulting in fluctuations in the speed and liquid level of the meniscus significantly greater than those of the meniscus on the left of the mold. This asymmetric flow was caused by the collision and bounce of bubbles with the bottom of the SEN and the characteristics of turbulence. Moreover, the asymmetric distribution was exacerbated and even dominated by the coalescence, breakup, and bounce of bubbles, especially at the bottom of the SEN and the meniscus.
To quantitatively analyze how various bubble-breakup models affect the flow field, Figure 13 compares the predicted average speed and the measured velocity [35] at the center line of the meniscus. It can be seen that the speed predicted by the Stochastic breakup model is in the best agreement with the measured results, which is relative to the average bubble diameter distribution shown in Figure 9. The bubble-breakup model had a great influence on the speed distribution since the bubble-breakup model directly determined the bubble diameter and then influenced the distribution of the speed through the liquid–bubble interaction.
Figure 14 compares the distribution of average speed along the casting direction. The findings indicate that the area between the upper SEN outlet and the meniscus was where the bubble-breakup models had the greatest impact on the speed near the SEN. Figure 14a demonstrates how the speed within 12 mm of the meniscus gradually decreased using the Taylor, K-H, and Stochastic models. This is contrary to the bubble diameter distribution in Figure 10b. The effect of the bubble on the lowering of velocity is greater the larger the bubble’s diameter. It should be noted that the speed of liquid steel will increase when the bubble diameter rises to a certain point owing to the bubble’s floating.

5. Conclusions

This work suggests a validated numerical model that includes the LES model, VOF model, DPM, bubble-collision model, and bubble-breakup model. It was explained how the bubble-breakup model affected the spacing and size of dispersed bubbles and two-phase flow in the CC mold. The following is a succinct summary of the conclusions:
1.
The instantaneous asymmetrical distribution of bubbles and two-phase flow in a slab CC mold was successfully predicted. The two-way coupled DPM took into account the bubbles’ coalescence, breakup, and bounce.
2.
The distribution of the bubble diameter showed an obvious distinction between the different bubble-breakup models. The average bubble concentration near the SEN gradually increased from the Taylor model to K-H model and Stochastic model. However, the average bubble concentration was gradually increased from the 1/4 width of the mold to the narrow face.
3.
The predicted average bubble diameter in the entire mold was 0.253 mm, 0.440 mm, and 0.741 mm with the Taylor model, K-H model, and Stochastic model, respectively. The proportion of fine bubbles was overpredicted with the Taylor and K-H breakup model. The bubble-breakup model had a noticeable impact on the distribution of the speed due to the direct determination of the size distribution of bubbles.
4.
The Stochastic breakup model had the best agreement with the measured data when compared to the average bubble diameter and meniscus speed. Thus, the fully coupled LES model, VOF model, and DPM with the collision and breakup model is recommended to correctly calculate the two-phase flow and distribution of multi-size bubbles during the CC process.

Author Contributions

Conceptualization, Y.T., L.X., S.Q. and R.Z.; methodology, Y.T. and S.Q.; software, validation and writing—original draft preparation, Y.T.; formal analysis and visualization, Y.T.; investigation, Y.T.; resources and supervision, L.X., S.Q. and R.Z.; data curation, S.Q.; writing—review and editing, L.X. and S.Q.; project administration and funding acquisition, S.Q. and R.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 52074024).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Authors are grateful to the National Natural Science Foundation of China for support of this research, Grant No. 52074024. Thanks are extended to Lei Lui and Jiaqi Liu for the useful technical discussions at Central Iron and Steel Research Institute (CISRI), Beijing, China.

Conflicts of Interest

The authors declare no conflict of interests.

References

  1. Thomas, B.G.; Zhang, L. Mathematical modeling of fluid flow in continuous casting. ISIJ Int. 2001, 41, 1181–1193. [Google Scholar] [CrossRef]
  2. Chen, W.; Zhang, L.; Ren, Q.; Ren, Y.; Yang, W. Large eddy simulation on four-phase flow and slag entrainment in the slab continuous casting mold. Met. Mater. Trans. B 2022, 53, 1446–1461. [Google Scholar] [CrossRef]
  3. Zhang, L.; Thomas, B.G. State of the art in evaluation and control of steel cleanliness. ISIJ Int. 2003, 43, 271–291. [Google Scholar] [CrossRef] [Green Version]
  4. Thomas, B.G. Review on modeling and simulation of continuous casting. Steel Res. Int. 2018, 89, 1700312. [Google Scholar] [CrossRef]
  5. Liu, Z.-Q.; Qi, F.-S.; Li, B.-K.; Jiang, M.-F. Vortex flow pattern in a slab continuous casting mold with argon gas injection. J. Iron Steel Res. Int. 2014, 21, 1081–1089. [Google Scholar] [CrossRef]
  6. Bai, H.; Thomas, B.G. Turbulent flow of liquid steel and argon bubbles in slide-gate tundish nozzles: Part I. model development and validation. Metall. Mater. Trans. B 2001, 32, 253–267. [Google Scholar] [CrossRef]
  7. Sánchez-Pérez, R.; García-Demedices, L.; Ramos, J.P.; Díaz-Cruz, M.; Morales, R.D. Dynamics of coupled and uncoupled two-phase flows in a slab mold. Metall. Mater. Trans. B 2004, 35, 85–99. [Google Scholar] [CrossRef]
  8. Singh, V.; Dash, S.K.; Sunitha, J.S.; Ajmani, S.K.; Das, A.K. Experimental simulation and mathematical modeling of air bubble movement in slab caster mold. ISIJ Int. 2006, 46, 210–218. [Google Scholar] [CrossRef] [Green Version]
  9. Li, L.-M.; Liu, Z.-Q.; Li, B.-K. Modelling of bubble aggregation, breakage and transport in slab continuous casting mold. J. Iron Steel Res. Int. 2015, 22, 30–35. [Google Scholar] [CrossRef]
  10. Liu, Z.; Qi, F.; Li, B.; Cheung, S. Modeling of bubble behaviors and size distribution in a slab continuous casting mold. Int. J. Multiph. Flow 2016, 79, 190–201. [Google Scholar] [CrossRef]
  11. Liu, Z.; Qi, F.; Li, B.; Jiang, M. Multiple size group modeling of polydispersed bubbly flow in the mold: An analysis of turbulence and interfacial force models. Metall. Mater. Trans. B 2015, 46, 933–952. [Google Scholar] [CrossRef]
  12. Yu, H.; Zhu, M. Numerical simulation of the effects of electromagnetic brake and argon gas injection on the three-dimensional multiphase flow and heat transfer in slab continuous casting mold. ISIJ Int. 2008, 48, 584–591. [Google Scholar] [CrossRef] [Green Version]
  13. Jin, K.; Thomas, B.G.; Ruan, X. Modeling and measurements of multiphase flow and bubble entrapment in steel continuous casting. Metall. Mater. Trans. B 2016, 47, 548–565. [Google Scholar] [CrossRef]
  14. Cho, S.-M.; Thomas, B.G.; Kim, S.-H. Bubble behavior and size distributions in stopper-rod nozzle and mold during continuous casting of steel slabs. ISIJ Int. 2018, 58, 1443–1452. [Google Scholar] [CrossRef] [Green Version]
  15. Chen, W.; Zhang, L. Effects of interphase forces on multiphase flow and bubble distribution in continuous casting strands. Metall. Mater. Trans. B 2021, 52, 528–547. [Google Scholar] [CrossRef]
  16. Chen, W.; Ren, Y.; Zhang, L. Large eddy simulation on the two-phase flow in a water model of continuous casting strand with gas injection. Steel Res. Int. 2019, 90, 1800287. [Google Scholar] [CrossRef]
  17. Lopez, P.E.R.; Jalali, P.N.; Björkvall, J.; Sjöström, U.; Nilsson, C. Recent developments of a numerical model for continuous casting of steel: Model theory, setup and comparison to physical modelling with liquid metal. ISIJ Int. 2014, 54, 342–350. [Google Scholar] [CrossRef] [Green Version]
  18. Chen, W.; Ren, Y.; Zhang, L.; Scheller, P.R. Numerical simulation of steel and argon gas two-phase flow in continuous casting using LES + VOF + DPM model. JOM 2019, 71, 1158–1168. [Google Scholar] [CrossRef]
  19. Chen, W.; Zhang, L.; Wang, Y.; Ren, Y.; Ren, Q.; Yang, W. Prediction on the three-dimensional spatial distribution of the number density of inclusions on the entire cross section of a steel continuous casting slab. Int. J. Heat Mass Transf. 2022, 190, 122789. [Google Scholar] [CrossRef]
  20. Wang, Y.; Fang, Q.; Zhang, H.; Zhou, J.; Liu, C.; Ni, H. Effect of argon blowing rate on multiphase flow and initial solidification in a slab mold. Metall. Mater. Trans. B 2020, 51, 1088–1100. [Google Scholar] [CrossRef]
  21. Yang, W.; Luo, Z.; Gu, Y.; Liu, Z.; Zou, Z. Simulation of bubbles behavior in steel continuous casting mold using an Euler-Lagrange framework with modified bubble coalescence and breakup models. Powder Technol. 2020, 361, 769–781. [Google Scholar] [CrossRef]
  22. Zhang, T.; Luo, Z.; Liu, C.; Zhou, H.; Zou, Z. A mathematical model considering the interaction of bubbles in continuous casting mold of steel. Powder Technol. 2015, 273, 154–164. [Google Scholar] [CrossRef]
  23. Chaudhary, R.; Thomas, B.G.; Vanka, S.P. Effect of electromagnetic ruler braking (embr) on transient turbulent flow in continuous slab casting using large eddy simulations. Metall. Mater. Trans. B 2012, 43, 532–553. [Google Scholar] [CrossRef]
  24. Real-Ramirez, C.A.; Carvajal-Mariscal, I.; Sanchez-Silva, F.; Cervantes-de-la-Torre, F.; Diaz-Montes, J.; Gonzalez-Trejo, J. Three-dimensional flow behavior inside the submerged entry nozzle. Metall. Mater. Trans. B 2018, 49, 1644–1657. [Google Scholar] [CrossRef]
  25. Hirt, C.W.; Nichols, B.D. Volume of fluid (VOF) method for the dynamics of free boundaries. J. Comput. Phys. 1981, 39, 201–225. [Google Scholar] [CrossRef]
  26. Smagorinsky, J. General circulation experiments with the primitive equations: I. The basic experiment. Mon. Weather. Rev. 1963, 91, 99–164. [Google Scholar] [CrossRef]
  27. Djebbar, R.; Roustan, M.; Line, A. Numerical computations of turbulent gas-liquid dispersion in mechanically agitated vessels. Chem. Eng. Res. Des. 1996, 74, 492–498. [Google Scholar]
  28. Tomiyama, A.; Tamai, H.; Zun, I.; Hosokawa, S. Transverse migration of single bubbles in simple shear flows. Chem. Eng. Sci. 2002, 57, 1849–1858. [Google Scholar] [CrossRef]
  29. O’Rourke, P.J. Collective Drop Effects on Vaporizing Liquid Sprays. Ph.D. Thesis, Princeton University, Princeton, NJ, USA, 1981. [Google Scholar]
  30. Taylor, G.I. The shape and acceleration of a drop in a high speed air stream. Sci. Pap. GI Taylor 1963, 3, 457–464. [Google Scholar]
  31. Reitz, R.D.; Beale, J.C. Modeling spray atomization with the Kelvin-Helmholtz/Rayleigh-Taylor hybrid model. At. Sprays 1999, 9, 623–650. [Google Scholar] [CrossRef]
  32. Apte, S.; Gorokhovski, M.; Moin, P. LES of atomizing spray with stochastic modeling of secondary breakup. Int. J. Multiph. Flow 2003, 29, 1503–1522. [Google Scholar] [CrossRef]
  33. O’Rourke, P.J. The TAB Method for Numerical Calculation of Spray Drop Breakup; SAE Technical Paper; SAE: Warrendale, PA, USA, 1987; No. 872089-1987. [Google Scholar]
  34. Ren, L.; Ren, Y.; Zhang, L.; Yang, J. Investigation on fluid flow inside a continuous slab casting mold using particle image velocimetry. Steel Res. Int. 2019, 90, 1900209. [Google Scholar] [CrossRef]
  35. Tian, Y.; Shi, P.; Xu, L.; Qiu, S.; Zhu, R. Numerical Modeling of Transient Two-Phase Flow and the Coalescence and Breakup of Bubbles in a Continuous Casting Mold. Materials 2022, 15, 2810. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic of multiphase flow in the mold during CC process.
Figure 1. Schematic of multiphase flow in the mold during CC process.
Materials 16 04666 g001
Figure 2. Computational domain and mesh distribution.
Figure 2. Computational domain and mesh distribution.
Materials 16 04666 g002
Figure 3. The procedure of the coalescence, breakup, and bounce of air bubbles in the CC mold.
Figure 3. The procedure of the coalescence, breakup, and bounce of air bubbles in the CC mold.
Materials 16 04666 g003
Figure 4. Three-dimensional distribution of spatial of air bubbles at t = 102.0 s: (a) Taylor model, (b) K-H model, and (c) Stochastic model.
Figure 4. Three-dimensional distribution of spatial of air bubbles at t = 102.0 s: (a) Taylor model, (b) K-H model, and (c) Stochastic model.
Materials 16 04666 g004
Figure 5. Comparison of bubbles’ average mass concentration with different bubble-breakup models: (a) Taylor model, (b) K-H model, and (c) Stochastic model.
Figure 5. Comparison of bubbles’ average mass concentration with different bubble-breakup models: (a) Taylor model, (b) K-H model, and (c) Stochastic model.
Materials 16 04666 g005
Figure 6. Effect of the breakup model on bubble average concentration at 170 mm below the meniscus.
Figure 6. Effect of the breakup model on bubble average concentration at 170 mm below the meniscus.
Materials 16 04666 g006
Figure 7. Mass concentration of air bubbles along the casting direction: (a) near the SEN, (b) 1/4 width of the CC mold, and (c) near the narrow face.
Figure 7. Mass concentration of air bubbles along the casting direction: (a) near the SEN, (b) 1/4 width of the CC mold, and (c) near the narrow face.
Materials 16 04666 g007aMaterials 16 04666 g007b
Figure 8. Effect of the bubble-breakup model on the diameter of bubbles in the CC mold.
Figure 8. Effect of the bubble-breakup model on the diameter of bubbles in the CC mold.
Materials 16 04666 g008
Figure 9. Variation of the diameter with time.
Figure 9. Variation of the diameter with time.
Materials 16 04666 g009
Figure 10. Effect of the bubble-breakup model on the average bubble diameter at different regions: (a) statistical interval, (b) distribution of the average bubble diameter.
Figure 10. Effect of the bubble-breakup model on the average bubble diameter at different regions: (a) statistical interval, (b) distribution of the average bubble diameter.
Materials 16 04666 g010
Figure 11. Effect of the bubble-breakup model on the removal position at the meniscus: (a) Taylor model, (b) K-H model, and (c) Stochastic model.
Figure 11. Effect of the bubble-breakup model on the removal position at the meniscus: (a) Taylor model, (b) K-H model, and (c) Stochastic model.
Materials 16 04666 g011
Figure 12. Distribution of the instantaneous flow filed at a certain moment using the Stochastic breakup model.
Figure 12. Distribution of the instantaneous flow filed at a certain moment using the Stochastic breakup model.
Materials 16 04666 g012
Figure 13. Comparison of the predicted and measured [35] average speed along the center line of the meniscus.
Figure 13. Comparison of the predicted and measured [35] average speed along the center line of the meniscus.
Materials 16 04666 g013
Figure 14. Distribution of the average speed along the casting direction: (a) 5 mm from the SEN, (b) 1/4 width of the CC mold, and (c) 5 mm from the narrow face.
Figure 14. Distribution of the average speed along the casting direction: (a) 5 mm from the SEN, (b) 1/4 width of the CC mold, and (c) 5 mm from the narrow face.
Materials 16 04666 g014
Table 1. Parameters of the simulation.
Table 1. Parameters of the simulation.
ParameterValue
Total length562 mm
Cross size510 mm × 50 mm
SEN outlet angle
SEN immersion depth40 mm
Casting speed0.425 m/min
Air flow rate90 mL/min
Density of water1000 kg/m3
Viscosity of water0.001 kg·m−1·s−1
Density of air1.225 kg/m3
Viscosity of air1.789 × 10−5 kg·m−1·s−1
Surface tension0.07197 N/m
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tian, Y.; Xu, L.; Qiu, S.; Zhu, R. Simulation of Spatial Distribution of Multi-Size Bubbles in a Slab Continuous-Casting Mold Water Model. Materials 2023, 16, 4666. https://doi.org/10.3390/ma16134666

AMA Style

Tian Y, Xu L, Qiu S, Zhu R. Simulation of Spatial Distribution of Multi-Size Bubbles in a Slab Continuous-Casting Mold Water Model. Materials. 2023; 16(13):4666. https://doi.org/10.3390/ma16134666

Chicago/Turabian Style

Tian, Yushi, Lijun Xu, Shengtao Qiu, and Rong Zhu. 2023. "Simulation of Spatial Distribution of Multi-Size Bubbles in a Slab Continuous-Casting Mold Water Model" Materials 16, no. 13: 4666. https://doi.org/10.3390/ma16134666

APA Style

Tian, Y., Xu, L., Qiu, S., & Zhu, R. (2023). Simulation of Spatial Distribution of Multi-Size Bubbles in a Slab Continuous-Casting Mold Water Model. Materials, 16(13), 4666. https://doi.org/10.3390/ma16134666

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop