Numerical Study of Binary Trickle Flow of Liquid Iron and Molten Slag in Coke Bed by Smoothed Particle Hydrodynamics

In the bottom region of blast furnaces during the ironmaking process, the liquid iron and molten slag drip into the coke bed by the action of gravity. In this study, a practical multi-interfacial smoothed particle hydrodynamics (SPH) simulation is carried out to track the complex liquid transient dripping behavior involving two immiscible phases in the coke bed. Numerical simulations were performed for different conditions corresponding to different values of wettability force between molten slag and cokes. The predicted dripping velocity changes and interfacial shape were investigated. The relaxation of the surface force of liquid iron plays a significant role in the dripping rate; i.e., the molten slag on the cokes acts as a lubricant against liquid iron flow. If the attractive force between the coke and slag is smaller than the gravitational force, the slag then drops together with the liquid iron. When the attractive force between the coke and slag becomes dominant, the iron-slag interface will be preferentially detached. These results indicate that transient interface morphology is formed by the balance between the momentum of the melt and the force acting on each interface.


Introduction
An efficient trickle flow of high-temperature melts in coke beds is necessary in ironmaking blast furnaces, as it ensures a smooth continuous process, which is a prerequisite for high productivity in operations that use lower amounts of reducing agents. The trickle flow is driven by the gravitational force that acts on the liquid iron and molten slag (two immiscible liquids), and the dripping behavior is balanced by viscous and interfacial forces. When the volume of the high-strength coke, which acts as a spacer in the lower part of the blast furnace, is reduced, the resistance to the melt passage increases, which can cause instability. It is, therefore, essential to maintain a desired trickle flow through the coke bed to prevent clogging problems that may occur if the pressure drop in the bed becomes excessive or

Basic Formulation of SPH
From the basic idea of SPH, it is evident that fluid motion is represented using a set of particle motion equations. A kernel function is introduced as an integral interpolation to solve a differential equation. In other words, the formulation is based on an interpolation scheme that allows the estimation of a vector or scalar function f at position r in terms of the values of the function at the discretization points.
In Equation (1), r denotes arbitrary coordinates, r denotes particle positions, V denotes the volume, W is the smoothing kernel function, and h is the radius of influence. The summation of function f can be replaced with a summation over particles only within the distance h from r i owing to the compactness; thus, W r ij , h = 0 when r ij > h. The kernel must possess a form symmetrical to r ij = 0. Here, i is the particle index, and j is the index of the neighboring particle around i. The kernel has at least a continuous first derivative and must satisfy the normalization condition as W r ij , h dr = 1.
Within the h → 0 limit, the kernel is required to reduce to a Dirac delta function δ r ij . Wendland's kernel can be applied to prevent having various kernel artifacts in a multiphase system [32]: where r ij 0 is the interparticle distance corresponding to the initial conditions, q = r ij /h, and it is assumed that h = 1.05 r ij 0 [27]. The gradient form of Equation (1) can be represented by using the divergence theorem as follows: When applying this approximation to dispersed fluids, the discontinuities in the density distribution of the fluids become significant, which increases numerical errors. Nonuniform distribution or insufficiency of particles in the regions near the interface results in a significant fluctuation of pressure. The moving least squares (MLS) method is useful for approximating the function to solve this problem [33], which is related to the numerical fluctuations in the pressure at the nearby interface. The pressure is a function of the local density, and thus, the smooth density field of a bulk phase leads to a continuous pressure distribution. The MLS method improves the mass-area-density consistency and filters small-scale pressure oscillations, as described briefly in the following section.

Density Approximation
The SPH formulation can be transformed into a particle-based format to express the mass-area-density consistency process. The density of the particles was expressed in terms of the sum of the kernel functions of the N particles present within the radius of influence as follows: where the subscripts i and j denote the particle indices, m is the mass, and ρ is the local density around the particle. Further, m j is the mass of particle j. Therefore, the kernel function around particle i can be discretized using the following equation, which is derived from Equations (1)- (3): The gradient of f i can be represented as expressed in Equation (6): MLS involves a first-order consistent gradient approximation, which allows pressure smoothing, and its first derivative values are obtained using the method for the homogenous bulk phase mentioned above. The method of the least square interplant with constraint condition (CLS) represents an improved scheme, leading to a more accurate approximation than the MLS method around the sampling points [33,34]. In the CLS method, the particles can be made to directly represent a physical quantity by extending the MLS method in the one-dimensional error space for multiple dimensions. In the three-dimensional space, the CLS method approximates the values of various physical parameters around the particles based on the following equation: where, x, y, and z are the coordinates of the sampling points and a 0 , a 1 , a 2 , and a 3 are the undetermined coefficients. If x = x i , y = y i , and z = z i , then ρ i = a 0 at particle i. Refer to a previous report for parameters determining the procedure [26].

Fluid Motion Equation and Discretization
The governing equations for a weakly compressible viscous flow are based on the relationship between the velocity of sound and the flow density under adiabatic conditions, as well as the Navier-Stokes Equations: where v is the fluid velocity, p is the pressure, c is the velocity of sound, µ is the viscosity, and F s is the interfacial force. Subsequently, Equation (8) can be formulated for each particle as follows: where Π is the artificial viscosity term, which is usually added to the pressure gradient term to help in diffusing sharp variations in the flow and dissipate the energy of the high-frequency term [35].
To determine the time derivative of pressure from Equation (7), Tait's equation of state can generally be used [28]: where γ (= 7.0) is the adiabatic exponent and ρ 0 is the true density value of the material. Considering the balance of the time step and the incompressible behavior of the artificial compressible fluid, an optimal value for c must exist.

Interfacial Force Model
Considering the interfacial force F s , the interparticle potential force is defined using the space derivative of potential E r ij . F s is localized at the liquid interface by applying it to the liquid elements in the transition region of the interface. The force per unit area F s i is then converted into force per unit volume using the expression [36]: where σ i is the surface tension or interfacial tension of particle i. The Fowkes hypothesis is considered in calculating the interfacial force on the multiphase boundary [37]. The Fowkes hypothesis explains that, in a system in which two immiscible liquid phases (liquid iron and molten slag) are in contact, the elements present at the two-phase interface are subject to forces. At the interface between liquid iron and molten slag, liquid iron interface elements receive the attractive force σ m equivalent to the "surface tension" of liquid iron and the dispersion force σ D from molten slag. The force acting on the interface elements of the molten slag can be described similarly. Hence, the interfacial tension σ ms is expressed as follows: This simple hypothesis indicates that the unknown dispersion force and interfacial tension can be calculated explicitly by applying the surface tension as the input and the interfacial tension of the two liquid phases in contact as the conditions. An immiscible blend of liquid iron and molten slag contacting the coke plate is considered, as illustrated in Figure 1. In terms of the tension balance on the solid-gas-liquid triple line in which liquid iron, coke, and gas are in contact with one another, the surface tension σ m of the liquid iron, surface tension σ c of the solid phase, and the solid-liquid interfacial tension σ mc are assumed to be balanced by the contact angle θ m . In other words, Young's equation reflects a horizontal balance among the interfacial tensions: σ m cos θ m + σ mc = σ c . Here, the unknown solid surface tension and the solid-liquid interfacial tension are eliminated from Young's and Fowkes' equations to obtain the following equation: Processes 2020, 7, x FOR PEER REVIEW 5 of 14 elements in the transition region of the interface. The force per unit area 〈 〉 is then converted into force per unit volume using the expression [36]: where is the surface tension or interfacial tension of particle i. The Fowkes hypothesis is considered in calculating the interfacial force on the multiphase boundary [37]. The Fowkes hypothesis explains that, in a system in which two immiscible liquid phases (liquid iron and molten slag) are in contact, the elements present at the two-phase interface are subject to forces. At the interface between liquid iron and molten slag, liquid iron interface elements receive the attractive force σm equivalent to the "surface tension" of liquid iron and the dispersion force from molten slag. The force acting on the interface elements of the molten slag can be described similarly. Hence, the interfacial tension is expressed as follows: This simple hypothesis indicates that the unknown dispersion force and interfacial tension can be calculated explicitly by applying the surface tension as the input and the interfacial tension of the two liquid phases in contact as the conditions. An immiscible blend of liquid iron and molten slag contacting the coke plate is considered, as illustrated in Figure 1. In terms of the tension balance on the solid-gas-liquid triple line in which liquid iron, coke, and gas are in contact with one another, the surface tension of the liquid iron, surface tension of the solid phase, and the solid-liquid interfacial tension are assumed to be balanced by the contact angle . In other words, Young's equation reflects a horizontal balance among the interfacial tensions: Here, the unknown solid surface tension and the solid-liquid interfacial tension are eliminated from Young's and Fowkes' equations to obtain the following equation: Equation (13) indicates that is determined by the surface tension of the liquid phase and the dispersion force acting between the different phases. The dispersion force is explicitly defined by this equation, and the static contact angle can be calculated using the potential interparticle model. Furthermore, considering other triple lines, such as that existing between the molten slag, the solid, and the gas, and that between the two liquid phases and the solid, the following equation is obtained: Equation (13) indicates that θ m is determined by the surface tension of the liquid phase and the dispersion force acting between the different phases. The dispersion force is explicitly defined by this equation, and the static contact angle can be calculated using the potential interparticle model. Furthermore, considering other triple lines, such as that existing between the molten slag, the solid, and the gas, and that between the two liquid phases and the solid, the following equation is obtained: In Equation (14), θ s is the contact angle between the molten slag and solid plate, and θ ms is the contact angle between the two liquid phases and the solid plate. This equation indicates that the liquid iron-molten slag-solid contact angle θ ms is represented by the contact angles θ m and θ s .

Physical Properties of Liquid Iron and Molten Slag
The condition of the lower part of the blast furnace was examined to determine the physical properties of liquid iron (ρ m , µ m , and σ m ) and molten slag (ρ s , µ s , and σ s ). Assuming that the molten iron at the lower part of the blast furnace has a sufficiently low oxygen concentration and is in a carbon-saturated state, the molten iron can be regarded as chemically stable at 1773 K. Hence, it can be assumed as follows: ρ m = 6800 kg/m 3 , µ m = 0.01 Pa s, and σ m = 1.25 N/m [38]. The contact angle between carbon-saturated iron and coke is reported to exceed 120 • . However, the physical properties of molten slag vary over a wide range. The typical slag composition in the dripping zone depends on the CaO-SiO 2 -Al 2 O 3 -MgO system at 1773 K; the basicity (CaO/SiO 2 mass %) ranges from 0.7-2.0 and decreases as the slag descends in the furnace. When the slag contains MgO, the reduction of MgO may proceed preferentially over that of SiO 2 , thereby influencing the wetting behavior of the molten slag on the coke substrate [8]. In this study, the liquid phase composition was considered as the simplest 40 mol% CaO-40 mol% SiO 2 -20 mol% Al 2 O 3 -0 mol% MgO slag, in which the physical property data were assumed in a single phase at a constant temperature of 1773 K. Since the physical properties of molten slag exhibit large deviations among the reported values, typical datasets from the same research group are adopted [39,40]: ρ s = 2574 kg/m 3 , µ s = 0.53 Pa s, and σ s = 0.49 N/m. According to a previous study, there is no agreement on the contact angle between the CaO-SiO 2 -Al 2 O 3 -based molten slag and the carbonaceous material. On one hand, the contact angle decreases when CaO/SiO 2 1 [11,41]; on the other hand, the contact angle remains constant at approximately 160 • in the case of graphite, regardless of the basicity [10]. These differences are thought to be due to the different compositions of the carbonaceous material. The change in these contact angles might be due to the reduction reaction of SiO 2 between the carbonaceous material and molten slag, as well as the formation of SiC as the interfacial product because an initial contact angle of 45 • between SiC and the slag was reported [11]. Therefore, in this study, the contact angle between the molten slag and coke was defined by taking θ m = 160 • for low-reactivity interface and θ m = 30 • for high-reactivity interface [8,10].

Calculation Condition
Detailed packed-bed-structure digital data consisting of multiple coke samples was constructed to simulate the lower part of the blast furnace. Assuming that the coke distribution was just above the raceway, representative coke samples with an average equivalent spherical volume diameter D = 0.0247 m was selected [42]. By using a three-dimensional scanning technique [26], 100 pieces of representative coke three-dimensional surface shape dataset were numerically obtained. About 300,000 surface points on each coke sample were obtained with a minimum resolution of 0.43 mm. The obtained coordinates were converted to standard triangulated language, and the surface shape was polygonal with a triangular mesh. The natural packed structure comprising these particles can be obtained by a DEM-based scheme. The basic format of the DEM is to track spherical particles. Multi-sphere (MS) DEM is a method utilizing a DEM contact force model that is expanded to handle the motion of freely shaped solids. It arrays spherical particles and expresses complex shapes to enable intuitive mounting. The position and rotation angles of each coke sample were determined by using a pseudorandom number, and the packed bed structure was determined by the MS-DEM simulation of a box-type container with 0.12-m sides, similar as a previous report [43]. Next, two immiscible liquids with a 0.12-m width and 0.02-m height were placed immediately above the packed bed to simulate the liquid iron-molten-slag trickle flow in the coke-packed bed. Position of liquid iron and molten slag was determined by using a pseudorandom number. The volume ratio was set as 8:2 to achieve the conditions similar to that of the actual operation. In SPH simulations for multi-phase flow that include a gas, the pressure differential becomes large between liquid and gas, and, consequently, the pressure gradient becomes excessive, thus making convergent calculation difficult. Therefore, in this research, a gas phase is assumed to exist in spaces where particles do not exist. The calculation domain is depicted in Figure 2. The well-known Courant-Friedrichs-Lewy condition was applied to determine dt. In this study, a calculation particle diameter d p of 1.00 mm was adopted as a constant value in all calculation processes. Thus, the following values were determined: dt = 1.0 × 10 −5 s, the analysis time is 10 s, and the number of total particles is 1,402,280. All the programs were coded by the author. Each computer code was written in Fortran 90/95 and compiled and executed by an Intel Fortran compiler on a Windows system. The CPU used in this work was Intel Core i7-7820 × (3.6 GHz, 8 cores). A basic parallelized algorithm was applied using a multi-core processor and OpenMP. Simulation time was 1 week or more.
Processes 2020, 7, x FOR PEER REVIEW 7 of 14 Therefore, in this research, a gas phase is assumed to exist in spaces where particles do not exist. The calculation domain is depicted in Figure 2. The well-known Courant-Friedrichs-Lewy condition was applied to determine dt. In this study, a calculation particle diameter dp of 1.00 mm was adopted as a constant value in all calculation processes. Thus, the following values were determined: = 1.0 × 10 s, the analysis time is 10 s, and the number of total particles is 1,402,280. All the programs were coded by the author. Each computer code was written in Fortran 90/95 and compiled and executed by an Intel Fortran compiler on a Windows system. The CPU used in this work was Intel Core i7-7820 × (3.6 GHz, 8 cores). A basic parallelized algorithm was applied using a multi-core processor and OpenMP. Simulation time was 1 week or more.  Figure 3 shows the temporal change in the three-dimensional distributions of the liquid iron (blue color) and molten slag (green color) as a result of the model calculation. The packed bed structure formed by irregularly shaped cokes is geometrically complex, and the liquid iron flow through the continuous void drops in the form of strings or droplets. From the viewpoint of localized flow behavior, several regions do not receive any flow, and it is assumed that liquid iron flows only through the limited region. However, the flow behavior of liquid iron is influenced by the presence of slag with different coke wettability, as shown in Figures 3b,c. As shown in Figure 3b, for a coke surface with poor wettability containing molten slag, the molten slag drops along with the liquid iron. However, the coke surface with good wettability shown in Figure 3c retains molten slag at the upper part of the packed bed owing to a high attractive force acting against the molten slag.  Figure 3 shows the temporal change in the three-dimensional distributions of the liquid iron (blue color) and molten slag (green color) as a result of the model calculation. The packed bed structure formed by irregularly shaped cokes is geometrically complex, and the liquid iron flow through the continuous void drops in the form of strings or droplets. From the viewpoint of localized flow behavior, several regions do not receive any flow, and it is assumed that liquid iron flows only through the limited region. However, the flow behavior of liquid iron is influenced by the presence of slag with different coke wettability, as shown in Figure 3b,c. As shown in Figure 3b, for a coke surface with poor wettability containing molten slag, the molten slag drops along with the liquid iron. However, the coke surface with good wettability shown in Figure 3c retains molten slag at the upper part of the packed bed owing to a high attractive force acting against the molten slag.   Figure 4a, liquid iron almost has the same distribution at t = 0.20 s, irrespective of the presence of slag and the change in θs. Conversely, at t = 1.20 s, the liquid iron with slag exists at a lower height compared to "without slag" iron. In the region shown in (a-1), liquid iron is concentrated regardless of θs, and in the region shown in (a-2), liquid iron is concentrated when θs = 30°. In Figure 4b, molten slag exists at a lower height when θs = 160° than that at 30°. The attractive force due to the wettability of the coke surface prevents slag from dripping in the direction of gravity. This change in wettability affects the flow form of the liquid iron shown in (a-1) and (a-2). It is necessary to determine the time change of each melt velocity to consider the dripping mechanism leading to this result. The average velocity in the gravitational direction per unit time is given, considering the temporal change in the center of gravity of each phase, as follows [44]:  Figure 4 shows the volume distribution of liquid iron and molten slag in the height direction at t = 0.20 and 1.20 s. In Figure 4a, liquid iron almost has the same distribution at t = 0.20 s, irrespective of the presence of slag and the change in θs. Conversely, at t = 1.20 s, the liquid iron with slag exists at a lower height compared to "without slag" iron. In the region shown in (a-1), liquid iron is concentrated regardless of θs, and in the region shown in (a-2), liquid iron is concentrated when θs = 30 • . In Figure 4b, molten slag exists at a lower height when θs = 160 • than that at 30 • . The attractive force due to the wettability of the coke surface prevents slag from dripping in the direction of gravity. This change in wettability affects the flow form of the liquid iron shown in (a-1) and (a-2). It is necessary to determine the time change of each melt velocity to consider the dripping mechanism leading to this result. The average velocity in the gravitational direction per unit time is given, considering the temporal change in the center of gravity of each phase, as follows [44]:   [45] and also appears to follow the Darcy-type equation. From a detailed perspective, however, the liquid iron tends to increase the dripping rate because of the presence of slag. Moreover, the velocity of liquid iron is affected by the wettability between the slag and coke. The angle θs = 30° yields the maximum liquid-iron velocity within range (A); in contrast, θs = 160° represents the maximum velocity in region (B). This reversal occurs at t = 0.31 s. As observed from Figure 4, the liquid iron dripping rate is higher than that of molten slag. Due to the difference in density between the two melts, the gravitational force of liquid iron was 2.7 times higher than that of molten slag. As molten slag has a high viscosity and low density, it appears that it prevents liquid iron from descending. However, when θs = 30°, the coke bed surface is wet by slag immediately and facilitates the sliding of the liquid iron. With time, this "lubrication" by molten slag is lost, because the liquid iron moves below the molten slag. When θs = 160°, the effect of wetting becomes weak, and the molten slag drips along with the liquid iron. Thus, the "lubrication" between the molten slag and liquid iron remains on the coke surface. Since the molten slag drops slowly, owing to the 0.37 times density and 53 times viscosity coefficient of the liquid iron, these observations appear as slight differences.  Figure 5 displays the average velocity of each melt. At t > 0.5 s, the liquid iron almost reaches a steady state and almost corresponds to the flow at the lower part of the blast furnace presented by Sugiyama et al. [45] and also appears to follow the Darcy-type equation. From a detailed perspective, however, the liquid iron tends to increase the dripping rate because of the presence of slag. Moreover, the velocity of liquid iron is affected by the wettability between the slag and coke. The angle θs = 30 • yields the maximum liquid-iron velocity within range (A); in contrast, θs = 160 • represents the maximum velocity in region (B). This reversal occurs at t = 0.31 s. As observed from Figure 4, the liquid iron dripping rate is higher than that of molten slag. Due to the difference in density between the two melts, the gravitational force of liquid iron was 2.7 times higher than that of molten slag. As molten slag has a high viscosity and low density, it appears that it prevents liquid iron from descending. However, when θs = 30 • , the coke bed surface is wet by slag immediately and facilitates the sliding of the liquid iron. With time, this "lubrication" by molten slag is lost, because the liquid iron moves below the molten slag. When θs = 160 • , the effect of wetting becomes weak, and the molten slag drips along with the liquid iron. Thus, the "lubrication" between the molten slag and liquid iron remains on the coke surface. Since the molten slag drops slowly, owing to the 0.37 times density and 53 times viscosity coefficient of the liquid iron, these observations appear as slight differences. The variation in the wettability between the coke and molten slag modifies the interfacial area between the four-phase iron-slag-coke-gas by the following mechanism. When some droplets are dispersed in the system under isothermal conditions, the Helmholtz free energy, = ∑ , of the system increases. Since the potential energy corresponding to the flow of the liquid iron droplets through the packed bed is equal for each case, the kinetic energy corresponding to the secondary droplet formation may be the same for all the conditions. However, the total energy in the system is not conserved due to the viscous damping from this calculation, but it is useful to clarify the effect of wettability between the coke and slag on the interfacial area of each phase. As our interest is focused on the effect of wettability between the slag and coke on the interfacial energy of each phase, the interfacial area for each phase should be estimated. The main advantage of the SPH method is its rapid prediction of the interface area from its initial condition and the interface-judged particle number .

Results
is geometrically determined from the initial conditions. In this study, the free surface of the liquid iron and slag are 0.04992 and 0.01248 m 2 , respectively, and the initial iron-slag interface area is 0.02110 m 2 . See Section 2 and a previous report [46] for the counting procedure of n. Figure 6 depicts the time change of the estimated interfacial area of each phase. At t = 0.1 s, the melt penetrates the coke bed, and both interfacial areas significantly vary simultaneously. In Figure 6a, for "without slag" liquid iron, the initial iron-gas surface area is significantly different from that of "with slag", but the area becomes almost equal to other conditions at t = 0.1 s. The case of θs = 160° sometimes indicates a larger iron-gas surface area than that of "without slag", because the iron slides down first and forms a contact interface at a lower height. Beyond t = 0.1 s, "without slag" iron has the lowest iron-gas surface area. This behavior can be explained as follows. As shown in Figure 6b, since the liquid is not affected by molten slag in "without slag" liquid iron, the iron-coke interfacial area significantly increases compared to that of "with slag" at t > 0.1 s; in other words, the iron-coke interface decreases with the presence of the molten slag. From Figure 6c, although the iron-slag interface area has a larger The variation in the wettability between the coke and molten slag modifies the interfacial area between the four-phase iron-slag-coke-gas by the following mechanism. When some droplets are dispersed in the system under isothermal conditions, the Helmholtz free energy, F = σ ij A i , of the system increases. Since the potential energy corresponding to the flow of the liquid iron droplets through the packed bed is equal for each case, the kinetic energy corresponding to the secondary droplet formation may be the same for all the conditions. However, the total energy in the system is not conserved due to the viscous damping from this calculation, but it is useful to clarify the effect of wettability between the coke and slag on the interfacial area of each phase. As our interest is focused on the effect of wettability between the slag and coke on the interfacial energy of each phase, the interfacial area for each phase should be estimated. The main advantage of the SPH method is its rapid prediction of the interface area A from its initial condition A 0 and the interface-judged particle number n.
A 0 is geometrically determined from the initial conditions. In this study, the free surface of the liquid iron and slag are 0.04992 and 0.01248 m 2 , respectively, and the initial iron-slag interface area is 0.02110 m 2 . See Section 2 and a previous report [46] for the counting procedure of n. Figure 6 depicts the time change of the estimated interfacial area of each phase. At t = 0.1 s, the melt penetrates the coke bed, and both interfacial areas significantly vary simultaneously. In Figure 6a, for "without slag" liquid iron, the initial iron-gas surface area is significantly different from that of "with slag", but the area becomes almost equal to other conditions at t = 0.1 s. The case of θs = 160 • sometimes indicates a larger iron-gas surface area than that of "without slag", because the iron slides down first and forms a contact interface at a lower height. Beyond t = 0.1 s, "without slag" iron has the lowest iron-gas surface area. This behavior can be explained as follows. As shown in Figure 6b, since the liquid is not affected by molten slag in "without slag" liquid iron, the iron-coke interfacial area significantly increases compared to that of "with slag" at t > 0.1 s; in other words, the iron-coke interface decreases with the presence of the molten slag. From Figure 6c, although the iron-slag interface area has a larger value at θs = 160 • than at θs = 30 • , the slag-coke interface area exhibits the opposite trend, as observed in Figure 6d. Since the coke surface with good wettability is covered with slag, it will likely increase the contact area between the free surface of the slag and iron-slag in a static state. In the case of θs = 160 • , the slag drops with the liquid iron due to gravity. This dripping occurs because the attractive force acting between the coke and slag is less than the gravitational force. However, for θs = 30 • , the attractive force acting between the coke-slag interface dominates, and the iron-slag interface is preferentially detached. Here, the effect of momentum in the direction of gravity appears. Since a "quasi-stable" interface is formed by the balance between the melt momentum and the force acting on each interface, the state predicted by the equilibrium theory is not always achieved during dripping [47]. This result indicates that static hold-up of the molten slag may promote the smooth dripping of liquid iron, as variations in the wettability between the coke and slag transiently change the flow field of the liquid iron. It emphasizes the distinctive interfacial features that can depend on interactions with the wettability of nonuniform packed cokes and the transient behavior of two melts. In the future, the calculation region of this model can be expanded and the model may be applied to continuous dripping, and the application of the model can also be extended to the entire bottom of the blast furnace.
value at θs = 160° than at θs = 30°, the slag-coke interface area exhibits the opposite trend, as observed in Figure 6d. Since the coke surface with good wettability is covered with slag, it will likely increase the contact area between the free surface of the slag and iron-slag in a static state. In the case of θs = 160°, the slag drops with the liquid iron due to gravity. This dripping occurs because the attractive force acting between the coke and slag is less than the gravitational force. However, for θs = 30°, the attractive force acting between the coke-slag interface dominates, and the iron-slag interface is preferentially detached. Here, the effect of momentum in the direction of gravity appears. Since a "quasi-stable" interface is formed by the balance between the melt momentum and the force acting on each interface, the state predicted by the equilibrium theory is not always achieved during dripping [47]. This result indicates that static hold-up of the molten slag may promote the smooth dripping of liquid iron, as variations in the wettability between the coke and slag transiently change the flow field of the liquid iron. It emphasizes the distinctive interfacial features that can depend on interactions with the wettability of nonuniform packed cokes and the transient behavior of two melts. In the future, the calculation region of this model can be expanded and the model may be applied to continuous dripping, and the application of the model can also be extended to the entire bottom of the blast furnace.

Conclusions
Numerical simulations using a multi-phase SPH method were performed on two immiscible high-temperature melts dripping through a coke-packed bed during an ironmaking blast furnace process, namely, a binary trickle flow in which liquid iron and molten slag dripped simultaneously. The advantage of this method is the direct estimation of the transient three-dimensional behavior of multi-phase flow based on body forces considering the force acting between different phases, including complex dispersed phases. The present detailed analysis helps in the rationalization of the unexplained transient behavior of liquid iron and molten slag presented by previous experiments. In particular, the effect of coke surface wettability for molten slag on binary trickle flow was investigated.
The flow behavior of liquid iron is influenced by the presence of slag in coke surfaces with different coke wettability. If the coke surface has poor wettability and contains molten slag, the molten slag drops along the molten iron. However, a coke surface with good wettability retains molten slag owing to a higher attractive force acting between the coke and slag rather than an interfacial force of attraction acting between the iron and slag. Although the molten slag dripping rate is lower than that of the liquid iron because of high viscosity and low density, the coke bed surface is wet by slag immediately and facilitates the sliding of the liquid iron. These results demonstrate that the static hold-up of molten slag may promote the smooth dripping of liquid iron. As predicted by novel findings presented in this study, when the low reducing agent ratio operation is realized, increases in the amount of liquid iron and molten slag passing through a unit volume of coke bed finally influence the amount of stagnant melt behavior in the coke bed through a mechanism in which the descent velocity of each melt depends on coke-slag wettability. Through this analysis, the mechanism underlying the simultaneous trickle flow of liquid iron-molten slag behaviors, as described in conventional research, can be explained from a unified point of view.