1. Introduction: A Review of Previous Analytical/Numerical Models for Slag Infiltration
It is known that molten slag offers lubrication through infiltration into the mold–strand gap during continuous casting [
1,
2]. The infiltrated liquid slag forms a thin film between solidified shell and solid slag and functions as a lubricant in between to avoid solid friction. The mechanism of slag infiltration has been investigated by numerous researchers due to its importance on heat transfer [
3], mold friction [
4], oscillation mark formation [
5,
6], and so on. However, it has not reached a consensus yet in several issues although the mechanism has been studied from diverse points of view including physical modeling [
7,
8,
9,
10], mathematical modeling [
4,
6,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20], and plant trials [
2,
9,
21,
22,
23,
24].
The first argument is when the slag infiltration takes place during the mold oscillation cycle. The timing of slag infiltration is an important factor for designing mold oscillation settings (stroke, frequency, and wave shape) as higher consumption can be achieved once the timing is known by extending the period of positive consumption through setting modification (e.g., non-sinusoidal mold oscillation) [
21]. Typically, slag infiltration is discussed with concepts called Positive Strip Time (PST) and Negative Strip Time (NST): NST is the period when the mold wall descends faster than the solidified shell, and the remaining period in the cycle becomes PST [
2].
Most of the modeling results show that slag consumption mainly occurs during NST [
9,
12,
14,
15,
16,
18,
19,
25]. It is quite clear from a fluid mechanics perspective that the slag flow is driven by solidified shell and mold wall in the narrow gap. Therefore, slag consumption is expected to reach its peak when both walls are dragging the slag downward during NST. Also, a pumping effect from slag rim is expected to enhance the consumption during NST [
25]. Numerous modeling approaches have been attempted mathematically from 1D analytical solutions [
9,
11,
12,
19] to 2D models [
15,
16,
18] and full 3D numerical model [
14] and showed corresponding results. Also, results from physical modeling using oil [
7,
8] agree with this conclusion.
However, statistical analysis of plant data seems to disagree with this hypothesis by showing a stronger correlation between consumption and PST rather than NST [
2,
9,
21,
23]. Based on this proportionality between PST and consumption, a non-sinusoidal oscillation mode has been proposed and successfully applied in high-speed continuous casting [
21]. Some lab-scale experiments using liquid metal [
10] and 2D numerical models [
6,
17] have been able to reproduce a similar trend that slag consumption during PST is greater than NST, while the study mainly focused on describing the phenomenon. Thus, a reasonable explanation for the different conclusions between modeling and reality is missing to date. Also, no clear answer has been given for the mechanism of higher slag consumption with non-sinusoidal mold oscillation so far despite its success in high-speed casting production.
Another controversy is regarding the slope of slag film channel. According to the lubrication theory, the liquid slag film requires a slope to open a gap between solidified shell and solid slag [
26]. Kajitani et al. claimed that the gap should diverge into the casting direction to reproduce the same behavior of real casters from their cold experiment (i.e., slag consumption decreases with higher casting speed or higher slag viscosity in real casters) [
7,
8]. Physically, a diverging mold–strand gap is possible to occur when a mold taper is not enough to compensate for the thermal shrinkage of solidified shell [
14]. However, mathematical models considering heat transfer show that the liquid film thickness is more likely to decrease gradually from the meniscus to the end of liquid lubrication length [
4,
6]. Also, following numerical works showed that the trends from real casters can be achieved with a converging liquid film as well [
18,
20]. The microstructure of slag film samples from real casters seems to back up the converging film based on the thickness of glassy slag layer [
1]. Recently, the author has developed a 3D numerical model that can estimate the size of mold–strand gap based on thermal shrinkage, mold taper, and ferrostatic pressure [
14]. Results from the model show that the size of mold–strand gap is determined by how much the thermal shrinkage is compensated by mold taper (defined as Contact Index). In other words, the gap profile varies locally depending on the steel grade (e.g., thermal expansion coefficient) and mold design as well as heat transfer. However, a detailed study of the slope effect on slag infiltration has not been implemented as the gap profile and liquid film thickness are results of the model calculation rather than input parameters.
In this article, the author aims to investigate these controversial issues from a fluid mechanics perspective. Although several advanced numerical models have been proposed up to date, the complexity of the models obscures clarification of these issues. Therefore, a simple but reliable 1D analytical model for the slag flow in mold–strand gap is developed using the superposition of analytical solutions for wall-driven flows (Couette flow and Stokes’ second problem typed flow) and pressure-driven flow (Poiseuille flow) [
27]. This model can handle transient effects of the flow during arbitrary mold oscillation including non-sinusoidal modes that have not been considered in previous analytical models. With this new feature, the timing of slag infiltration, slope effect of slag film, and mechanism of slag infiltration with non-sinusoidal mold oscillation are studied.
2. Mechanism of Slag Infiltration
Among the several roles of casting power during continuous casting, providing lubrication between mold and solidified steel shell is crucial since it affects the casting process and surface quality directly. The formation of slag film in the process is described in
Figure 1a. The liquid and solid slag supplied from the slag pool and powder bed are infiltrated into the gap between mold and shell while the shell is withdrawn into the casting direction. It is known that mold oscillation is the key to keep the stable infiltration of slag film during the withdrawal.
Figure 1b shows more details of the slag film in the gap. As shown in the figure, the infiltrated liquid slag film is the only liquid phase which can offer lubrication between solid phases. Therefore, deep penetration of the liquid slag film as well as higher slag consumption is preferred in the casting process.
The prediction of slag infiltration requires two components, slag film thickness and velocity profile of the slag film (
Figure 1b). As a multidisciplinary approach is required for the estimation of slag film thickness (or mold–strand gap), analytical models typically require the thickness as an external input. Several methods have been proposed to calculate the film thickness separately such as using a vibration theory [
11,
12], Bikerman equation [
19], coupling with heat transfer [
4,
13], and so on. Despite analytical models cannot offer the whole picture of slag infiltration, it gives a useful insight for slag film pressure. In addition to the role as a lubricant to avoid solid–solid contact, the liquid slag film is expected to generate positive pressure against compression between shell and mold (resulting from the thermo-mechanical interaction between thermal shrinkage, mold taper, and ferrostatic pressure) to open the gap [
14] as shown in
Figure 1b. Physically, the liquid film needs to be squeezed into the flow direction to generate positive pressure against compression. This implies that the profile of liquid film thickness is likely to have a slope rather than being parallel into the casting direction, but it does not indicate the sign of inclination because the mold wall oscillates vertically. In other words, mold oscillation allows continuous casting to adapt to both scenarios that the liquid film channel (1) diverges into the casting direction [
7,
8,
11,
12,
13,
19] or (2) converges into the casting direction [
4,
14,
18,
20].
In the case of converging liquid film into the casting direction, positive pressure is generated mostly during the Negative Strip Time (NST) or downstroke of mold oscillation more precisely according to the lubrication theory (
Figure 2a): The slag flow dragged by both walls (shell and mold) into the narrowing gap builds up positive pressure and it widens the gap. This is ideal for the slag infiltration because the liquid slag flowing into the casting direction reaches its peak when the gap widens during the NST. On the other hand, negative pressure develops during the upstroke
(
Figure 2b) when the slag flow rate
becomes negative (i.e.,
). Naturally, this leads to the closing of mold–strand gap.
Compared to the converging film, positive pressure develops during the Positive Strip Time (PST) or upstroke of mold oscillation in the scenario of diverging liquid film (
Figure 3b): The slag flow is squeezed in the channel when the slag flow rate
becomes negative. For the negative flow rate
, the flow driven upward by the mold wall needs to be greater than the flow driven downward by the solidified shell. Therefore, the mold wall velocity must be higher than the casting speed to create positive pressure in this scenario. The diverging liquid film is not as favorable as the converging liquid film in terms of infiltration because the liquid film shrinks during the NST when the positive slag flow (
> 0) is maximized due to the negative pressure developed in the gap (
Figure 3a). When the flow is squeezed by approaching two walls in both scenarios (e.g., downstroke with converging film in
Figure 2a or upstroke with diverging film in
Figure 3b), the positive peak pressure in the film channel pushes the slag away to both ends of the channel that correspond to the meniscus and end of the liquid film channel. Pressure gradients formed by the positive pressure generate a flow against the wall-driven-flow in the wider side of gap (i.e., Adverse Pressure Gradient), but promote the flow in the narrower side (i.e., Favorable Pressure Gradient). On the other hand, the generated negative pressure during the upstroke with converging film (
Figure 2b) or the downstroke with diverging film (
Figure 3a) sucks the slag from both ends. Accordingly, the opposite direction of pressure-driven flow is generated compared to the cases with positive pressure.
In the slag film channel, there are two different kinds of driving forces that infiltrate the liquid film into the gap: (1) Shell motion (i.e., casting speed) and mold oscillation (i.e., wave shape, stroke, and frequency), (2) gravity and pressure gradient in the gap. From the first driving forces, a wall-driven flow is expected: The steady motion of solidified shell generates a flow such as a Couette flow while the oscillating motion of mold wall drives a Stokes’-2nd-problem-typed flow. When a constant viscosity is assumed, a linear velocity profile develops between the oscillating mold velocity and the steady shell velocity if the viscous effect is dominant, otherwise, some non-linear profile appears depending on the film thickness, kinematic viscosity of slag, and mold frequency. This transient effect results in a transition of the velocity profile due to the oscillatory motion of mold wall (
Figure 4).
In addition to the wall-driven flow, the gravity and pressure gradient generate a pressure-driven flow (so-called, Poiseuille flow). This type of flow typically shows a parabolic velocity profile with its peak in the middle such as the velocity profiles shown in
Figure 2 and
Figure 3. As explained above with the figures, the pressure gradients are developed in the gap by the continuity effect in cases that the film channel has a slope (e.g., diverging or converging film): The film pressure varies locally to satisfy the constant flow rate
Q in the film so that the flow driven by pressure gradients makes up for the change of thickness. Therefore, the magnitude of pressure gradient is determined by how much the film thickness decreases or increases between inlet and outlet (i.e.,
in
Figure 2 and
Figure 3).
It is crucial to consider these wall-driven and pressure-driven flows together to predict the slag behavior accurately. For example, from the fact that the thickness of liquid slag film decreases near the meniscus (i.e., converging film), positive pressure develops during NST but negative pressure during PST. Therefore, some counter-intuitive behavior such as positive and higher consumption during PST could be explained by periodic suction (during PST) and overflow (during NST) if the magnitude of pressure-driven flow is comparable to the wall-driven flow.
In the current study, the analytical model focuses on the transient velocity and pressure of slag film with a certain slope. The slope is calculated from a force balance between liquid film pressure and compression in the mold–strand gap. Therefore, the model can predict both scenarios (converging and diverging film channel) based on the compression applied to the mold–strand gap. The ferrostatic pressure is used for the compression under the assumption that the ferrostatic pressure fully transmits to the solidified shell, which is a possible scenario in the real process [
14]. The slag film thickness is calculated from measurement data of casting powder consumption. Therefore, the film thickness has a spatial variation based on the calculated slope but static in time.
3. Description of the Model
As the slag flow in mold–strand gap can be described by a combination of Couette flow
, Stoke’s second problem typed flow
, and Poiseuille flow
, the velocity profile of slag film
is obtained by the superposition of analytical solutions as follows:
where,
Figure 5 describes the schematic of simulation domain with its coordinate system.
The superposition principle (Equation (1)) is valid since the governing equation is linear mathematically. The transient displacement
and velocity of mold wall
are obtained from the equations below [
28] based on oscillation settings such as amplitude
and angular frequency
:
The coefficients and are fitting parameters for non-sinusoidal mold oscillation. The equations converge to sinusoidal oscillation when .
As the analytical solution
is for sinusoidal oscillation (
), Fourier expansion is implemented to approximate a non-sinusoidal wave to a sum of multiple sinusoidal waves in cases of non-sinusoidal mold oscillation:
In this study, ten cosine and ten sine functions are used for the approximation. An advantage of using the Fourier expansion is that any kind of periodic function can be approximated into a sum of sinusoidal functions. In real scenarios with an industrial mold oscillator, a conversion of a sampled mold velocity from the oscillator into a mathematical form
using Equation (11) is required in cases that a periodic function for the mold oscillation is not available. After the approximation through Fourier expansion, the analytical solution for
with each sinusoidal oscillation is superposed (i.e.,
) to obtain the flow driven by non-sinusoidal mold oscillation. The pressure gradient
in
is obtained separately from the equation derived from the Reynolds equation below [
11]:
Here, a static film thickness h which has a spatial variation in
x (casting direction) is assumed as follows:
where,
is the slope of the film,
is the lubrication length, and
is the slag film thickness at the middle of lubrication length (
x =
L/2). The pressure gradient term includes the continuity effect from slope and pressure difference between inlet
and outlet
of the slag film channel. In cases of no slope (i.e., parallel slag film and h=const), the pressure gradient converges simply to
. The inlet pressure
is calculated based on the depth of slag pool, and the outlet pressure
is assumed to be zero. Equations (1)–(8) are solved with boundary conditions below:
At
y = 0 (contact with solid slag)
At
y = 1 (contact with solidified shell)
The static motion of solidified shell with casting speed drives a linear velocity profile across the film as shown in Equation (2), while the gravity and pressure gradient produce a parabolic flow with its peak at the center of film thickness as described in Equation (3). The oscillatory motion of mold wall generates an unsteady flow which profile can vary between linear and non-linear depending on the relation between momentum transfer, film thickness, and frequency. The final flow obtained from the summation of three flows (Equation (1)) representing the viscous effect , pressure gradient effect and transient effect result in diverse velocity profiles depending on the contribution of each flow.
The flow chart shown in
Figure 6 explains the algorithm of analytical model. Firstly, the oscillating mold wall velocity is obtained based on an industrial oscillation setting. In cases of non-sinusoidal mold oscillation, the fitting parameters (
and
in Equations (10) and (11)) are tuned to get a certain modification ratio
or targeted PST. After that, the mold oscillation velocity described by Equation (11) is approximated to a summation of sinusoidal waves (Equation (12)) through Fourier expansion. Next, based on the externally given slag film thickness
h and an initially guessed slope
s, the steady flow (
is calculated from the analytical solution (Equation (2)). Then, the unsteady flows
) are calculated with a certain time step. The pressure gradient term is obtained based on the transient mold wall velocity
as well as slope of the channel (Equation (13)). At every time step, all the flow solutions are summed to obtain the final velocity profile (
). Transient slag consumption and force from film pressure are obtained based on the transient pressure and velocity profiles. After that, a time-averaged force from film pressure is compared to compression from ferrostatic pressure (
, and the loop is iterated to find a proper slope
until the time-averaged force balances the compression.
4. Validation of the Model
The analytical model is tested with measurement data from plant trials in Shin et al. [
2]. Predicted slag consumption is compared to the measured powder consumption for the validation. The analytical model requires a liquid slag film thickness
as an input. Several assumptions are made to estimate the thickness from the measurement data. Firstly, the solid slag film is assumed to be attached to the mold wall. In other words, liquid slag only contributes to slag consumption. Secondly, the net slag consumption per cycle from
is assumed to be zero from the fact that it drives a periodic flow. Under these assumptions, the slag consumption can be calculated by the equations below from the analytical solutions for
(Equation (2)) and
(Equation (3)):
where the
is the mass flow rate of slag film
, the
is the slag consumption per steel surface area
. There are three terms contributing to the slag consumption
:
term from steady shell motion
,
term from gravity,
term from pressure gradient. One thing worth mentioning here is the relationship between consumption and casting speed
. As shown in Equation (19), the slag consumption
seems increasing with casting speed
due to the contribution of
(
term) but decreasing with viscosity by the other terms. On the other hand, the relation in inverse-proportion between slag consumption and casting speed or viscosity appears in the definition of
by the terms associated with gravity and pressure gradient. However, the prediction of slag consumption is not as explicit as it appears in the equations since the film thickness
h is a function of casting speed
. Therefore, it is not a complete solution unless the relationship between casting speed
and film thickness
h is given by another model, although the analytical solution can predict the same consumption trend to real casters such as less consumption
with higher casting speed
or viscosity
to some extent. The formula for slag consumption (Equation (20)) can be simplified further to
when the slag film is very thin (
h << 1). This means that the contribution of pressure-driven flows (gravity and pressure gradient) diminishes as the film becomes thinner and the consumption is determined by solidified shell motion (
term from
). From the simplified formula
, the liquid slag film thickness can be estimated by
This simplified equation is used to calculate the liquid film thickness
(in Equation (16)) from the casting powder consumption in Shin et al. [
2].
Table 1 shows the operating conditions used in the plant trials [
2] and measured powder consumptions. Calculations for 20 different cases including three oscillation modes (sinusoidal and non-sinusoidal oscillation with two modification ratios,
) are compared to the measurement data. The lubrication length of liquid slag film is assumed to
L = 0.1 m. The same slag properties used in Shin et al. [
2] (slag density
, slag viscosity
) are applied in the calculations.
Figure 7 visualizes mold displacement and velocity for Case A2 (sinusoidal), A3 (
and A4 (
calculated by Equations (10) and (11). It shows a clear trend from
Figure 7b that the NST becomes narrower but the descending velocity becomes faster as the modification ratio increases.
With the calculated mold wall velocities (
Figure 7b) as a boundary condition, slag flows in mold–strand gap for Case A2 and A4 are compared in
Figure 8. The estimated liquid slag film thicknesses from Equation (21) for A2 and A4 are 173
and 189
respectively. Here,
y = 0 is the mold wall or solid slag layer fully attached to the mold wall, and
y = 1 is the solidified shell. The calculated velocity profiles at the center of film channel (
x =
L/2) are almost linear during the oscillation cycle. This means that the influences of gravity and pressure gradient are imperceptible in the thin film thickness used in the calculations as assumed in the derivation of Equation (21). Namely, the formula to calculate liquid film thickness from slag consumption
is valid with thickness
h < ~200 μm, which falls into a typical range of liquid film thickness observed in continuous casting [
22]. Also, any transient effects such as the transition of velocity profile (
Figure 4) are not clearly observed in the results and the entire velocity profile responds quickly to the mold oscillation. Thus, the viscous effect dominates the overall slag flow due to its thin film thickness.
Figure 9 shows the difference of velocity profiles between top (
x = 0) and bottom (
x =
L) of the film channel during the oscillation cycle in Case A2.
Compared to the velocity profile obtained at the center (
x = L/2), they show parabolic profiles by the continuity effect as described in
Figure 2 and
Figure 3. Because of the converging film (
s = −4.75 × 10
−4) which slope is calculated from the force balance through iterations (
Figure 6), positive pressure develops during the downstroke period
and it pushes the slag flow to both ends (e.g.,
t = 3
T/8 or
t = 4
T/8). This continuity or pressure gradient effect promotes the flow at the bottom of slag channel (
x =
L) but interrupts the flow at the top (
x = 0) by bending the velocity profile to parabolic shapes with opposite directions depending on the location. The same mechanism is applicable for explaining the flow during the upstroke period
with negative pressure (e.g.,
t = 0).
The transient slag consumption and force from slag pressure for Case A2, A3, and A4 are shown in
Figure 10. Both profiles show a similar trend with regards to the modification ratio: The distribution of both curves becomes sharper as the modification ratio increases, which corresponds to the mold wall velocity in
Figure 7. Also, both reach their peak during the NST but the lowest negative consumption and force during the PST are obtained with the highest modification ratio (Case A4,
). The calculated net consumptions for A2, A3, and A4 are 0.2285, 0.2219, and 0.2486
, which agree well to the measurement data in
Table 1. On the contrary, it shows a small discrepancy in terms of the timing of peak. The slag consumption shows a delay of its peak as the modification ratio increases, while the force responds to the mold velocity precisely so that the maximum force is obtained at
t =
T/2 when the mold descending velocity is its peak. It is suspected that the transient effect discussed in
Figure 4 is responsible for the small delay in consumption. More details of this effect are discussed in
Section 5.
Figure 11 shows the estimated liquid film thickness and slope of the film from the analytical model and comparison of the consumption to the measurement data. The results are plotted with the modification ratio to visualize the effect of non-sinusoidal mold oscillation with Case A1–A10.
The results show a corresponding trend that the liquid film thickness, slope of the film, and slag consumption increase with the modification ratio. All the estimated slopes for slag film are negative to generate a positive force against compression. From the relation that the slag film pressure decreases with film thickness but increases with slope magnitude (Equations (13)–(15)), a steeper slope is obtained with a thicker slag film to generate the same magnitude of force against compression. The calculated consumption also increases with the modification ratio and agrees well with the measurement data.
Figure 12 shows the results for Case B1–B11 to visualize the effect of casting speed by fixing the modification ratio to 24%.
As expected, the calculated slag consumption decreases with casting speed and the same trend is observed between film thickness and slope in the results for Case A1–A10. The calculated slag consumptions show that they match better in cases with higher casting speed because the estimation of slag film thickness becomes more accurate with thinner slag film (
when
h << 1), which can be obtained with higher casting speeds as shown in
Figure 12a. Namely, other factors such as gravity, pressure gradient, and transient effect start contributing to the slag consumption as the slag film gets thicker. Thus, an advanced model to estimate slag film thickness will be required to improve the accuracy of cases with thick liquid films (
).
6. Timing of Slag Consumption
The validated analytical model is now applied to the meniscus region to investigate the timing of slag consumption. Instead of using Equation (21) and iterations to estimate the film thickness and slope, they are assumed to
= 2 mm with slope
s = −0.0033, which are externally given based on previous 2D model calculations [
29]. The Fourier number is used as an indicator to quantify the relationship between momentum transfer, slag film thickness, and mold oscillation frequency discussed in
Figure 4. The Fourier number is defined as follows:
where,
is the kinematic viscosity of slag (
),
is the characteristic time,
is the characteristic length. Here, the period of mold oscillation (
) and slag film thickness (
) are used for the
and
respectively. Note that each Fourier number represents numerous cases that have different combinations of slag viscosity, mold oscillation setting, and slag film thickness. In this parametric study, the kinematic viscosity is used to control the Fourier number. For the mold oscillation setting, conditions from Case A4 (
Table 1) are applied.
Figure 14 shows the normalized slag consumption with different Fourier numbers and the transient velocity profiles for
Fo = 0.027. The slag consumption is normalized to understand the general trend qualitatively.
From the consumption profile in
Figure 14a, a delay of peak is observed as the Fourier number decreases: The peak of slag consumption comes after the peak of descending mold velocity (at period = 0.5). A transition of velocity profile due to the slow response of slag flow to mold oscillation generates this kind of delay as shown in
Figure 14b near the mold wall (
y < 0.3). On the other hand, the force from slag pressure does not show any delay regardless of Fourier number and it responds precisely to the mold oscillation as shown in
Figure 10 and
Figure 13. Also, the negative consumptions observed during the PST in high Fourier numbers turn to positive and the contribution of PST to the total slag consumption becomes higher as the Fourier number decreases. The transient velocity profiles for
= 0.027 show a separation of regions where mold oscillation dominates (
y < 0.3), and shell motion dominates (
y > 0.3). The mold oscillation dominant region decreases as the Fourier number decreases and the stable positive consumption is obtained from the shell motion dominant region. Consequently, the fluctuation of overall consumption reduces over time, and it flattens the consumption pattern. This phenomenon explains the discrepancy between previous results from numerical models. Some works showed positive consumption during the entire oscillation cycle [
17,
29] while others showed periodic variation between negative and positive consumption [
14,
15,
18]. According to the results from this parametric study, the consumption pattern depends on the conditions which can be represented by the Fourier number. A thicker slag film with higher mold frequency and lower slag viscosity limit the influence of mold oscillation over the entire film thickness, and the region driven by steady shell motion allows stable positive infiltration as shown in
Figure 14b. Accordingly, this results in higher consumption during the PST.
In addition to the transient effect represented by Fourier number, the influence of pressure gradient effect on timing of slag consumption is investigated with different thickness ratios of slag film channel between inlet
and outlet
.
Figure 15 shows the transient consumption during the mold oscillation cycle with different film thickness ratios.
The contribution of pressure-driven flow induced by the continuity effect (discussed in
Section 2) increases with the thickness ratio and it turns out to delay the peak of consumption from the peak of descending mold velocity (at period
= 0.5). The results show that the peak is shifted to outside of NST when the thickness ratio becomes greater than
(in
Figure 15a). The transient velocity profiles for
in
Figure 15b visualize parabolic flow patterns against the mold wall velocity: The adverse pressure gradient developing near the inlet in the converging film (e.g.,
Figure 2a) dominates the flow and generates a backflow (e.g., velocity profile at
in
Figure 15b). This backflow against the mold wall motion deviates the consumption pattern from the mold oscillation (i.e., the delay of the consumption peak). Similar behavior is reported from 2D numerical model results recently [
30]. On the contrary, no delay is observed for the force from slag pressure like the parametric study for transient effect with Fourier number.
Therefore, this implies that there is a time difference between slag consumption and mold–strand gap behavior near the meniscus, which aligns with infiltration mechanisms discussed in previous studies [
8,
10,
12]. The current parametric studies seem to reveal that the time difference (lagging or phase shift) is attributed to the transient effect and pressure gradient effect, according to a fluid mechanics perspective. Compared to the corresponding behavior in thin slag films with sinusoidal oscillations (e.g., Case A2 in
Figure 10) that the force and consumption both reach their peak at the maximum descending velocity of mold wall (period
), the slag flow starts showing some lagging behavior compared to its pressure as the transient effect (represented by
) or pressure gradient effect (represented by
) increases near the meniscus. Since the mold–strand gap opens and closes based on the slag pressure which precisely responds to the mold oscillation, it naturally gives diverse consumption patterns (timing for peak consumption, ratio of consumption during PST and NST, etc.) depending on the extent of mismatch between the slag pressure and slag flow. Therefore, the consideration of time difference between slag and gap behavior will be important for the optimization of mold oscillation. A coupling with a model for transient mold–strand gap will be necessary to investigate this phenomenon further in the future.