Optimization Design of Submerged-Entry-Nozzle Structure Using NSGA-II Genetic Algorithm in Ultra-Large Beam-Blank Continuous-Casting Molds

To achieve uniform cooling and effective homogenization control in ultra-large beam-blank molds necessitates the optimization of submerged-entry-nozzle (SEN) structures. This study employed computational fluid dynamic (CFD) modeling to investigate the impact of two-port and three-port SEN configurations on fluid flow characteristics, free-surface velocities, temperature fields, and solidification behaviors. Subsequently, integrating numerical simulations with the non-dominated sorting genetic algorithm II (NSGA-II) and metallurgical quality-control expertise facilitated the multi-objective optimization of a three-port SEN structure suitable for beam-blank molds. The optimized parameters enabled the three-port SEN to deliver molten steel to the meniscus at an appropriate velocity while mitigating harmful effects such as SEN port backflow, excessive surface temperature differences, and shell thickness reduction due to fluid flow. The results indicated that the three-port SEN enhanced the molten-steel flow pattern and mitigated localized shell thinning compared with the two-port SEN. Additionally, the optimized design (op2) of the three-port SEN exhibited reduced boundary layer separation and superior fluid dynamics performance over the initial three-port SEN configuration.


Introduction
To meet the growing demands for high strength, toughness, and high reliability, beam blanks are evolving towards larger integral types.However, the complex cross-section of ultra-large beam blanks and the size effects of large billets lead to uneven moltensteel flow in the transition zone between the web and the flange, resulting in uneven solidification.This unevenness can result in steel breakout accidents and longitudinal cracks [1].Additionally, the single-port straight submerged entry nozzle (SEN) commonly employed in current beam-blank continuous-casting molds often allows molten steel to penetrate deeply, thereby hindering the removal of non-metallic inclusions through flotation and preventing efficient slag-powder melting at the meniscus [2][3][4][5].Many product defects in continuous casting correlate with molten-steel flow within the mold [6], which are significantly influenced by SEN structure, thereby affecting inclusion transport, meniscus status, solute transport, and solidification [7,8].To ensure exceptional metallurgical quality during the casting of ultra-large beam-blank molds, a suitable new SEN type was designed.
Previous designs of beam-blank-mold SENs have primarily focused on port number, bore diameter, angle, and immersion depth.Xu et al. [9] demonstrated the impacts of different port numbers on fluid flow, heat transfer, and solidification through water modeling experiments and mathematical models.Additionally, they recommended a three-port SEN to maintain optimal meniscus conditions and consistent shell thickness at the mold exit, which reduces breakout risks from flange depressions and fillet ruptures.Furthermore, during the continuous casting of the beam-blank mold, Chen and Yang et al. [4,10,11] examined the effects of various SEN structural parameters on vortex formation, impact depth, and free-surface velocity of molten steel during continuous casting.Simulation results highlighted the bore diameter, outlet angle, immersion depth, and horizontal-port angles as crucial factors influencing molten-steel flow.Also, optimized SEN parameters according to mass calculations were recommended to enhance beam-blank cleanliness.
Several mold studies have been conducted through mathematical simulation.To elucidate fluid flow patterns inside the mold and the backflow phenomenon, as well as its controlling parameters, Garcia-Hernandez et al. [12] conducted a mathematical simulation of the fluid dynamics in the SEN and mold.Hibbeler et al. [13] developed a mathematical model for simulating the temperature profile and stress/deformation field inside the solidifying shell of a beam-blank continuous-casting mold.The model sheds light on the mechanism of crack formation in the filet region, owing to a combination of a thinner shell, excessive surface temperature differences, and mechanical stress.
Recent advancements have utilized mathematical models and optimization strategies, such as heuristic search, neural networks, genetic algorithms, and knowledge bases, which have been shown to effectively optimize product quality.Chakraborti [14][15][16] utilized genetic algorithms for optimizing continuous-casting speed through process variable adjustments.At a steel factory in Slovenia, Santos [17], developed a practical parameteroptimization approach employing a genetic algorithm and numerical simulation for the spray cooling zone during the continuous-casting process.
This study established a computational fluid dynamic (CFD) model for the solidification process of molten steel in an ultra-large beam-blank mold (900 × 510).The study investigated how two-port and three-port SEN configurations affect transient fluid flows, heat transfer, and solidification.Upon building on the research results from two SEN simulations, an interactive model integrating numerical simulations of three-port SENs with the non-dominated sorting genetic algorithm II (NSGA-II) multi-objective optimization algorithm was proposed.A knowledge base for metallurgical quality control, incorporating constraints on the process and casting quality in beam-blank molds, was connected to the algorithm.Structural constraints were also incorporated into the knowledge base to integrate design limits for the SEN.The genetic algorithm utilized the knowledge to explore spatial parameter settings and determine optimal three-port SEN structure parameters.This approach ensures compliance with SEN requirements and guarantees high productivity as well as quality in beam-blank continuous-casting molds.

Design of the SENs and CFD Models 2.1. Design of the SENs
This study designed two-port and three-port SENs and conducted numerical simulations to investigate their influence on metallurgical behavior in beam-blank continuouscasting molds.The outlet angle and immersion depth of these SENs were fixed at 4 • and 110 mm, respectively.The geometric schematic and installation methods of SENs, as well as the schematic of the beam-blank transverse section, are shown in Figure 1.

Basic Assumptions of the CFD Models
Given the complexity of the fluid flow, heat transfer, and solidification of molten steel in the mold, the following presumptions were made to develop the mathematical model: (1) The molten steel was treated as an incompressible Newton fluid.
(2) The mold's taper and oscillation do not significantly affect fluid flow patterns.
(3) The latent heat associated with solid phase transition was negligible compared with the latent heat of solidification.

Basic Assumptions of the CFD Models
Given the complexity of the fluid flow, heat transfer, and solidification of molten steel in the mold, the following presumptions were made to develop the mathematical model: (1) The molten steel was treated as an incompressible Newton fluid.
(2) The mold's taper and oscillation do not significantly affect fluid flow patterns.
(3) The latent heat associated with solid phase transition was negligible compared with the latent heat of solidification.

Fluid Flow Model
Two types of fluid flow were considered within the cast strand: interdendritic flow in the mushy zone and turbulent flow in the liquid core [18].The standard k ε − model was employed to simulate the turbulence effects of fluid flow inside the cast strand [19].The following are expressions of the governing equations for fluid flow.
Continuity equation: where ρ is the fluid density ( Momentum equation: Turbulent kinetic energy k and turbulent dissipation rate ε equation:  Continuity equation: where ρ is the fluid density (kg m −3 ) and u is the velocity (m s −1 ).Momentum equation: Turbulent kinetic energy k and turbulent dissipation rate ε equation: where P represents the static pressure (Pa); µ e f f is the effective viscosity, determined by summing the laminar viscosity µ and turbulent viscosity µ t , (kg m −1 s −1 ); g denotes the gravity acceleration (m s −2 ); and G k signifies the turbulence kinetic energy resulting from the mean velocity gradient (m −2 s −2 ).The empirical constants used in the low Reynoldsnumber model were The mushy zone (partially solidified) was treated as a porous medium employing the enthalpy-porosity technique (EPT).The porosity of every cell was determined by the liquid percentage within the cell.Where there is total solidification, there is zero porosity, which results in the disappearance of velocity in those areas [6,20].Additional source terms S p , S k and S ε were added to the momentum conservation equation, turbulent kinetic energy k equation, and turbulent dissipation rate ε equation, respectively, accounting for the mushy zone's decreased porosity: where f l is the liquid fraction; ξ denotes the small constant, set as 0.0001 to prevent floating point division by zero; A represents the mushy zone constant (1 × 10 8 [21]); and u c signifies the casting speed (m s −1 ).

Heat Transfer and Solidification Model
The EPT approach was also applied to simulate the heat transfer and solidification of molten steel [22].
Energy conservation equation: where λ e f f is the effective thermal conductivity, determined by summing the thermal conductivity λ and the turbulent thermal conductivity λ t , (W m −1 K −1 ), and T is the temperature (K).
The following describes the total enthalpy: where h re f is the reference enthalpy (J kg −1 mol −1 ); T re f is the reference temperature (K); C p denotes the specific heat (J kg −1 K −1 ); and L represents the latent heat (J kg −1 ).The energy conservation equation introduced a source term S h to explain the liquidsolid transition and the creation of the mushy zone: Equation ( 14) was used to calculate the liquid fraction: where T l and T s are the liquidus and solidus temperatures of the molten steel (K), respectively.The effective length of the mold was 780 mm.To prevent potential backflow at the mold exit, the computational domain was extended to 1430 mm.A polyhedral mesh structure divided the entire computational area, with local grid refinement applied to the SEN to ensure model accuracy.The domain adjacent to the mold side was enhanced with a mesh boundary layer to more accurately simulate the evolution of the solidified shell.Figure 2a depicts the geometrical model of the two-port SEN computational domain, while Figure 2b shows the mesh partition scheme for the mold.The entire mold contained ~200,000 polyhedral cells.
The physical characteristics of Q235B carbon structural steel, crucial for fluid flow, heat transfer, and solidification in the mushy zone, are presented in Table 1.Temperaturedependent variations in density, specific heat, thermal conductivity, and viscosity were calculated using the JMatpro 14.3 software (Figure 3).The physical characteristics of Q235B carbon structural steel, crucial for fluid flow, heat transfer, and solidification in the mushy zone, are presented in Table 1.Temperaturedependent variations in density, specific heat, thermal conductivity, and viscosity were calculated using the JMatpro 14.3 software (Figure 3).

Boundary and Solution Method
Specific boundary conditions were implemented as follows: (1) At the SEN inlet, a velocity-inlet boundary condition was implemented.The inlet velocity was determined by the billet section and casting speed (Equation ( 15)).
(3) The initial and inlet temperatures of the molten steel were determined by summing the liquidus temperature and the superheat.
(4) In the mold region, the top surface of the mold was adiabatic with zero shear stress.The heat flux across SEN walls was zero, and both mold and SEN walls were considered stationary and nonslip.The average heat flux applied to the mold walls is shown in Equation ( 18).
(5) The outlet assumes fully developed flow conditions, with normal gradients of all variables set to zero.

Boundary and Solution Method
Specific boundary conditions were implemented as follows: (1) At the SEN inlet, a velocity-inlet boundary condition was implemented.The inlet velocity was determined by the billet section and casting speed (Equation ( 15)).
(3) The initial and inlet temperatures of the molten steel were determined by summing the liquidus temperature and the superheat.
(4) In the mold region, the top surface of the mold was adiabatic with zero shear stress.The heat flux across SEN walls was zero, and both mold and SEN walls were considered stationary and nonslip.The average heat flux applied to the mold walls is shown in Equation ( 18).
(5) The outlet assumes fully developed flow conditions, with normal gradients of all variables set to zero.
where S out and S SEN are the areas of the mold outlet and SEN inlet (m 2 ), respectively; D SEN is the hydraulic diameter of the SEN (m); ρ w is the cooling water density (kg m −3 ); C p is the specific heat of the cooling water (J kg −1 K −1 ); Q w is the cooling water flow (m 3 s −1 ); and ∆T water is the temperature difference in the cooling water between the inlet and outlet of the copper tube (K).S e f f represents the effective contact area between the molten steel and mold walls (m 2 ).
The flow and heat transfer calculations incorporate the casting speed, u cast = 0.75 m s −1 , superheat of molten steel, T superheat = 25 K, cooling water flow, Q w = 8800 L min −1 , and temperature difference between the inlet and outlet of the copper tube, ∆T water = 4 K.
Transient calculations were performed using the ANSYS Fluent software, and a pressure-based solver was employed to solve the governing equations.The algorithm for pressure-velocity coupling is known as SIMPLEC.The discretization scheme for the pressure algorithm was PRESTO, while a second-order upwind discretization scheme was applied for momentum.Convergence criteria were met when the energy residual was ≤ 10 −6 , and other residuals were ≤ 10 −4 .

Fluid Flow
To visually illustrate the flow patterns in the mold for two-port and three-port SENs, Figure 4 depicts the velocity vectors on the x-z plane (y = 0 m) and the plane formed by x-z at 70 • .Radial-flow SENs exhibited a "self-braking" effect that enhanced flow stability in the web region [9].Molten steel from two SEN ports was directed toward the web center.Upon collision, a component of the flow increased spirally to the surface while the other half spirally descended, thereby creating the well-known "double roll" pattern [23].The upward stream contributed to an active meniscus with increased temperature.
Figure 4a, depicts the flow pattern of the two-port SEN near the flange.Upon impact with the initially solidified shell of the flange at a high velocity of 1.3 m/s, the SEN jet flow divided into two longitudinal circulations.A strong upper circulation centered at 0.04 m below the meniscus close to the free surface induced significant surface fluctuations.The circular flow below the SEN, centered at 0.3 m below the meniscus and close to the flange, did not facilitate the effective removal of non-metallic inclusions near the mold center [24].Contrastingly, Figure 4b illustrates the injection of molten steel at lower velocities toward the mold corners with the three-port SEN, owing to the longer distance between the nozzle outlets and mold corners.The centers of upper and lower circulations, located at 0.06 m and 0.3 m below the meniscus, respectively, exhibited significantly weaker upper circulations and a shift in the impact zones of circular flow toward the mold center, thereby facilitating meniscus stability and inclusion removal.Transient calculations were performed using the ANSYS Fluent software, and a pressure-based solver was employed to solve the governing equations.The algorithm for pressure-velocity coupling is known as SIMPLEC.The discretization scheme for the pressure algorithm was PRESTO, while a second-order upwind discretization scheme was applied for momentum.Convergence criteria were met when the energy residual was

Fluid Flow
To visually illustrate the flow patterns in the mold for two-port and three-port SENs, Figure 4 depicts the velocity vectors on the x-z plane (y = 0 m ) and the plane formed by x-z at 70  .Radial-flow SENs exhibited a "self-braking" effect that enhanced flow stability in the web region [9].Molten steel from two SEN ports was directed toward the web center.Upon collision, a component of the flow increased spirally to the surface while the other half spirally descended, thereby creating the well-known "double roll" pattern [23].The upward stream contributed to an active meniscus with increased temperature.Figure 5 shows the predicted maximum velocity on the free surface, which was obtained 20 m m below the meniscus, for both SEN types.The molten-steel flow conditions Figure 5 shows the predicted maximum velocity on the free surface, which was obtained 20 mm below the meniscus, for both SEN types.The molten-steel flow conditions below the meniscus significantly influenced free-surface velocity.The two-port SEN exhibited a notably higher rate of change in free-surface velocity near the flange compared with the three-port SEN.A maximum free-surface velocity of 0.2 m/s at the flange center increased the risk of steel/slag interface damage and slag entrainment [25].Additionally, the velocity distribution on the free surface of the three-port SEN was more uniform compared with the two-port SEN, thereby effectively preventing the formation of a dead zone and the freezing of molten steel [26].Thus, comprehensively considering flow characteristics in both models, the three-port SEN exhibited superior metallurgical effects.
below the meniscus significantly influenced free-surface velocity.The two-port SEN exhibited a notably higher rate of change in free-surface velocity near the flange compared with the three-port SEN.A maximum free-surface velocity of 0.2 m/s at the flange center increased the risk of steel/slag interface damage and slag entrainment [25].Additionally, the velocity distribution on the free surface of the three-port SEN was more uniform compared with the two-port SEN, thereby effectively preventing the formation of a dead zone and the freezing of molten steel [26].Thus, comprehensively considering flow characteristics in both models, the three-port SEN exhibited superior metallurgical effects.

Temperature Field and Solidification
Figure 6 illustrates the temperature distribution along the casting direction with two distinct types of SEN.As the molten-steel flow enters the impacting zone, the surface temperature at the center of the two-port SEN flange increased from 1731 K at 45 m m to 1757 K at 192 m m below the meniscus.Upon exiting this zone, the surface tempera- ture decreased.The temperature variation observed with the three-port SEN exhibited similar trends compared with the two-port SEN.Both exhibited localized high-temperature zones.These findings indicated the close relationship between temperature distribution and the flow pattern of molten steel within the mold.Despite variations in different SEN configurations, temperature changes in the impact zone remained consistent.The direct impingement of high-temperature molten steel from SEN ports on the solidification front resulted in localized maximum superheat, with heat dissipating through longitudinal circulation.Notably, the circumferential temperature distribution of the three-port SEN remained more uniform across the web center, flange center, and flange tip center, thereby reducing the possibility of cracks and defects [1].

Temperature Field and Solidification
Figure 6 illustrates the temperature distribution along the casting direction with two distinct types of SEN.As the molten-steel flow enters the impacting zone, the surface temperature at the center of the two-port SEN flange increased from 1731 K at 45 mm to 1757 K at 192 mm below the meniscus.Upon exiting this zone, the surface temperature decreased.The temperature variation observed with the three-port SEN exhibited similar trends compared with the two-port SEN.Both exhibited localized high-temperature zones.These findings indicated the close relationship between temperature distribution and the flow pattern of molten steel within the mold.Despite variations in different SEN configurations, temperature changes in the impact zone remained consistent.The direct impingement of high-temperature molten steel from SEN ports on the solidification front resulted in localized maximum superheat, with heat dissipating through longitudinal circulation.Notably, the circumferential temperature distribution of the three-port SEN remained more uniform across the web center, flange center, and flange tip center, thereby reducing the possibility of cracks and defects [1].
Figure 7 shows the liquid fraction at different cross-sections under both types of SEN.In this study, a liquid fraction of 0.33 defined the solidification front.Under the two-port SEN, the thickness of the solidified shell at the mold-outlet (z = 0.78 m) web center, flange center, and fillet measured 11.75 mm, 12.05 mm, and 8.607 mm, respectively.For the threeport SEN, these values were 11.6 mm, 11.48 mm, and 10.238 mm, indicating a maximum thickness change reduction of 2.081 mm compared with the two-port SEN.This difference was attributable to the fact that ports of the three-port SEN were farther from the mold wall, thereby resulting in lower molten-steel impact velocities compared with the two-port SEN.Consequently, the homogeneity of the solidified shell was less affected by shock pressure from jet flows against the mold wall.Thus, the solidification behavior of the three-port SEN was more acceptable within the mold, potentially reducing the risk of liquid breakout [25].This difference was attributable to the fact that ports of the three-port SEN were farther from the mold wall, thereby resulting in lower molten-steel impact velocities compared with the two-port SEN.Consequently, the homogeneity of the solidified shell was less affected by shock pressure from jet flows against the mold wall.Thus, the solidification behavior of the three-port SEN was more acceptable within the mold, potentially reducing the risk of liquid breakout [25].The multi-objective optimization model integrates the NSGA-II algorithm with CFD models.The NSGA-II algorithm, which encompasses optimization strategies and a set of constraints, was employed to determine the optimization parameters for the SEN struc-

Multi-Objective Optimization Model for Three-Port SEN Structure
The multi-objective optimization model integrates the NSGA-II algorithm with CFD models.The NSGA-II algorithm, which encompasses optimization strategies and a set of constraints, was employed to determine the optimization parameters for the SEN structure and to integrate new input variables into the CFD model.The CFD model simulates the flow, heat transfer, and solidification behavior, with the results being fed back to the NSGA-II algorithm.Given the metallurgical phenomena of these two types of SEN (Section 3), the three-port SEN was selected as the focus of optimization for this study.

Knowledge Base Representation for Metallurgical Quality Control
The criteria below are structured in terms of the flow and thermal behavior of molten steel in the mold, thereby drawing cause-and-effect relationships from the literature.
Hot molten steel with high momentum impinging against the solidified shell resulted in the thinning of the shell, thereby increasing the risk of deformation and breakout under pressure.It was crucial to maintain a shell thickness at the mold exit above a specific minimum value, preventing breakout due to liquid ferrostatic pressure.The minimum acceptable value was ~12 mm.
S m ≥ 12 (2) The velocity of the free surface (V s ) [28,29].An optimal velocity range for molten steel at the free surface was crucial as too low a velocity may result in freezing near the meniscus, impeding proper melting of mold flux.However, an excessively high value can lead to excessive level fluctuation and slag entrainment.An acceptable free-surface velocity range was 0.1-0.4m/s.0.1 ≤ V s ≤ 0.4 (3) The effective area fraction at the port (β) [12,30,31].The effective area fraction represents the flow fraction with a positive directional velocity component at the port area, ensuring it exits the SEN domain effectively.Maximally, the design of the SEN minimized or eliminated backflow at each port area, enhancing the utilization of the effective outlet area, thereby reducing jet velocity and achieving a more stable flow pattern.
(4) The temperature difference between the fillet and inner-flange tip corner at the mold exit (∆T) [32,33].
In continuous-casting mold for beam blanks, the fillet and inner-flange tip corners exhibited the highest and lowest wide surface temperatures, respectively.The temperature difference between them served as an indicator of hot surface temperature uniformity.To prevent, longitudinal cracks at the fillet, the SEN design should minimize the temperature difference at the mold exit.

Definition of the Optimization Issues
Multi-objective optimization design was conducted on the primary parameters of the three-port SEN, guided by relevant performance indicators from the metallurgical quality-control knowledge base.The aim was to obtain a superior SEN structure, enhanced flow, and heat transfer performance.The objective includes maximizing the effective area fraction at the port and minimizing the temperature difference between the fillet and the inner-flange tip corner.Additionally, constraints were placed on the shell thickness at the mold exit and the velocity of the free surface to prevent significant damage.The design variables include the bore diameter (t 1 ), the outlet angle (t 2 ), the immersion depth (t 3 ), and the port's horizontal angles (t 4 ) of the three-port SEN.They ranged from 45 mm to 50 mm for t 1 , from 2 • to 15 • for t 2 , from 110 mm to 126 mm for t 3 , and from 104 • to 124 • for t 4 , respectively.The size of the square port was calculated according to the aspect ratio of 2:1 and equal-area principles relative to the SEN inlet.The optimization issue can be expressed mathematically as follows:

Surrogate Models and Error Metrics
Generally, formulating objective functions for complex nonlinear issues can be challenging.As an effective alternative, the polynomial response surface method (PRSM) was employed to predict response values for unknown points according to the information from known points.This approach facilitated design optimization with relatively low computational cost [34].The PRSM model utilizes a quartic polynomial function in subsequent optimization procedures, defined as Design of experiments (DOE), including full factorial design, Latin hypercube, central composite design, and Taguchi orthogonal array, were employed to generate sampling points for determining PRSM coefficients [35].To construct a high-precision surrogate model, 30 training points were selected in the design space using Latin hypercube sampling.Furthermore, to evaluate the accuracy of the model, five additional validation points distributed over the design domain were randomly generated.The accuracy of the surrogate model was validated using R-square (R 2 ), root-mean-square error (RMSE), maximum absolute percentage error (MAPE), and relative error (RE) (Equations ( 21)-( 24)) [36].
where n denotes the number of sample points; y i and ŷ are the exact function value of the ith sample point and the corresponding value of the surrogate model; and y represents the average value of y i .The values of these performance metrics revealed the predictive accuracy and capability of the surrogate models at new points.The overall accuracy of the approximate model was estimated using R 2 and RSME, while the local accuracy was calculated using MAPE.Generally, the surrogate model is deemed more accurate when R 2 approaches unity, and RMSE and RAPE values are smaller.

Multi-Objective Optimization Algorithm and Procedure
The NSGA-II algorithm, known for its efficiency in multi-objective optimization through non-dominated sorting, was adopted [34,37].In this study, NSGA-II was employed to obtain Pareto fronts of two objectives: ∆T and β.Details of NSGA-II parameter definitions are summarized in Table 2. Furthermore, Figure 8 shows a flowchart illustrating the multi-objective optimization process for the three-port SEN structure.

Error Analysis of Surrogate Models
In this section, the relationship between design variables and performance indicators was established using the DOE and PRSM.Each set of training points was inputted into the Fluent software to calculate their flow and heat transfer processes.Subsequently, multi-objective optimization was conducted (Figure 8) after each simulation.Ultimately, 28 out of the 30 sets of sampling points were reserved as effective data to establish the approximate model for each performance indicator.Equations ( 25)-( 28) derive quartic polynomial functions for m S , s Table 3 presents the modeling accuracies of the PRSM functions established according to the simulation results of the Latin hypercube design, thereby explicitly comparing the performance of these surrogate models.It was evident that the 2 R value was ~1, and the values of RAAE and MAPE were also within 0.2.Additionally, the RE errors of the CFD results and surrogate models at the five verification points were all less than 0.2 (Figure 9).In this study, despite the quartic polynomial not achieving optimal modeling ac-

Error Analysis of Surrogate Models
In this section, the relationship between design variables and performance indicators was established using the DOE and PRSM.Each set of training points was inputted into the Fluent software to calculate their flow and heat transfer processes.Subsequently, multi-objective optimization was conducted (Figure 8) after each simulation.Ultimately, 28 out of the 30 sets of sampling points were reserved as effective data to establish the approximate model for each performance indicator.Equations ( 25)-( 28) derive quartic polynomial functions for S m , V s , β, and ∆T versus the three-port SEN structural parameters x 1 , x 2 , x 3 , and x 4 .Among x 1 , x 2 , x 3 , and x 4 were the normalization processing results of t 1 , t 2 , t 3 , and t 4 (Equation ( 29)).
Table 3 presents the modeling accuracies of the PRSM functions established according to the simulation results of the Latin hypercube design, thereby explicitly comparing the performance of these surrogate models.It was evident that the R 2 value was ~1, and the values of RAAE and MAPE were also within 0.2.Additionally, the RE errors of the CFD results and surrogate models at the five verification points were all less than 0.2 (Figure 9).In this study, despite the quartic polynomial not achieving optimal modeling accuracy, it met the requirements for subsequent analysis.(25) x

Multi-Objective Optimization Results
Figure 10 illustrates the Pareto fronts for optimizing the three-port SEN structure using the NSGA-II algorithm.Notably, two goals-β and ∆T-contradicted one another, as an increase in β invariably increased ∆T.Therefore, the optimal design requires balancing practical engineering needs.Upon focusing solely on β, the optimization of point 1 (op1)'s design features a bore diameter of 45 mm, an outlet angle of 11 • , an immersion depth of 110 mm, and port horizontal angles of 114 • .The actual value β, derived using the CFD model, was 0.9162, and was slightly higher compared with the surrogate model prediction of 0.9106.Conversely, given the single objective ∆T, the op3 design was obtained with a bore diameter of 46 mm, an outlet angle of 9 • , an immersion depth of 111 mm, and port horizontal angles of 112 • .In this design, the values of β and ∆T were 0.8570 and 156.9 K, respectively.
Materials 2024, 17, x FOR PEER REVIEW 16 of 21  To balance β and T Δ , a compromise solution of op2 design was obtained based on the "distance minimization" criterion (Equation (30) [35]).

Multi-Objective Optimization Results
( ) where u is the number of objective functions,  4 presents the surrogate model predictions with the CFD model results to validate their accuracy.It was observed that the maximum relative error of β was <4.26% and the error of T Δ was <12.45% between CFD results and surrogate models.This further indicated that optimization results based on the PRSM model achieved a high level of accuracy.To balance β and ∆T, a compromise solution of op2 design was obtained based on the "distance minimization" criterion (Equation (30) [35]).
where u is the number of objective functions, u = 2, and F k i is the objective function i at the optimization point k.
The optimized SEN featured a bore diameter of 45 mm, an outlet angle of 10 • , an immersion depth of 110 mm, and port horizontal angles of 114 • .In comparison with the initial three-port SEN, β increased by 9.20% and ∆T decreased by 45.12%, and S m and V s , met metallurgical standards.Additionally, Table 4 presents the surrogate model predictions with the CFD model results to validate their accuracy.It was observed that the maximum relative error of β was <4.26% and the error of ∆T was <12.45% between CFD results and surrogate models.This further indicated that optimization results based on the PRSM model achieved a high level of accuracy.To quantify the backflow phenomenon, a line (Figure 11) parallel to the SEN axis was drawn through the center of the port to analyze velocity magnitudes (Figure 12).Overall, as molten steel rebounded after impinging at the bottom of the SEN, the velocity gradually decreased from the lower to upper sides, eventually leading to backflow.The effective area fraction of SEN increased by 9.20% post-optimization, indicating enhanced steel-flow delivery with reduced backflow compared with the initial design.Given this phenomenon, it was observed that jets exhibited a degree of dynamic distortion, thereby necessitating the further fluid dynamic analysis of the SEN [12,26].To quantify the backflow phenomenon, a line (Figure 11) parallel to the SEN axis was drawn through the center of the port to analyze velocity magnitudes (Figure 12).Overall, as molten steel rebounded after impinging at the bottom of the SEN, the velocity gradually decreased from the lower to upper sides, eventually leading to backflow.The effective area fraction of SEN increased by 9.20% post-optimization, indicating enhanced steel-flow delivery with reduced backflow compared with the initial design.Given this phenomenon, it was observed that jets exhibited a degree of dynamic distortion, thereby necessitating the further fluid dynamic analysis of the SEN [12,26].Given these findings and principles from boundary layer theory, it was observed that boundary layer separation occurred when molten steel reached the upper internal side of the port.This separation created low-pressure zones with significant velocity changes, resulting in high kinetic energy dissipation (Figures 13 and 14).The initial three-port SEN exhibited spatial inhomogeneity with a large low-pressure zone at the upper side of the Given these findings and principles from boundary layer theory, it was observed that boundary layer separation occurred when molten steel reached the upper internal side of the port.This separation created low-pressure zones with significant velocity changes, resulting in high kinetic energy dissipation (Figures 13 and 14).The initial three-port SEN exhibited spatial inhomogeneity with a large low-pressure zone at the upper side of the port (Figure 13).Conversely, the op2 design exhibited fewer low-pressure zones in the backflow area.Additionally, the exiting flow of the optimized SEN was homogeneously distributed with reduced energy dissipation (Figure 14).Given these findings and principles from boundary layer theory, it was observed that boundary layer separation occurred when molten steel reached the upper internal side of the port.This separation created low-pressure zones with significant velocity changes, resulting in high kinetic energy dissipation (Figures 13 and 14).The initial three-port SEN exhibited spatial inhomogeneity with a large low-pressure zone at the upper side of the port (Figure 13).Conversely, the op2 design exhibited fewer low-pressure zones in the backflow area.Additionally, the exiting flow of the optimized SEN was homogeneously distributed with reduced energy dissipation (Figure 14).Given these findings and principles from boundary layer theory, it was observed that boundary layer separation occurred when molten steel reached the upper internal side of the port.This separation created low-pressure zones with significant velocity changes, resulting in high kinetic energy dissipation (Figures 13 and 14).The initial three-port SEN exhibited spatial inhomogeneity with a large low-pressure zone at the upper side of the port (Figure 13).Conversely, the op2 design exhibited fewer low-pressure zones in the backflow area.Additionally, the exiting flow of the optimized SEN was homogeneously distributed with reduced energy dissipation (Figure 14).The op2 design effectively mitigated backflow at the port outlet, optimizing the use of the entire outlet area.It exhibited superior fluid dynamics performance, ensuring more uniform and stable molten-steel flow in the mold, thereby enhancing the continuous-casting processes.Thus, optimizing the SEN structure through a multi-objective approach proved beneficial for enhancing performance.

Conclusions
Upon conducting numerical simulations and multi-objective optimization modeling, this study investigated the influence of different SENs on fluid flow, heat transfer, and solidification in an ultra-large beam-blank continuous-casting mold.A new SEN structure design was developed based on these findings.The analysis of simulation results led to the following conclusions: (1) The maximum jet velocity of the three-port SEN was reduced to 1.1 m/s, while the two-port SEN was 1.3 m/s.Additionally, owing to the longer distance between the three-port SEN and the mold wall, the change rate in free-surface velocity was significantly smaller compared with the two-port SEN, and the velocity distribution was more uniform.
(2) The three-port SEN ensured a more consistent mold-outlet thickness distribution compared with the two-port SEN and exhibited a more acceptable surface temperature trend along the casting direction.This design reduces the risk of breakout due to skin rupture at the fillet.
(3) The op2 design of the three-port SEN demonstrated an enhanced effective area fraction at the port compared with the pre-optimized SEN, thereby reducing boundary layer separation and mitigating low-pressure areas and high dissipation to a considerable extent.Furthermore, the optimized SEN exhibited enhanced fluid dynamics and achieved a more stable flow pattern.
This study underscored the effectiveness of the new SEN design in enhancing casting quality and operational stability in ultra-large beam-blank continuous-casting processes.Please note that surrogate models are not limited to PRSM; exploring alternative approaches is recommended.For more information on methods that can be employed, please refer to citation [38].

2. 3 .
Governing Equations of the CFD Models 2.3.1.Fluid Flow Model Two types of fluid flow were considered within the cast strand: interdendritic flow in the mushy zone and turbulent flow in the liquid core [18].The standard k − ε model was employed to simulate the turbulence effects of fluid flow inside the cast strand [19].The following are expressions of the governing equations for fluid flow.

Figure 2 .
Figure 2. Schematic of (a) geometry model and (b) mesh structures.

Figure 2 .
Figure 2. Schematic of (a) geometry model and (b) mesh structures.

S 2 mC 1 J
are the areas of the mold outlet and SEN inlet ( ), respectively;SEN Dis the hydraulic diameter of the SEN ( m ); w ρ is the cooling water density ( is the specific heat of the cooling water ( kg K ); w Q is the cooling water flow (

−
m s ); and water T Δ is the temperature difference in the cooling water between the inlet and outlet of the copper tube ( K ).eff S represents the effective contact area between the molten steel and mold walls ( 2 m ).The flow and heat transfer calculations incorporate the casting speed, , and temperature difference between the inlet and outlet of the copper tube,

3 .
Comparison of Fluid Flow, Heat Transfer, and Solidification for Two-Port SEN and Three-Port SEN Figure4a, depicts the flow pattern of the two-port SEN near the flange.Upon impact with the initially solidified shell of the flange at a high velocity of 1.3 m/s , the SEN jet flow divided into two longitudinal circulations.A strong upper circulation centered at 0.04 m below the meniscus close to the free surface induced significant surface fluctuations.The circular flow below the SEN, centered at 0.3 m below the meniscus and close to the flange, did not facilitate the effective removal of non-metallic inclusions near the mold center[24].Contrastingly, Figure4billustrates the injection of molten steel at lower velocities toward the mold corners with the three-port SEN, owing to the longer distance between the nozzle outlets and mold corners.The centers of upper and lower circulations, located at 0.06 m and 0.3 m below the meniscus, respectively, exhibited significantly weaker upper circulations and a shift in the impact zones of circular flow toward the mold center, thereby facilitating meniscus stability and inclusion removal.

Figure 4 .
Figure 4. Comparison of velocity vectors for (a) two-port SEN and (b) three-port SEN.

Figure 4 .
Figure 4. Comparison of velocity vectors for (a) two-port SEN and (b) three-port SEN.

Figure 5 .
Figure 5.Comparison of maximum free-surface velocity for two-port SEN and three-port SEN.

Figure 5 .
Figure 5.Comparison of maximum free-surface velocity for two-port SEN and three-port SEN.

Figure 6 .
Figure 6.Comparison of surface temperature variation along the casting direction for two-port SEN (red lines) and three-port SEN (blue lines).

Figure 7
Figure 7 shows the liquid fraction at different cross-sections under both types of SEN.In this study, a liquid fraction of 0.33 defined the solidification front.Under the two-port SEN, the thickness of the solidified shell at the mold-outlet (z = 0.78 m ) web center, flange center, and fillet measured 11.75 m m , 12.05 m m , and 8.607 m m , respectively.For the three-port SEN, these values were 11.6 m m , 11.48 m m , and 10.238 m m , indicating a maximum thickness change reduction of 2.081 m m compared with the two-port SEN.This difference was attributable to the fact that ports of the three-port SEN were farther from the mold wall, thereby resulting in lower molten-steel impact velocities compared with the two-port SEN.Consequently, the homogeneity of the solidified shell was less affected by shock pressure from jet flows against the mold wall.Thus, the solidification behavior of the three-port SEN was more acceptable within the mold, potentially reducing the risk of liquid breakout[25].

Figure 6 . 21 Figure 7 .
Figure 6.Comparison of surface temperature variation along the casting direction for two-port SEN (red lines) and three-port SEN (blue lines).Materials 2024, 17, x FOR PEER REVIEW 11 of 21

Figure 7 .
Figure 7.Comparison of the liquid fraction at various cross-sections for (a) two-port SEN and (b) three-port SEN.

Figure 9 .
Figure 9. Relative error of the CFD results and surrogate models at five verification points.Figure 9. Relative error of the CFD results and surrogate models at five verification points.

Figure 9 .
Figure 9. Relative error of the CFD results and surrogate models at five verification points.Figure 9. Relative error of the CFD results and surrogate models at five verification points.

Figure 10 illustrates
Figure 10 illustrates the Pareto fronts for optimizing the three-port SEN structure using the NSGA-II algorithm.Notably, two goals-β and T Δ -contradicted one another, as an increase in β invariably increased T Δ .Therefore, the optimal design requires bal- ancing practical engineering needs.Upon focusing solely on β , the optimization of point 1 (op1)'s design features a bore diameter of 45 m m , an outlet angle of 11  , an immersion depth of 110 m m , and port horizontal angles of 114  .The actual value β , de- rived using the CFD model, was 0.9162, and was slightly higher compared with the surrogate model prediction of 0.9106.Conversely, given the single objective T Δ , the op3 de- sign was obtained with a bore diameter of 46 m m , an outlet angle of 9  , an immersion depth of 111 m m , and port horizontal angles of 112  .In this design, the values of β

2 u= , and k i F
is the objective function i at the optimization point k.The optimized SEN featured a bore diameter of 45 m m , an outlet angle of 10  , an immersion depth of 110 m m , and port horizontal angles of 114  .In comparison with the initial three-port SEN, β increased by 9.20% and T Δ decreased by 45.12%, and m S and s V , met metallurgical standards.Additionally, Table

Figure 11
Figure 11 illustrates the velocity profiles at the outlet port of the initial three-port SEN and the op2 design, examining the recirculation effect.It was observed that fluid recirculated backward at the upper edge of the port, thereby indicating an incomplete utilization of the entire outlet area.

Figure 11
Figure 11 illustrates the velocity profiles at the outlet port of the initial three-port SEN and the op2 design, examining the recirculation effect.It was observed that fluid recirculated backward at the upper edge of the port, thereby indicating an incomplete utilization of the entire outlet area.

Figure 11 .
Figure 11.Velocity profiles at the central symmetrical plane of (a) initial three-port SEN and (b) op2 design.

Figure 11 .
Figure 11.Velocity profiles at the central symmetrical plane of (a) initial three-port SEN and (b) op2 design.

Materials 2024 , 21 Figure 12 .
Figure 12.Velocity magnitude along a line from upper to lower sides of the SEM port center.

Figure 12 .
Figure 12.Velocity magnitude along a line from upper to lower sides of the SEM port center.

Figure 12 .
Figure 12.Velocity magnitude along a line from upper to lower sides of the SEM port center.

Figure 13 .
Figure 13.Pressure contours at the center symmetrical plane of (a) initial three-port SEN and (b) op2 design.

Figure 13 .
Figure 13.Pressure contours at the center symmetrical plane of (a) initial three-port SEN and (b) op2 design.

Figure 12 .
Figure 12.Velocity magnitude along a line from upper to lower sides of the SEM port center.

Figure 13 .
Figure 13.Pressure contours at the center symmetrical plane of (a) initial three-port SEN and (b) op2 design.

Figure 14 .
Figure 14.Kinetic energy dissipation contours at the center symmetrical plane of (a) the initial three-port SEN and (b) op2 design.

Table 1 .
Physical parameters of the Q235B.

Table 2 .
Parameters for the NSGA-II algorithm.

Table 3 .
Accuracies of the PRSM surrogate models.

Table 3 .
Accuracies of the PRSM surrogate models.

Table 4 .
Comparison of optimal and initial designs for the three-port SEN.

Table 4 .
Comparison of optimal and initial designs for the three-port SEN.