Simulation of Steel Casting Thermal Behavior

A deep dive into the governing physics, numerical methods, and real-world engineering tools that model heat transfer, solidification, and micro structure formation in continuous steel casting — from first principles to factory-floor automation.

Simulation of Steel Casting Thermal Behavior
Continuous Casting • Stefan Problem • Solidification Physics

Start with the Single
Governing Idea:
Energy + Latent Heat

Every continuous casting simulation ultimately begins with one governing principle: energy conservation. Heat flows through the material according to classical conduction laws, but at the solid-liquid interface an additional phenomenon emerges. As liquid metal transforms into solid, enormous quantities of latent heat are released. The interaction between heat transport and phase transformation creates a moving boundary problem known as the Stefan problem, which forms the foundation of modern solidification modeling.

Q
The Fundamental Casting Equation

Heat Flows.
Metal Solidifies.
Interfaces Move.

Unlike conventional heat transfer problems, continuous casting requires simultaneous prediction of both temperature evolution and the constantly moving solid-liquid interface responsible for shell growth.

Four Core Physical Concepts

Fourier Conduction
Stefan Condition
Latent Heat Release
Free Boundary Motion
1
Heat Transport Physics

Fourier-Type Conduction

Within both the fully liquid steel and the growing solid shell, thermal energy transport is governed by Fourier heat conduction. This equation describes how temperature changes over space and time as energy flows from hotter regions toward cooler regions.

The complexity arises because liquid steel and solid steel possess different thermal properties. Thermal conductivity, density, and specific heat vary across the phase boundary, requiring simulation algorithms to update material properties continuously as solidification progresses.

Properties That Change Across The Interface

Thermal Conductivity
Density
Specific Heat
2
Interface Energy Balance

The Stefan Condition

At the moving solid-liquid boundary, a special energy balance known as the Stefan condition must be satisfied. This relationship links interface velocity directly to the difference in heat flux arriving from the liquid side and leaving through the solid side.

The resulting energy imbalance cannot disappear. Instead, it is consumed by the phase transformation process itself, driving movement of the solidification front and controlling shell growth throughout continuous casting.

Stefan Energy Balance

Heat Flux In
−
Heat Flux Out
→
Interface Movement
3
Dominant Energy Source

Why Latent Heat Dominates

During solidification, steel releases enormous quantities of latent heat. For typical steel grades, latent heat is approximately 270 kJ/kg, making it one of the largest energy contributions in the entire casting system.

This energy release strongly influences shell growth, temperature gradients, mold heat flux, and cooling behavior. Any simulation that neglects latent heat or models it inaccurately will significantly underestimate solidification rates and shell thickness development.

270 kJ/kg

Typical Steel Latent Heat

The latent heat released during phase transformation can equal or exceed the sensible heat contribution over large temperature intervals, making it a dominant factor in continuous casting simulations.

Latent Heat Controls

Shell Growth
Cooling Rate
Thermal Gradient
Solidification Time
4
Mathematical Complexity

Interface as a Free Boundary

Classical heat-transfer problems typically involve fixed boundaries whose positions are known in advance. Continuous casting is fundamentally different. The location of the solid-liquid interface is unknown and must be determined simultaneously with the temperature field.

This moving interface is described as a free boundary. Every change in temperature influences interface position, while every movement of the interface alters the temperature field. The two quantities remain tightly coupled throughout the entire solution process.

Coupled Solution Process

Temperature Field
↔
Interface Position
↔
Shell Thickness
One Equation

Governs The Entire Process

Continuous casting simulation begins with energy conservation, but evolves into a coupled solidification problem where heat flow, latent heat release, and interface motion must

Numerical Challenge

Why “Solidification” Breaks Simple Simulations

The elegance of the Stefan formulation conceals a profound numerical challenge: once the interface position becomes an unknown, the problem becomes inherently nonlinear.

∞
The Root of Nonlinearity

The Circular Dependency

In a standard conduction problem, the computational mesh is fixed and temperatures are solved at known node locations. In a moving-boundary problem, the interface position changes with every time step — and its position depends on the very temperature field you are trying to solve.

Temperature
Thermal Field
Interface
Moving Boundary
Time Step
New Position

This circular dependency makes explicit time-stepping schemes unstable, while implicit schemes require careful treatment to avoid oscillation or non-convergence. Mesh-tracking approaches can be accurate but become computationally expensive and prone to mesh distortion over large displacements.

Numerical Remedies

Two Ways to Handle the Phase Change

Method 01

The Enthalpy Method

H

Instead of tracking the interface explicitly, the governing equation is reformulated in terms of total enthalpy, including both sensible and latent heat contributions.

Key Transformation
Sharp phase boundary → finite-temperature mushy zone → fixed-domain nonlinear PDE

The approach allows standard solvers to handle the problem, although some interface sharpness is sacrificed.

Method 02

Effective Specific Heat

Cp*

An alternative, often used in commercial codes, artificially increases the specific heat over the solidification temperature range to absorb the latent heat contribution.

Numerical Sensitivity

A poorly selected temperature interval or overly coarse time step can introduce significant error and may cause latent heat release to be missed between consecutive iterations.

Industrial Codes

Why Generic Thermal Solvers Are Not Enough

These numerical challenges explain why industrial casting simulators require significantly more sophisticated code architecture than generic thermal solvers.

GRID RESOLUTION
Must resolve phase-change behavior accurately.
TIME STEPPING
Must capture latent heat release without skipping critical changes.
CONVERGENCE
Criteria must account for nonlinear phase-change physics.
PROPERTY DATA
Material properties must be interpolated appropriately across the phase range.
!

Plausible Does Not Always Mean Accurate

Poorly designed codes may produce visually plausible temperature contours while fundamentally misrepresenting shell growth rates — with potentially dangerous consequences for process-control decisions.

Physical Reality

The Mushy Zone Reality

Real steels do not solidify at a single temperature. Solid and liquid coexist across a temperature range as a two-phase mixture.

Liquid Mushy Zone Solid
Solid + Liquid coexist
Viscosity
Influenced by the fraction solid.
Permeability
Changes as the two-phase structure develops.
Thermal Conductivity
Varies throughout the mushy region.
Mechanical Strength
Evolves with the increasing fraction solid.
Thermodynamic Coupling

Connecting Temperature to Fraction Solid

Capturing the mushy zone correctly requires coupling the thermal solver to thermodynamic databases such as IDS or Thermo-Calc, which compute equilibrium or non-equilibrium fraction solid as a function of temperature and alloy composition.

Temperature
→
Alloy Composition
→
Fraction Solid
The Modeling Shift

From Sharp Boundary to Physical Reality

SHARP STEFAN
Explicit moving interface
MUSHY ZONE
Solid and liquid coexist
ENTHALPY
Fixed-domain nonlinear formulation

Accurate solidification simulation requires numerical methods designed around the nonlinear physics of phase change — not simply a conventional heat-conduction solver.

Mold/Shell Interaction

The Thermal Detail That Damages Parts

Taper: Too Much vs. Too Little

Mold taper is designed to follow shell contraction. Excessive taper squeezes the shell prematurely, causing cracks and buckling. Insufficient taper opens air gaps, reducing heat transfer and risking catastrophic breakouts as liquid steel bursts through weak shells.

Interface Physics: Flux, Distortion, Funnel Effects

Optimal taper must account for copper mold distortion, flux layer thickness and viscosity, and funnel geometry in thin-slab casting. Funnel molds demand steeper taper adjustments due to longer contact lengths and complex geometry.

Casting Speed Reshapes Everything

Higher casting speeds shift heat-flux profiles, reduce shell thickness at mold exit, and demand speed-dependent taper optimization. Even deviations of tenths of a millimeter outside the optimized taper window significantly increase defect rates.

The shell/mold thermal interaction is a micron-scale phenomenon with meter-scale consequences: a flux gap of 0.1–0.3 mm can halve heat flux and double shell temperature, directly triggering cracks and breakouts.

Online Casting Control • CastManager • IDS • Tempsimu

Online Era: Fast Coupled
Heat Transfer + Solidification
+ Microstructure

The ultimate objective of casting simulation is not simply offline analysis but active process control. Advances in computational speed, validated physical models, and industrial-scale data integration have enabled simulation systems capable of supporting real-time operational decisions. Modern online-capable environments combine thermal modeling, solidification prediction, microstructure evolution, and cooling optimization into a unified framework that continuously guides production toward defined quality targets.

AI
Online Process Optimization

Predict.
Optimize.
Control.

Modern casting software ecosystems continuously transform thermal data into solidification predictions, quality metrics, and process-control recommendations before defects have an opportunity to form.

Three Core Digital Platforms

CastManager
IDS
Tempsimu
1
Real-Time Thermal Modeling

CastManager: Heat Transfer Intelligence

CastManager is designed specifically for continuous casting environments where calculations must be completed fast enough to support operational decisions. The software computes thermal behavior from mold entry through secondary cooling and up to the final solidification point along the strand.

Its primary challenge is balancing computational speed with physical accuracy. Advanced numerical optimization techniques allow the system to produce accurate shell thickness, mold heat-flux, and temperature predictions while remaining suitable for real-time control environments.

CastManager Computes

Strand Temperature
Mold Heat Flux
Shell Thickness
Thermal Fields
2
Metallurgical Prediction

IDS: Solidification & Microstructure

IDS (Interdendritic Solidification) converts thermal predictions into metallurgical outcomes. Using alloy chemistry and local cooling conditions, it calculates fraction-solid evolution, solidification paths, partition coefficients, and equilibrium or non-equilibrium transformation behavior.

The software transforms temperature histories into practical quality indicators, allowing engineers to predict casting performance before defects become observable in production.

IDS Predicts Quality Risks

Segregation
Porosity
Hot Cracking
Grain Structure
Thermal ↔ Metallurgical

Two-Way Coupling

Cooling conditions influence solidification behavior, while evolving solid fractions modify local thermal properties. IDS and CastManager continuously exchange information through this coupled relationship.

3
Secondary Cooling Engineering

Tempsimu: Cooling Optimization

Tempsimu focuses on the secondary cooling region where most solidification is completed. Before engineers implement changes on the caster, the software evaluates proposed modifications such as spray-zone redesign, nozzle adjustments, altered water distributions, roll configuration updates, or revised cooling-zone boundaries.

By evaluating these changes virtually, engineers can identify undesirable thermal consequences before operational deployment, reducing production risk and protecting product quality.

Tempsimu Evaluates

Nozzle Layout
Roll Configuration
Water Distribution
Cooling Zones
Integrated Optimization Loop

Three Tools. One Decision Engine.

CastManager supplies thermal predictions, IDS converts those predictions into microstructural and quality metrics, and Tempsimu validates cooling strategies. Together they create a continuously improving optimization loop.

Online Control Architecture

CastManager
→
IDS
→
Tempsimu
→
Control Action

Continuous Casting

From Equations to Reality: Continuous Casting Is Not Just Heat

While the Stefan problem captures essential solidification physics, a real continuous casting machine demands far more. High-fidelity simulation couples turbulent fluid flow, thermal transport, and solidification into a computationally demanding multi-physics problem.

The Real Casting Problem

Three Physics, One Process

Continuous casting cannot be represented accurately by thermal conduction alone. The mold region is governed by interacting flow, heat transfer, and phase-change phenomena.

FLOW
Turbulent Motion
Jet impingement & recirculation
HEAT
Thermal Transport
Heat flux & temperature fields
SOLID
Shell Formation
Growth & structural stability
01 — Fluid Dynamics

Full 3D Navier–Stokes Coupling

∇·u

Liquid steel flowing through the submerged entry nozzle (SEN) and into the mold pool is turbulent, highly inertial, and thermally stratified. High-fidelity studies solve the full three-dimensional, time-dependent Navier–Stokes equations coupled to the energy equation.

What the Flow Model Reveals
Jet Impingement Recirculation Zones Surface Velocities Inclusion Transport

Turbulence models such as k-ε, LES, or hybrid approaches capture the flow structures that directly influence shell formation and inclusion transport.

Q
02 — Thermal Boundary

Mold Heat Flux: The Critical Boundary

At the point where the liquid steel jet strikes the solidifying shell — the jet impingement zone — local heat fluxes reach extraordinary magnitudes. These extreme values drive rapid initial shell formation and influence whether the nascent shell can withstand ferrostatic pressure as it exits the mold.

Near-Field Peak Heat Flux
800–1,300 kW/m²

Approximate peak range under near-field conditions, depending on casting conditions and spatial grid resolution.

03 — Computational Scale

Grid Resolution Demands

Steep thermal and velocity gradients near the jet impingement zone and solidification front demand very fine computational grids. Industrial-scale molds can therefore push simulations into high-performance computing territory.

Coarser regions Refined critical zones
600–1600
mm mold width range
700–1200
mm mold height range

Coarse grids systematically under-resolve peak heat flux values, potentially leading to non-conservative shell-thickness predictions.

The Critical Insight

Why Flow Cannot Be Ignored

A purely conductive thermal model of the mold pool dramatically misrepresents the temperature distribution. Turbulent convection by the steel jet homogenizes much of the liquid pool, elevating temperatures in recirculation zones and maintaining superheat deeper into the strand than conduction alone would predict.

Conduction Only

Simplified Thermal Field

Omits turbulent convection and can over-predict shell thickness while underestimating superheat penetration.

Flow + Thermal

Coupled Reality

Captures turbulent mixing, recirculation, superheat transport, and their influence on shell formation.

Turbulent Jet → Recirculation → Superheat Transport → Shell Formation
Practical Implications

Simulation Becomes a Process-Design Tool

Coupled flow-thermal simulations reveal how SEN geometry, mold electromagnetic stirring configurations, and casting speed influence the process. This enables virtual prototyping of new nozzle geometries before costly physical trials.

JET ANGLE
Alters impingement patterns
PORT AREA
Influences flow distribution
SUBMERGENCE
Changes jet interaction
Result
Virtual prototyping → Faster development → Improved first-cast success
Final Insight

Continuous Casting Is a Coupled Physics Problem

Accurate simulation must connect turbulent flow, heat transfer, and solidification. Only by resolving their interaction can engineers predict shell formation, superheat penetration, heat-flux distribution, and the process response of a real continuous casting machine.

What's Your Reaction?

like

dislike

love

funny

angry

sad

wow